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Preface 


Although SPICE has been with us since 1971, use of the program dur¬ 
ing the past 10 years has exploded because of the personal computer. 
IBM’s introduction of the desktop personal computer and the 
widespread availability of SPICE, both in the academic and commer¬ 
cial form, have made analog circuit simulation available to almost 
every circuit designer. 

As a designer, the author struggled for many years with SPICE 
because, although the program holds an unlimited potential for aiding 
in circuit design, practical use of the program is limited by nonconver¬ 
gence failures, timestep control errors, and numeric integration fail¬ 
ures. These types of problems either stop the simulator or introduce 
error in the simulation output. Because of these failures, simulation 
results are either difficult to obtain or contain numerous inaccuracies. 

Although these types of errors occur almost regularly, very little has 
been written about how to correct these simulation ills. The motivation 
for this text came from the lack of information available to designers 
on how to overcome these obstacles of simulation. Obtaining accurate 
high-quality simulation results requires learning how to overcome 
these simulation failures. 

Someday, mathematicians and programmers will be able to produce 
a circuit simulator which always converges on the correct solution and 
always produces accurate results. Until then, as SPICE users, we must 
learn to overcome these limitations. 


Book Organization 

This book is written as both a tutorial and as a reference. 

As a tutorial, the text includes a bound copy of the RSPICE™ circuit 
simulator and RGRAPH™ graphical postprocessor. Chapters 3, 4, 5, 
and 6 contain example circuits which should be simulated with 
RSPICE/RGRAPH. These circuits are integrated into the text and 
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demonstrate the common failure mechanisms of SPICE as well as the 
actions required to correct the failure. The disk which accompanies the 
book contains all of the circuit files, organized by chapter, mentioned in 
the text. 

As a reference, the book includes systematic, step-by-step proce¬ 
dures to correct common simulator failures. The book also includes 
multiple examples of simulation failures. These examples can be used 
to identify similar failures from the reader's simulations. 


RSPICE/RGRAPH Disk 

The accompanying RSPICE/RGRAPH disk contains the RSPICE cir¬ 
cuit simulator and the RGRAPH graphical postprocessor. The disk also 
contains the circuit examples used throughout this text. An installa¬ 
tion program automates the task of creating the proper directory struc¬ 
ture, copying the files, and assigning the proper directory pointers. 


Requirements 

To install the programs and files, you will need approximately 2 MB of 
free space on your hard disk. The install program creates the following 
directories and subdirectories. 

\RCGV33 

\RCGV33\RSPICE 

\RCGV33\CH2 

\RCGV33\CH3 

\RCGV33\CH4 

\RCGV33\CH5 

\RCGV33\CH6 

\RCGV33\DEM0 

The RSPICE and RGRAPH files are then copied to the appropriate 
directories. 


Installation 

We'll assume you're installing the software from a floppy disk. To 
install the programs, simply type: 

X:instal1 X: Y: 

Where X is the floppy disk drive designator and Y is the hard disk drive 
designator. The installation program will complete the remaining 
installation chores. 
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Summary of Chapter Content 

Chapter 1 is a brief overview of circuit simulation. Included in this 
chapter is a historical review of the development of CANCER, SPICE1, 
SPICE2, and SPICE3. The chapter continues with how computer sim¬ 
ulation can be an aid in the design process and what is needed to sim¬ 
ulate circuits productively and effectively. 

Chapter 2 looks at how SPICE works, and shows how SPICE con¬ 
structs the system equations from the elements of the input file, how 
SPICE iterates to a solution, how SPICE computes the DC bias, DC 
sweep, AC frequency sweep, and transient analysis solutions. Under¬ 
standing the fundamental algorithms of SPICE is crucial to under¬ 
standing how to correct the common simulator failure mechanisms. 

Chapter 3 examines the problem of nonconvergence. Nonconver¬ 
gence is one of the most common and most frustrating problems facing 
simulation users. But most nonconvergence problems can be overcome 
simply by using the options and controls found within SPICE. Chapter 
3 looks at the causes of nonconvergence and presents systematic, step- 
by-step procedures which will eliminate (most) nonconvergence simu¬ 
lation failures. 

Chapter 4 is the first of two chapters which focus on improving the 
accuracy of the transient analysis. In Chapter 4, the three numeric 
integration routines of SPICE are examined. Numeric integration fail¬ 
ures introduce error into the simulation result, and SPICE does not 
have a routine which warns users of integration failure. As simulation 
users, we must learn to detect integration failure from the simulation 
result and take the appropriate corrective action. 

Chapter 5 finishes the discussion on transient analysis by looking at 
the timestep control algorithms in SPICE. Like numeric integration, 
SPICE cannot detect timestep control failures. Users must learn to 
look for timestep failures in the simulation output. Chapter 5 exam¬ 
ines different timestep control failures and illustrates the proper cor¬ 
rective action. 

Chapter 6 examines the 34 options within SPICE. Each of the 
options are presented with a detailed explanation and suggested set¬ 
tings. 

Appendix A compares the default option settings of several popular 
vendor-offered circuit simulators. 

Throughout this text the program SPICE is referenced many times. 
For this book, SPICE refers specifically to Berkeley SPICE2G.6. 


Ron M. Kielkowski 
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Chapter 



What Is SPICE? 


SPICE, and specifically SPICE2G.6, is an engineering design tool 
which falls into the category of general purpose analog circuit simula¬ 
tors. Just some of the tools which fall into this same category include 
SPICE3, PSPICE, HSPICE, IS_SPICE, MICRO-CAP IV, RSPICE, 
SABER, and many others. 

Computer simulation can be a powerful supplement to traditional 
design techniques. In most design strategies, simulation can be an 
aid in the initial design development, during the breadboarding 
phase, and during debugging and diagnostic phases. For some cir¬ 
cuits, initial design theories must be tested before circuit design 
begins. With simulation, circuit blocks may be represented as behav¬ 
ioral elements and simulated in a functional form. Behavioral ele¬ 
ments allow designers to test circuit theory without the time involved 
in developing transistor and component-level descriptions of each cir¬ 
cuit function. 

For many circuits, breadboarding is impossible because of excessive 
circuit complexity, layout-specific parasitic effects, or, as in the case of 
integrated circuits, both effects. For these types of circuits, simulation 
*nay be the only avenue to investigate the circuit behavior before build- 
mg a working prototype. 

For most circuits, component-value variation will have a direct effect 
on circuit performance and product yield. With simulation, designers 
can effectively predict the performance of a circuit as one or more cir¬ 
cuit variables are changed. 

For all of these reasons, computer simulation and SPICE are playing 
an increasingly important role in electronic circuit design. 
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matically vary circuit components as multiple simulations are per¬ 
formed. Monte Carlo and worst-case analyses enhance a user’s ability 
to measure circuit performance over a number of changing circuit 
parameters. 


Evaluating circuit vulnerability 

Computer simulation gives the designer a chance to test a circuit 
for vulnerability prior to being assembled. Sneak-path conditions, 
signal race conditions, and power dissipation/power overload condi¬ 
tions can all be quickly and easily investigated through the use of 
simulation. While many of these conditions could be examined on a 
breadboard, when circuit simulation models are forced beyond their 
performance limits, they don’t smoke, smolder, smell, fuse, or explode 
like their breadboard equivalents. Simulation allows a designer to 
examine the circuit without the risk of damage to the circuit or the 
designer! 

Failure analysis/diagnosis/prognosis 

As a failure-analysis tool, circuit simulation has limits. The most obvi¬ 
ous limit is the lack of ability to predict layout-dependent parasitic 
behaviors. SPICE simulates components as if they were ideal, non¬ 
interacting elements. For example, a common problem with CMOS 
(complementary metal-oxide-silicon) circuits is the phenomenon 
known as latch-up. Latch-up occurs when an N-channel and P-channel 
transistor pair are fabricated in close proximity. Because of the p- 
and n-type silicon layers in the transistor pair, a pair of parasitic bipo¬ 
lar transistors is also formed. Latch-up occurs when the parasitic 
bipolar transistors begin to conduct current. If an N-channel and 
P-channel transistor pair are entered in SPICE, the latch-up effect will 
not be observed unless the parasitic bipolar transistors are included in 
the netlist. Figure 1.3 illustrates the parasitic bipolar transistors 
formed in P-well CMOS circuits. Trying to use SPICE to predict this 
type of failure is futile. 


Substrate P-wel! 

Source Gate Drain contact Drain Gate Source contact 

_xi_/ ^_ n _rx_ 

N+ J [ N+ J 

P-well 

_ y 

Parasitic PNP 

N-Substrate 

^Hjure 1.3 Parasitic bipolar transistors formed in a P-well CMOS inverter. 









2 Chapter One 


What Is a SPICE-Compatible Simulator? 

Almost all circuit simulators read a circuit file which describes the 
types of elements in the circuit and the type of analysis to be per¬ 
formed. After processing the circuit file, the simulator performs the 
desired analysis and generates the output in either tabular or graphi¬ 
cal form. 

Of all the analog circuit simulation tools available, the overwhelm¬ 
ing majority of them are SPICE-like or SPICE-compatible. SPICE-like 
means a simulator is capable of producing an analysis result similar 

the SPICE result for a given circuit, although the SPICE-like simu¬ 
lator may not be able to read a standard SPICE circuit file. SPICE- 
compatible means a simulator will read a SPICE circuit file, perform 
the desired analysis and produce the output result in standard 
SPICE2G.6 form. Table 1.1 illustrates the SPICE compatibility of sev¬ 
eral vendor-offered simulators. 

The Birth of SPICE 

Ironically, the man most responsible for the development of SPICE 
wasn’t interested in circuit design or circuit simulation. In 1968, Ron 
Rohrer, a junior faculty member at the University of California, Berke¬ 
ley, and an authority on circuit optimization , was assigned to teach a 
class in network synthesis. But Rohrer was convinced that synthesis 
was of little practical use. Rohrer was much more interested in circuit 
optimization and optimization techniques. Circuit optimization involves 
performing multiple simulations on a circuit to study the change in out¬ 
put with a change in one or more circuit variables. Up to this time, 
Rohrer’s work with optimization had been stifled by the lack of a fast, 
efficient circuit simulator. Because of the unavailability of simulation 
tools which were appropriate for use in circuit optimization, Rohrer 
decided to replace the course on circuit synthesis with a class on circuit 
simulation. Ultimately, Rohrer was hoping to develop a new circuit sim¬ 
ulator for his work in optimization. 1 


TABLE 1.1 Vendor-Offered Simulator Compatibility 


SPICE-compatible 

SPICE-like 

Hspice 

Saber 

Pspice 

Micro-Cap II, III 

IS_Spice 


Rspice 


Micro-Cap IV 


Accusim 


Analog Workbench 
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cancer 

In the class on circuit simulation, Rohrer and a dozen students assem¬ 
bled a nonlinear circuit simulator which was to become the foundation 
for SPICE. This first simulator was christened CANCER 2 (Computer 
Analysis of Nonlinear Circuits Excluding Radiation) by a student 
named Larry Nagel. 

CANCER was capable of DC operating point, DC sweep, transient 
sweep, and AC frequency sweep analyses, the same analysis types 
available today in SPICE. CANCER simulations might include resis¬ 
tors, capacitors, inductors, and two types of nonlinear devices, junction 
diodes and bipolar transistors. Diodes were modeled with the Shockley 
diode equations, and bipolar transistors were modeled with the 
Ebbers-Moll transistor equations. While CANCER was an outstanding 
tool for its day, the program’s “lifetime” was limited because the solu¬ 
tion routines could handle no more than 400 components and/or 100 
circuit nodes. 

CANCER wasn’t the only simulator available in 1970. 3 Other pro¬ 
grams of the time included IBM’s ECAP and ECAP-II. BIAS 4 and 
SLIC 5 (Simulator for Linear Integrated Circuits) originated from 
early work on CANCER. Autonetics, a division of Rockwell, developed 
TRAC 6 (Transient Radiation Analysis by Computer Program) in 
1968. The TRAC program was the basis for Motorola’s TIME and 
MTIME 7 programs. TRAC was also the foundation for Berkeley’s 
SINC program. 

SPICE1 

During the early 1970s, Nagel continued improving the CANCER pro¬ 
gram. In 1971, this improved version, named SPICE 1 (Simulation Pro¬ 
gram with Integrated Circuit Emphasis ), was released into the public 
domain. Because Berkeley distributed the program with almost no 
charge, SPICE 1 quickly became an industry standard simulation tool. 
SPICE 1 offered several improvements over CANCER. The bipolar 
transistor model was changed to the Gummel-Poon 8 model. JFET and 
MOSFET devices were added with the Shichman-Hodges 9 model. 
SPICE 1 also offered a new approach to modeling known as macromod¬ 
eling. With macromodels, engineers could describe portions of a circuit 
m the form of relocatable circuit templates (subcircuits). 

During this time, the rapid development of the integrated circuit 
mdustry fueled the work on SPICE. In many ways, integrated circuits 
were very different from board-level circuits. Many IC problems could 
not be examined with traditional design techniques. Because of this, 
computer simulation of integrated circuits proved to be an invaluable 
design tool. 
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Because of the push from the integrated circuits industry, however, 
many of the algorithms in SPICE were optimized for these types of cir¬ 
cuits. This is a limitation which has lived on to this day. While SPICE 
can be used to simulate almost any type of electronic circuit, board- 
level circuits and circuits using inductors or transformers require 
extra care because many of the SPICE algorithms may not be well 
suited to these types of circuits. 10 

SPICE2 

The next major release of the program came in 1975 with the intro¬ 
duction of SPICE2. SPICE2 offered significant improvements over 
SPICE 1 or CANCER, including a new equation formulation for voltage- 
defined elements (inductors and voltage-controlled voltage and current 
sources). Both the accuracy and speed of transient analysis were 
improved with the development of two dynamic timestep control algo¬ 
rithms and a multi-order implicit integration scheme. SPICE2 device 
model improvements were added to keep pace with the changing 
device technologies. From 1975 through 1983, Berkeley continued 
improving and upgrading the SPICE2 program. In 1983, SPICE2 ver¬ 
sion G.6 (SPICE2G.6) was released to the public domain. SPICE2G.6 
was to be the last fortran version of SPICE released by the university. 
(SPICE2G.6 is still available from Berkeley.) 

SPICE3 

CANCER and both generations of SPICE were written in fortran 
source code. With the increasing use of Unix-based workstations, 
Berkeley made the decision to rewrite the SPICE2 program in C. The 
new C version of the program was known as SPICE3. 

SPICE3 was to be a superset of SPICE2, including all of the analysis 
types and device models of SPICE2 as well as new features such as 
improved device models, voltage- and current-controlled switches, 
pole-zero analysis, and a graphical postprocessor for viewing simula¬ 
tion results. Unfortunately, much of the enthusiasm (and research 
money) available in the early ’70s had disappeared, and the over¬ 
whelming task of converting some 22,000 lines of fortran code to C was 
left to a handful of students. Because of the enormity of the task, the 
first release of SPICE3A.1 contained dozens of bugs and coding errors. 
Worse yet, somewhere in the translation between fortran and C, much 
of the SPICE2G.6 functionality was lost. SPICE3 became a program 
which was not backward-compatible with SPICE2G.6. As of this writ¬ 
ing, SPICE3F.2 is the latest release of SPICES. Although Berkeley 
continues to improve the program, SPICE3 is not yet completely 
backward-compatible with SPICE2. 



What Is SPICE? 5 


Some of the features which have not been written into SPICE3 
include the ability to run multiple circuit netlists from a single file, 
polynomial-controlled dependent sources, polynomial capacitors and 
inductors, on-line resistor temperature coefficients, circuit topology 
*»rror checking, and temperature sweep analyses. 11 While work-arounds 
are available for most of these features, SPICE2 users may need to 
rewrite old circuit netlists to make them compatible with SPICE3; this 
includes rewriting old macromodels to make them SPICE3-compatible. 
Because of the hundreds of SPICE2 macromodels available today, and 
because of the large installed base of SPICE2-compatible simulators, 
the lack of backward compatibility is probably the overriding reason 
SPICE3 has not replaced SPICE2 as the industry-standard circuit 
simulator. 

Although it may not have all of the features of SPICE2, SPICE3 offers 
several technical advantages to SPICE2. SPICE3 is written in modular 
C code which is easier to modify than the fortran of SPICE2. SPICE3 
demonstrates superior convergence characteristics to the SPICE2 algo¬ 
rithms. In rewriting the SPICE3 device models, several SPICE2 
errors were discovered and corrected. As Berkeley continues adding 
enhancements and SPICE2 compatibility to the program, the day 
will come when SPICE3 will replace SPICE2. But since the vast major¬ 
ity of us are using a SPICE2-based simulator, this text will focus on 
SPICE2. 

From 1980 up through today, several vendor-offered versions of 
SPICE appeared. Some of the better-known simulators which got their 
start in this time frame include Meta-Software’s HSPICE, IntuSoft’s 
ISJ3PICE, Spectrum Software’s MICRO-CAP, and MicroSim’s 
PSPICE. All of these were developed from the original SPICE2 frame¬ 
work. Figure 1.1 illustrates the functional relationship between 
SPICE2G.6 and the vendor-offered simulators. While many other 
SPICE-based programs exist, these four represent the best-known 
simulators. Although a handful of vendors offer products based on 
SPICE3, the overwhelming majority of SPICE-like simulators are still 
based on SPICE2G.6. 

Why Simulate? 

Why should you simulate your circuits? What are the benefits of simu¬ 
lation? Too many times you hear the answer, “Circuit simulation 
replaces breadboarding.” Nothing could be further from the truth! 
Simulation does not replace the breadboard; simulation complements 
the breadboard . Many things can be learned from the breadboard that 
cannot be learned from simulation, and simulation may reveal many 
things not readily learned in the lab. 
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Figure 1.1 Graphical representation of SPICE-compatible simulators. 
(.Reprinted from Successfully Simulating Circuits With SPICE. Used with 

permission .) 


Although thousands of reasons can be given for simulating circuit 
designs, the majority of them fall into one of the following four cate¬ 
gories. 

Verifying design theories 

Simulation offers the ability to quickly test circuit design theories 
before a single wire is soldered or mask fabricated. Verifying the theory 
of a design can be done at different levels, behavioral models, macro¬ 
models, and circuit component (transistor) models. 

Behavioral models express circuit blocks in mathematical relation¬ 
ships. In behavioral form, simulating an entire system response 
becomes realizable. Behavioral blocks are easy to construct, and simu¬ 
late hundreds or thousands of times faster than a circuit at the compo¬ 
nent level. Behavioral blocks represent the highest level of the 
simulation hierarchy. 

Macromodels express circuit blocks in the form of simplified equiva¬ 
lent circuits. Macromodel circuit elements may include real circuit 
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compo nen te and ideal circuit components (dependent voltage and cur¬ 
rent sources). Macromodels simulate slower than behavioral models 
but much faster than transistor-level models. Because of their simpli¬ 
fied representation, macromodels do not give the fine detail or accu¬ 
racy of a transistor-level simulation. 

A transistor-level simulation is simply a simulation of the compo¬ 
nents in a circuit. This type of simulation offers the most accuracy of 
any analysis but often with a stiff run-time penalty. A transistor-level 
simulation represents the lowest level of the simulation hierarchy. 

In verifying circuit theories, designers often begin simulation with 
behavioral models. As the design progresses and more accuracy is 
required, macromodels replace the behavioral models. If even more 
accuracy is required, transistor-level models may replace the macro¬ 
models. In many cases, a simulation may be made up of behavioral 
models, macromodels, and transistor-level models, all in a single run. 
The ability to simulate at these different levels offers users a fast, effi¬ 
cient means to test the theory of a design. 


Circuit performance and yield analysis tool 

Simulation allows the designer to quickly test a circuit over a variety 
of conditions including temperature variations, element value varia¬ 
tions, and power supply variations. Once a circuit netlist has been 
established, circuit parameters (including temperature, element val¬ 
ues, or power supply levels) may be altered and resimulated. The abil¬ 
ity to alter circuit parameters offers designers a fast, efficient means of 
testing the circuit operation under a variety of operating conditions. 

In the real world, as opposed to the simulation world, resistors, 
capacitors, inductors, transistors, and other electrical component val¬ 
ues vary from part to part. Even under the tightest process control, 
most common electrical components have part-to-part variations of 
1 percent to 5 percent or more. Designers recognize this and try to 
develop designs which are insensitive to these variations. If the yield of 
a circuit depends on the component value variations, simulation may 
be useful in estimating circuit yields. 

Circuit simulation offers users a fast, efficient means to test a cir- 
^t s performance characteristics as one or more of the component val¬ 
ues change. By performing multiple simulations and observing the 
output of each simulation, component variation limits can be estab¬ 
lished to guarantee a functional circuit. Figures 1.2a and b illustrate 
bow multiple simulations may be used to characterize the performance 

°f a circuit. 

Many of the vendor-offered SPICE packages now offer Monte Carlo 
®ud worst-case analyses. Monte Carlo and worst-case analyses auto- 
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While SPICE will not predict layout-dependent parasitic effects 
before they are known, once these effects are known and the parasitic 
elements entered in the netlist, SPICE can be used to help determine 
a solution to the problem. In the example of the CMOS latch-up, once 
the parasitic bipolar transistor has been identified as the cause of 
the failure, a designer may discover that adjusting the beta (Hfe) of 
the transistors may reduce or eliminate the latch-up effect. Once a 
parasitic element has been identified as the cause of a problem, sim¬ 
ulation can be a fast, efficient means to determine how to correct the 
problem. 

These are just some of the reasons that so many designers are turn- 
i„g to simulation. Simulation doe, not replace the breadboard. Often, 
the type of information you obtain from the simulation is not the same 
as the information gathered from a real prototype, and the converse is 
also true. Used together, the breadboard and the circuit simulator form 
an indispensable design aid. 

Which Simulator Is Best? 

This is a question without an answer. Asking which simulator is 
best is like asking which op-amp is best, which bipolar transistor is 
best, or which solder is best. You can order a copy of IntuSoft’s 
IS_SPICE for $99, or you can order a Cray version of Meta-Software’s 
HSPICE for $120,000. Obviously there are differences between the 
two programs but there are also several similarities. Figure 1.1 illus¬ 
trates the functional relationship between SPICE2G.6 and several 
of the vendor-offered simulators. Functionally, most of the SPICE- 
compatible simulators are supersets of SPICE, which means that a 
circuit netlist which runs in SPICE2G.6 will also run in HSPICE, 
IS_SPICE, MICRO-CAP IV, or PSPICE, but the converse is not true. 
Circuits which use HSPICE-specific functions and features will not 
run in SPICE2G.6, nor will they run in PSPICE, MICRO-CAP, or 
IS_SPICE, 

The question you should ask is, “Which simulator is best suited for 
my design needs?” Of course, SPICE2G.6 can still be ordered from 
Berkeley for $150 (source code only). But many of the vendor-offered 
programs offer advantages and capabilities over and above Berkeley 
SPICE. Some of the features added in the vendor-offered tools include 
improved device models, Monte Carlo and worst-case analysis, pole- 
zero analysis, network analysis, scaleable parameters and functions, 
automatic optimization, device model libraries, and hotline support. 
While many of the added features of one vendor-offered tool may over¬ 
lap another, many features are unique to any SPICE-compatible 
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simulator. Look for the features which will best meet your design 
needs. Only you can determine which simulator is best for the way you 
design. 

A word about price 

As with most free-market commodities, the price you pay for a simula¬ 
tor depends on the features and capabilities added to the program, the 
number of copies a vendor expects to sell (personal computer software 
usually sells for less than workstation software because the PC market 
is much larger), and the speed of the machine on which you run the 
software. The last reason may seem a little odd. Many software com¬ 
panies have decided that running a simulator on a faster machine 
means obtaining the results sooner, or being more productive with 
their software. Because of this, many software vendors adjust the price 
of the tools based on the speed of the computer. Shop around. You may 
find you can save a bundle depending on the type of computer you’re 
using. 

What Do You Need to Simulate? 

In the most basic form, simulation requires a text editor to create the 
input files and a simulation program. A minimum machine recommen¬ 
dation is a 386 PC (with the coprocessor). As far as a simulation pro¬ 
gram and editor, the MS-DOS editor (EDIT) is an easy way to create 
your first input files. 

If you are using a PC for simulation, the simulator you start with 
should be a 32-bit version of SPICE. SPICE uses 64-bit math through¬ 
out the program. Many of the early PC versions of SPICE were 16-bit 
programs. To move a single-node voltage or branch current required 
four memory transfers. A 32-bit program moves data almost twice as 
fast as a 16-bit version of SPICE; because of this, a 32-bit version of 
SPICE executes much quicker than a 16-bit version. While the 16-bit 
versions of SPICE will run on any PC (from the pcXT through 486s), 
the 32-bit programs will only run on 386sx, 386dx, 486, and higher 
CPUs. The RSPICE and RGRAPH programs which accompany this 
text are 32-bit programs. Ask your simulation vendor if their program 
is a 32-bit or 16-bit version. 

But the preceding paragraph describes only the basics. If you start 
^th just a simulator and editor, you will quickly realize that accurate 
simulation requires accurate component models and accurate device 
Qtodel parameters. One of the most important features any of the 
vendor-offered SPICE versions can add is a library of common electri¬ 
cal components. Models of bipolar transistors, diodes, op-amps, and 
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thousands of other devices can be stored in a library and called into the 
simulator as needed. Of course, you can develop your own device mod¬ 
els, but developing accurate component models and model parameter 
sets can be a costly, time-consuming process. 

Whether you obtain device models from a library or develop them 
yourself, before you can obtain accurate, reliable, meaningful results, 
you need to know how to use your simulator. Using your simulator 
means learning how to assemble input files in the syntax of SPICE, 
learning how to avoid the problems of nonconvergence, timestep con¬ 
trol, and numeric integration inaccuracies. A common trap new SPICE 
users fall into is assuming the simulation output is always correct. 
Often SPICE will fail to converge on a solution or, worse yet, will con¬ 
verge on the wrong solution! Numeric integration and timestep control 
problems can lead to inaccurate results. Even though the SPICE pro¬ 
gram is over 20 years old, we are still learning which numeric algo¬ 
rithms work well and which ones don’t. 11 

The final ingredient for high-quality simulations is user education. 
Learning how to overcome nonconvergence problems, learning to 
decide whether the simulation output is accurate or erroneous, and 
learning to diagnose simulation ills are the lessons which must be mas¬ 
tered to achieve accurate, efficient simulations. 

These three elements—a good simulator, accurate models, and good 
user skills—are like the legs of a tripod; no one element is more impor¬ 
tant than any other. Without a good simulator, accurate device models, 
and good user skills, simulation results are simply imaginary numbers. 

The obstacles of simulation 

While the preceding paragraphs proclaim how wonderful SPICE is, in 
truth, the average user will have many problems when using the pro¬ 
gram. SPICE is not a foolproof design aid. For many years, the author 
had the job of helping the engineers of a large design community learn 
to use SPICE. Although several SPICE-based packages were in use, 
including SPICE2G.6, many simulations still suffered from nonconver¬ 
gence problems, timestep control anomalies, and numeric integration 
errors. For many designers , the obstacles of simulation are the inability 
to both determine the cause of and correct an unexpected simulation 
result. 

This book is appropriately titled. Inside SPICE: Overcoming the 
Obstacles of Circuit Simulation. The obstacles of simulation are things 
which make the simulator fail, like nonconvergence; things which 
make the simulator produce the wrong answer, like timestep control 
problems; and things which make the simulator produce unexpected 
results, like numeric integration instabilities. Many designers don’t 
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realize SPICE can make mistakes, much less know how to correct a 
mistake if it occurs. These are the real obstacles to circuit simulation, 
imd these are the problems this text attempts to address, identify, and 
correct. 


Summary 

SPICE can be an indispensable design/analysis tool, but, like any tool, 
SPICE has limitations. Even the finest tools are ineffective in the 
hands of an inexperienced user. SPICE is a complicated design tool. 
You cannot simply toss a netlist and device models into a simulator and 
expect correct results. Like a craftsperson skillfully making a fine 
piece of furniture, a simulation user will carefully select the simulator 
options which are appropriate for a given circuit and analysis type. 
Accurate simulation results don’t happen by accident. Both the device 
models and the simulator options you select will determine the speed, 
the accuracy, and the efficiency of your result. 

This text focuses specifically on the algorithms of SPICE2G.6. Since 
the majority of SPICE-like simulators are based on the SPICE2G.6 
algorithms, almost everything presented in this text will apply 
directly to the SPICE-like simulators. While many of the vendors have 
added modifications to the original 2G.6 algorithms, none of the 
SPICE-like simulators has completely eliminated nonconvergence 
problems, timestep control problems, or numeric integration prob¬ 
lems. For these reasons, it is important for all SPICE users to under¬ 
stand why these problems arise and how they can minimized or 
eliminated. All of the solutions presented in this text are in the form 
of standard user-selectable SPICE controls, and most of the SPICE- 
like simulators have these same controls. Teaching SPICE users 
how to tailor the simulator’s control parameters to enhance the speed 
and accuracy of the simulation result is the goal of this text; the 
lessons learned here can be applied to every SPICE-like simulator 
available today. 
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Understanding Circuit Simulation 


The SPICE Engine 

SPICE has been available for over 20 years, yet many simulation users 
understand very little about how SPICE works. Many engineers know 
much more about the parameters of their automobile engine than they 
do about the parameters of the SPICE engine. Test yourself. How many 
cylinders does your engine have? How much horsepower does the 
engine produce? How many valves per cylinder do the heads have? How 
many forward gears are in the transmissions? Many of you will answer 
these without a moment’s hesitation, and very few of us would buy a car 
if we did not know these basic characteristics of the automobile. 

But how would you answer the same type of questions regarding 
your circuit simulator? How many solution algorithms does SPICE use 
to solve the circuit equations? How many timestep control algorithms 
does SPICE offer for transient analysis? Which numeric integration 
routine does SPICE use? What is the smallest voltage or current 
SPICE can resolve? What is the upper limit on resistance in SPICE? 
How would you do on these questions? 

In this chapter, these and other questions about the simulator will be 
answered. Here, the internal workings of SPICE will be exposed. You 
^iU learn how the simulator works, how it represents circuit elements 
ln the system matrix, and how the matrix is filled and solved. This 
chapter is going to look under the hood of a simulation vehicle, at the 

engine of SPICE. 

tinder the Hood 

Many people learn about automobiles by simply raising the hood, look- 
pointing, and asking questions. Unfortunately, software doesn’t 
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lend itself to this same type of troubleshooting. So not surprisingly, 
how circuit simulation works is not widely understood. 

Although it may seem so at times, there is nothing magical about 
SPICE. SPICE will not generate a single solution which could not be 
generated by paper and pencil calculations. The SPICE program sim¬ 
ply automates the calculations. 


The System Equations 

SPICE starts an analysis by writing a set of nodal equations which 
describe the elements in the circuit. Figure 2.1 illustrates a simple 
resistive circuit. The nodal equations for this circuit may be con¬ 
structed by summing the currents leaving each of the three circuit 
nodes, Vi, V 2 , and V 3 . Figure 2.2a and Eq. 2.1 describe the current leav¬ 
ing the Vy circuit node. 


-3 amps + ——— = 0 (2.1) 

5 ohms 

Figure 2.26 and Eq. 2.2 describe the current leaving the V 2 circuit 
node. 


Y *~ Y }- + - Yl -+ V 2-Vs _ o (2.2) 

5 ohms 10 ohms 5 ohms 

Figure 2.2c and Eq. 2.3 describe the current leaving the V 2 circuit 
node. 


Yl s YY l +- Yl -=0 (2.3) 

5 ohms 10 ohms 



Figure 2.1 A simple linear resistive circuit. 
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With simple algebraic manipulations, Eq. 2.1—2.3 may be written in 
the form of Eqs. 2.4-2.6, respectively. 

•2*V, - ,2*V 2 = 3 (2.4) 

-2*Vi + .5*V 2 - .2*V 3 = 0 (2.5) 

~.2*V 2 + ,3*V 3 = 0 (2.6) 



(c) 

2® ur# 2-2(a), (b), and (c) Summing the branch currents to 
Qeier ®me the nodal equations. 
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Equations 2.4-2.6 are the nodal system equations which describe the 
circuit of Fig. 2.1. These three equations have three unknown voltages 
V u V 2 , and V 3 . Gaussian elimination is one of the best-known methods 
for solving simultaneous equations and can be used if we rewrite Eqs. 
2.4-2.6 in matrix form. 


' .2 

-.2 

■ 

0 


■ Vi 


■ ■ 

3 

-.2 

.5 

-.2 

* 

V 2 

= 

0 

0 

* 

-.2 

.3 

■ 


.v 3 . 


,0, 


Equation 2.7 is the matrix form of the system equations. Applying 
forward-elimination to Eq. 2.7 yields Eq. 2.8. (Readers are encouraged 
to derive Eq. 2.8 from 2.7 just for the fun of it.) 


‘.2 

-.2 

■ 

0 


' v,‘ 


v 

0 

.3 

-.2 

* 

V 2 

= 

3 

, 0 

0 

.25 
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. v Sl 


3 

* J 


The value of voltages V u V 2 , and V 3 may be determined by applying 
back-substitution to Eq. 2.8 resulting in: 


y 3 = 3/.25 

= 12 volts 

(2.9) 

y _ [.2*(V 3 ) + 3] 

= 18 volts 

(2.10) 

2 .3 


y _ t.2 *(V 2 ) + 3] 

= 33 volts 

(2.11) 

1 .2 



(Author’s note: It took the author approximately 12 minutes to work 
through this example; it took SPICE exactly .09 seconds to calculate 
the same result!) 


Elements in the Matrix 

For many circuits, SPICE uses simple nodal analysis techniques to 
determine the circuit voltages. But the solution algorithm becomes 
much more complicated when nonlinear devices and charge-storage 
elements are included in the circuit. Nonlinear and charge-storage ele¬ 
ments must be reduced to simplified equivalent circuits before being 
entered in the system equations. These simplifications are required 
because the system matrices only accept linear I-V relationships. To 
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understand how SPICE develops the matrix, the linear I-V character¬ 
istic of each circuit element must be found. 

The matrix 

The system equations in SPICE are represented in a set of matrices. 
The equations characterize the linear representation between the volt¬ 
age and current for every element in the circuit. The system matrices 
are shown in Fig. 2.3. 

The voltage array is the solution array and represents the node volt¬ 
ages of the circuit. When performing a simulation, SPICE tries to 
determine the node voltage values which satisfy KirchhofFs voltage 
and current laws for the circuit. The solution voltages found in the sim¬ 
ulation output come directly from the voltage array. 

The current array is one of the known values of the circuit and rep¬ 
resents the independent branch currents generated from current 
sources and active devices. During simulation, SPICE determines the 
branch currents from the current source settings or from the previous 
voltages applied to the active device terminals. 

The conductance array is another known value of the circuit and rep¬ 
resents the linear relationship between voltage and current for every 
element in the circuit. The values in the conductance array are deter¬ 
mined from the elements of the circuit. Nonlinear elements such as 
diodes and transistors, and charge-storage elements such as capacitors 
and inductors, are represented in the conductance array by their linear 
equivalent circuits. The entries in the current and conductance arrays 
are used to determine the next set of solution voltages. 

Resistors 

Circuit resistors only appear in the conductance array. The value in the 
conductance array represents the inverse of the resistance value. Fig- 
tire 2.4 represents the I-V characteristic of a resistor and the resistor’s 
equivalent conductance value. The resistor’s value changes neither 
with applied voltage nor time, so the value of conductance which repre¬ 
sents the resistor remains constant for the duration of the simulation. 
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***** 2.3 The system matrices used in SPICE. 
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Nonlinear elements 

Representing diodes, transistors, and other nonlinear devices in the 
system equations is slightly more complicated than simple resistors. 
The reason for the complication is that nonlinear devices are defined 
by current-voltage relationships which change in response to the DC 
bias of the device. This is a problem for the system equations because 
the matrices represent linear relationships (G * V = /). The conduc¬ 
tance array represents the linear I-V relationship for every element in 
the circuit. 

The way around this problem is to break the nonlinear I-V relation¬ 
ships into many smaller linear approximations. Figure 2.5a illustrates 
a simple nonlinear function. Figure 2.55 illustrates the same function 
expressed as several linear pieces. Any nonlinear function may be 
expressed as a series of linear (piecewise linear) approximations. As 
the number of linear segments increases, the accuracy of the approxi¬ 
mation increases. SPICE uses this same technique to simulate nonlin¬ 
ear circuit elements. 

During the simulation, SPICE transforms the nonlinear elements 
into simple linear equivalent circuits. The linear equivalent circuit is 
determined from the voltage bias of the device. At each step in an anal¬ 
ysis, SPICE uses the linear equivalent circuit to represent the device 
as it is operating at that voltage bias . 

Figure 2.6a illustrates the I-V characteristics of a diode. If the solu¬ 
tion (voltage) is close to the voltage V b , the only part of the diode I-V 
characteristic of interest to the simulator is the I-V characteristics 
close to V b . In this case, the nonlinear characteristics of the diode may 
be represented by a straight-line approximation equal to the tangent of 
the diode characteristic. 

Figure 2.66 illustrates a straight-line model of the diode. The straight- 
line model is determined from the tangent of the diode I-V characteris- 
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Figure 2.5(a) and ( b) A nonlinear function and a piecewise 
linear approximation of the function. 


tics at the voltage bias V b . The straight-line model intersects the I d axis 
at Ieq gold has a slope equal to G d , the dynamic diode conductance. This 
straight-line model is also known as the linearized diode model. The 
technique of building the linear model is known as linearization. While 
the original nonlinear diode characteristics are expressed by Eq. 2.12, 
Eq.2 .13 represents the linearized diode equation. 


I d =Is 



JZ%_ 

NKT 



( 2 . 12 ) 


I d = G d *V d + I eq (2.13) 

Figure 2.7 illustrates an equivalent circuit representation of the lin¬ 
earized diode equation (2.13). The circuit in Fig. 2.7 is the equivalent 
circuit SPICE uses to represent diodes in the system equations. This 
equivalent circuit is known as the linear diode model. In SPICE, diodes 
QTe modeled by a conductance (resistance) in parallel with a current 
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source. During a simulation, the conductance value G d (1 !G d is also 
known as the small-signal resistance or the dynamic resistance) will be 
stored in the conductance array, and the equivalent current I eq will be 
stored in the current array. In turn, these values will be used to com¬ 
pute a new set of solution voltages. At each solution point, the diode 
current and conductance values are computed and stored within the 
system equations. 

When discussing the linear diode model, the distinction between the 
linear model and the small-signal diode model should be made. 
Although the linear diode model is similar to the small-signal diode 
model, the two are not the same. The linear model of the diode 
contains two elements, the dynamic conductance (G d ) of the diode and 
the equivalent current (J eq ). The small-signal diode model also con¬ 
tains two elements, the dynamic resistance (1 lg d ) and the small-signal 
capacitance ( c d ). The linear diode model may or may not contain the 
small-signal capacitance ( c d ) (depending on whether a transient anal¬ 
ysis is being performed), and the small-signal diode model does not 
contain the equivalent current ( I eq ). Figure 2.8a illustrates the linear 
diode model, whereas Fig. 2.86 illustrates the small-signal diode 
model. 

Transistors are represented by extending the linear diode model. For 
bipolar, JFET, and MOSFET transistors, the linear model is estab¬ 
lished by using the same linearization technique previously described 
for the diode. But in the case of transistors, the linearization process 
must account for every pair of device terminals. For the bipolar, the 
base-emitter, collector-emitter, and collector-base terminals are lin¬ 
earized. For the MOSFET, the gate-source, drain-source, source-drain, 
and source-bulk terminals are linearized. Once the linearized models 
are found, the conductance and equivalent current values are stored in 
the system equations. 



(a) (b) 


^Hiure 2.8(a) and (b) The linear equivalent circuit of a diode and the 
small-signal model of a diode. 
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Charge-storage elements 

While understanding the nonlinear devices is made much easier by 
drawing the I-V characteristics of the device, understanding the charge- 
storage elements is more difficult because the I-V characteristics can¬ 
not be drawn as a simple two-dimensional graph. For capacitors and 
inductors, the I-V characteristics change with the applied voltage, cur¬ 
rent, and time. But, like the linear and nonlinear devices, the capacitor 
and inductor values must be stored in the system equations. 

In a small-signal frequency sweep, inductors and capacitors repre¬ 
sent linear elements, but during a transient time sweep, capacitors 
and inductors possess a definite nonlinear I-V characteristic which is a 
function of time. To visualize the I-V characteristics of a capacitor 
would require a three-dimensional graph, with one axis representing 
voltage, one representing current, and one representing time. 

Figure 2.9 illustrates a simple capacitor circuit. At time T= 0, the 
switch closes and current flows onto the capacitor plates. Over a period 
of time, the capacitor voltage rises to the battery potential, and the cur¬ 
rent flowing into the capacitor decreases to zero. Figure 2.10a illus¬ 
trates the relationship of the capacitor current and time. Figure 2.106 
illustrates the relationship of the capacitor voltage and time. 

To visualize the I-V characteristics of the capacitor, imagine Figs. 
2.10a and b placed together on a common time axis. Figure 2.11a illus¬ 
trates the superposition of Figs. 2.10a and b on a common axis. Notice 
that the common axis in Fig. 2.11a and b is marked in discrete time 
intervals. The linear I-V relationship of the capacitor can be found by 
drawing a line between the capacitor current on the Y axis and the 
capacitor voltage on the X axis. The resulting line represents the linear 
I-V relationship of the capacitor at that moment in time. If you follow 
this procedure for each timepoint, you will witness the I-V relationship 
change over time. This effect is illustrated in Fig. 2.116. Figure 2.12 
illustrates Fig. 2.116 as viewed along the time axis. 


* 


Figure 2.9 A simple time- 
domain circuit. 
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figure 2.10(a) and ( b ) Capacitor current vs, time and capacitor voltage vs. time for Fig. 2.9 circuit. 

During transient analysis, the solution points are discrete moments 
in time. At each solution timepoint, the numeric integration routines 
determine the linear I-V relationship for capacitors and inductors. The 
linear relationship is expressed graphically in Fig. 2.13 and mathe¬ 
matically in Eq. 2.14. 

I C = G*V C + I eq (2.14) 

In SPICE, charge-storage elements are represented with simplified 
linear equivalent circuits. (For the capacitor, the circuit represents the 
current-voltage relationship expressed in Eq. 2.14.) These equivalent 
circuits are known as the companion models. Figures 2.14a and b illus¬ 
trate the companion model of the capacitor and inductor. The conduc¬ 
tance and current source value (or resistance and voltage source value 
for an inductor) of the companion model are recalculated at each new 
transient timepoint by the numeric integration algorithms. Once 
determined, the elements of the companion models are stored in the 
system matrices. At each new solution timepoint in a transient analy¬ 
sis, SPICE computes new companion models for each capacitor and 
inductor in the circuit. 

SPICE simulates the behavior of electrical components by repre¬ 
senting each element as a simple linear equivalent circuit. This is 
known as linearizing the circuit. As the simulation progresses, SPICE 
recalculates new linearized equivalent circuits for each nonlinear and 
charge-storage element in the circuit. As the circuit voltages and cur¬ 
rents change, the linearized models change to reflect the nonlinear 
behavior of the circuit. 

Matrix construction by inspection 

In the section on the system equations in this chapter, you saw a sim¬ 
ple circuit transformed into a set of nodal equations. The nodal equa¬ 
tions were then transformed into a set of system equations and put in 
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Figure 2.11(a) and ( b ) Capacitor current and voltage vs. time and the I-V 
characteristics of discrete timepoints. 


matrix form. SPICE basically uses this same procedure to build the 
system equations and circuit matrix. To save time, though, SPICE uses 
a shortcut to develop the system matrix. This shortcut allows SPICE to 
build the matrix as each element is read from the input file. This short¬ 
cut is known as matrix construction by inspection . Matrix construction 
by inspection builds the system matrices and identifies the location of 
each element in the matrices as soon as the nodes connected to the ele¬ 
ment are defined. 1 
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Figure 2.12 The I-V characteristics of a capacitor as viewed along the time axis. 


Matrix construction by inspection builds the system matrix with the 
help of predefined element templates. The templates describe the posi¬ 
tion in the matrix for the conductance and current values of a particu¬ 
lar device. The templates for the resistor, the current source, and the 
diode are shown in Fig. 2.15a, 6, and c. Every element which can be 
simulated in SPICE has an element template. Even the four-terminal 
MOSFET transistors have a predefined element template. The tem¬ 
plate for the MOSFET is shown in Fig. 2.16. The system matrix in 
SPICE expands as more elements are read from the input file. When 
SPICE reads the element from the input file, the element template 
describes the positions in the conductance and current arrays for the 



Figure 2.13 The linearized model 
for a capacitor voltage and cur¬ 
rent at time 71 
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(b) 


Figure 2.14(a) and ( b ) The companion models of a capacitor 
and an inductor. 

device. The template also determines whether the arrays need to be 
expanded to add additional nodes for new elements. By the time 
SPICE reads the last element from the circuit file, all of the circuit 
components have been assigned proper locations in the system matrix. 

Once the templates for an element are known, matrix construction 
by inspection may be applied to almost any circuit. As a simple exam¬ 
ple, examine Fig. 2.1. A SPICE netlist for the circuit shown in Fig. 2.1 
is shown in Fig. 2.17. 
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for Node A for Node B 



Node B Conductance Voltage Current 

array array array 

Figure 2.15(a) The template of a resistor. 
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Figure 2.15(6) The template of a constant current source. 
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Figure 2.15(c) The template of the diode equivalent circuit. 


As SPICE reads the current source 1 1 , the current source template 
assigns the node connections to the rows and columns in the matrices 
which correspond to the nodes of the source (node 1 and node 0). As 
prescribed by the template, SPICE creates an initial 2x2 system 
matrix as shown in Fig. 2.18a. SPICE continues by reading the resis¬ 
tor R lm Using the resistor template, the matrix is increased by one 
additional node as shown in Fig. 2.186, and the conductance values 
corresponding to the resistor are assigned to the conductance array. 
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figure 2.16 The template of a four-terminal MOSFET. 
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Figure 2.17 SPICE netlist for Fig. 2.1. 


The conductance of a 5-ohm resistor is .2 mhos, and the value of the 
resistor is placed in the conductance array as prescribed by the tem¬ 
plate. The R 2 resistor is added to Fig. 2.18c. Notice how overlapping 
template entries are added to existing entries. The R s resistor is added 
to Fig. 2.18c?, and R 4 is added to Fig. 2.18c. Figure 2.18c is a 4x4 matrix 
which corresponds to the four nodes in the circuit (nodes 1, 2, 3, and 
ground). But since ground is defined as the reference voltage and is 
equal to zero, the matrix row and column which correspond to ground 
may be eliminated. Figure 2.18/* illustrates the final matrix after the 
elimination of the ground row and column. Notice that Fig. 2.18/* is an 
exact match to Eqs. 2.4-2.6, and all this was achieved without writing 
a single nodal equation! 

SPICE uses matrix construction by inspection to initialize the sys¬ 
tem equation matrix. Matrix construction by inspection is a fast and 
efficient manner to construct the system equations. SPICE uses matrix 
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Figure 2.18(a) and ( b ) An example of matrix construc¬ 
tion by inspection. 
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Figure 2.1 B(c) and (cf) An example of matrix construction by 
inspection. 


construction by inspection to develop the system equations as the cir¬ 
cuit components are read from the input file. By the time the last ele¬ 
ment has been read from the file, the system matrices have been 
defined completely. 


The Matrix Solution 

Once the conductance and current arrays have been filled, SPICE 
must solve the node voltage values in the voltage array. SPICE has two 
solution algorithms, one for linear circuits and one for nonlinear cir¬ 
cuits. In this text a linear circuit is defined as one which contains only 
linear elements, voltage sources, and current sources. A linear element 
is defined as a passive element that has a linear voltage-current rela¬ 
tionship. 2 For a DC and transient analysis, linear elements include 
resistors and linear-dependent voltage and current sources. For an AC 
small-signal analysis, linear elements include resistors, capacitors, 
inductors, and linear-dependent voltage and current sources. 

Linear analyses 

If a circuit contains only linear elements, SPICE uses the computer 
equivalent to Gaussian elimination to solve the matrix. The algorithm 
SPICE uses is known as LU decomposition. For the purpose of under¬ 
standing, LU decomposition is an efficient computer method of per- 
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Figure 2AB(e) and (f) An example of matrix construe* 
tion by inspection. 


forming Gaussian elimination to solve a linear set of equations. 1 An 
example of Gaussian elimination was used earlier in this chapter in 
Eqs. 2.7-2.11. 


Nonlinear analyses 

SPICE uses Gaussian elimination to solve the circuit equations when 
only linear elements are present. But, with the introduction of one or 
more nonlinear elements, SPICE must use a different solver algo¬ 
rithm. The second solver is a nonlinear solution technique known as 
the Newton-Raphson algorithm. 

Iterating to a solution. If a circuit contains only linear elements, the cir¬ 
cuit equations are algebraic in nature. Equation 2.15 is an example of 
a simple algebraic equation. 


7x + x = 32 (2.15) 

Equation 2.15 may be solved with simple algebraic manipulations as 
shown in Eqs. 2.16 and 2.17. 

8x = 32 (2.16) 

x = 4 (2.17) 
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But with the introduction of one or more nonlinear elements in the 
circuit, the circuit equation becomes transcendental in nature. Equa¬ 
tion 2.18 is an example of a transcendental equation. 

ZrcOe) + * = 32 (2.18) 

Equation 2.18 cannot be solved with simple algebraic manipulations. 
Instead, transcendental equations are often solved with an iterative 
guessing technique. Begin by rearranging Eq. 2.18 into Eq. 2.19. 

fix) = Inix) + x - 32 (2.19) 

Construct a table as shown in Fig. 2.19a. Equation 2.19 may be 
solved by making a series of iterative guesses until, at the proper value 
of X, the f(x) result is equal to or very close to zero. The author’s suc¬ 
cession of guesses is shown in Fig. 2.196. 

When one or more nonlinear devices is introduced into the circuit, 
the system equations become transcendental equations and cannot 
be solved with simple Gaussian elimination. In this case, a nonlinear 
solution technique known as the Newton-Raphson algorithm is applied 
to the system matrix. The Newton-Raphson algorithm, a method of 
successive approximations, is an iterative approach to solving a set of 
nonlinear equations. This algorithm starts the iterative process with 
an initial guess and finds the solution through a series of successive 
guesses. 

Equation 2.20 is the Newton-Raphson formula. 

Xn+1= X"-f^) (2 - 20) 

A simple example of the Newton-Raphson formula might be helpful 
here. To apply the Newton-Raphson algorithm solution algorithm to 
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Figure 2.19(a) and ( b ) Iterating 
to the solution of a transcenden¬ 
tal equation. 
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Eq. 2.18, rewrite Eq. 2.18 in the form illustrated in Eq. 2.19. This 
will be the function F(x). The solution to this equation is the value of 
X which makes the function F(x) equal zero. Ib apply the Newton- 
Raphson algorithm, the derivative of Eq. 2.19 must be found. Equation 
2.21 represents the derivative F'(x). 

f'(x) = - + 1 (2.21) 

X 


The resulting Newton-Raphson formula is shown in Eq. 2.22. 


X n + 1 =X n - 


ln(Xn) +Xn - 32 



( 2 . 22 ) 


Starting the iterative process with an initial guess of X n = 1 results 
in the table of values shown in Fig. 2.20. Notice that after the fourth 
iterative value, the values of X n and X n + 1 remain constant. The 
Newton-Raphson iterative process begins with an initial guess and ter¬ 
minates when the difference between successive guesses falls to zero. 
When the Newton-Raphson algorithm has found the exact solution, the 
value predicted by any additional iterations remains unchanged. In 
Fig. 2.20, the fourth and fifth iterations are identical and match the 
solution calculated earlier in Fig. 2.196. 


Terminating the iterations. SPICE uses the Newton-Raphson algorithm 
to solve the circuit equations when one or more nonlinear devices are 
entered in the circuit. SPICE starts with an initial guess for every node 
voltage in the circuit and begins iterating. With each successive itera¬ 
tion, a new set of node voltages is predicted. The solution routine mon¬ 
itors the node voltage of the present iteration and the previous 
iteration value. Ideally, at the exact solution, the node voltage between 
successive iterations should be identical, or the difference between 
iterative voltage values should be zero. But, because of the way digital 
computers represent numbers, saying when two numbers are exactly 
equal can be difficult (because of round-off errors). Because of this dif- 
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Figure 2.20 The Newton-Raphson iterations on a 
transcendental equation. 
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ficulty, SPICE monitors the difference between iterative node voltage 
values and compares the difference with a predefined error tolerance. 
When the difference between iterative voltage values is less than the 
error tolerance for every node of the circuit, SPICE terminates the 
iterative process for that solution point. 

Nonconvergence in SPICE. In addition to the error tolerance limits, 
SPICE limits the total number of iterations each analysis type is 
allowed to process. If the iterative node voltages have not satisfied the 
error tolerance requirements before SPICE exceeds the iteration limit, 
SPICE aborts the simulation and proclaims the infamous nonconver¬ 
gence error message. 

Each different analysis type has a different limit on the number of 
iterations allowed. Each analysis iteration limit may be reset by the 
user with the .OPTIONS statement. Chapter 3 will explain how to set 
each of the iteration limits. 

Newton-Raphson example. The Newton-Raphson algorithm is an itera¬ 
tive process which allows SPICE to quickly determine the node voltage 
values of the circuit. To see how SPICE applies the Newton-Raphson 
algorithm to a circuit, look at the circuit shown in Fig. 2.21. 

The diode of the circuit is characterized by the I-V equation shown in 
Eq. 2.23. 


I d = lpA*[exp(40*V d ) -1] (2.23) 

From this, the nodal equation for this one-node circuit may be writ¬ 
ten by summing the currents leaving the V x node. The nodal equation 
is shown in Eq. 2.24. 


F(V d ) = 0 = -5 +-& +lpA*[exp(40*V r <i ) -1] (2.24) 


V", 
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Equation 2.24 forms our circuit function F(V d ). To apply the Newton- 
Raphson formula to this circuit, the derivative of Eq. 2.24 must be 
found. The derivative of Eq. 2.24 is shown in Eq. 2.25. 


-f ff - = ° + + 40pA*exp(40*V r rf ) (2.25) 


The iterative equation for this circuit is shown in Eqs. 2.26 and 2.27. 


X 


n + 


V _ WXn) 

J*jLrt _ 

F' (X n ) 


(2.26) 


y _ y _ -5 + . 5*V d + lpA*[exp(40*y d ) - 1] 
d + 1 d .5 + 40pA*exp(40*V d ) 


(2.27) 


Starting with an initial guess of the voltage at 1 volt, 14 iterations are 
required to reach .7291 volts. The iteration voltage values are shown in 
Table 2.1. Notice the 14th iteration in Table 2.1; both V d and V d + i are 
.7291. The iterative procedure was set to terminate when the voltage- 
to-voltage iterations were equal to within four significant digits. 

In this example, the initial voltage was chosen to be 1 volt. If a dif¬ 
ferent starting voltage had been used, the iterative procedure will still 
iterate to .7291 volts, but depending on the initial voltage, the iterative 
process may require more or less iterations. 

Table 2.2 illustrates the number of iterations required to iterate to 
.7291 volts for a number of different starting voltages. Notice that the 
number of iterations changes substantially with the initial guess. The 
Newton-Raphson algorithm converges quickly to a solution if the initial 
guess is close to the exact solution. The previous statement is important 


TABLE 2.1 Iterative Voltage Values for Diode Circuit 
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TABLE 2.2 Iterations Required vs. Starting Voltage 


Starting voltage # Iterations 


.6 V 

119 

.65 V 

22 

.7 V 

5 

IV 

14 

2 V 

54 

3 V 

94 


because if SPICE does not achieve convergence within the allowed 
number of iterations, the simulation will terminate with a nonconver¬ 
gence failure. 

In the disk included with this text, the ch_2 subdirectory contains a 
file named DIODE.EXE. This simple program contains the diode cir¬ 
cuit of Fig. 2.19 and the Newton-Raphson algorithm. The program 
prompts the user for a starting voltage, then prints the resulting iter¬ 
ative voltage values. This program was used to generate the entries in 
Table 2.2. Experiment with the DIODE.EXE program and different 
starting voltage values. Confirm the entries in Table 2.1. 

Example summary. SPICE uses the Newton-Raphson algorithm as the 
nonlinear solution engine. The iterative process begins with an initial 
guess for every node voltage in the circuit. The iterative process con¬ 
tinues until either the solution values between iterations are identical 
or nearly identical, or the number of iterations exceeds the allowed 
limit. If the latter occurs, SPICE prints a nonconvergence warning 
message and terminates the simulation. 

Because there is a limit on the number of iterations allowed to find 
the solution, selecting an initial voltage close to the exact solution may 
be very important in overcoming nonconvergence problems. Setting an 
initial node voltage will be discussed in detail in Chap. 3. 

SPICE Operation Overview 

SPICE begins an analysis by reading elements from the input file. 
Using matrix construction by inspection and a set of predefined element 
templates, the system equations are described in a set of linear matri¬ 
ces. Once the system matrix has been defined, SPICE performs the 
analyses described in the input file commands. 

DC operating point analysis 

Every analysis begins with a DC operating point calculation. The DC 
operating point calculation establishes the DC bias of the circuit. In 
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addition, the operating point voltages are used as the initial condition 
for all other sweep analyses. 

During the DC operating point analysis, only the DC currents and 
voltages are calculated. The circuit capacitors are modeled as ideal 
open circuits, and the circuit inductors are modeled as ideal short 
circuits. 

To calculate the DC operating point, the matrix entries representing 
voltage and current sources are set to their proper source values. With 
this complete, SPICE makes an initial guess of the node voltage for 
every other node in the circuit and stores the guess in the voltage array. 

When the voltage array has been filled with initial values, SPICE 
calls a routine known as LOAD. The LOAD routine uses the initial 
node voltage guess to calculate the equivalent current and conductance 
(linearized circuit models) for each nonlinear element in the circuit. 
(The LOAD routine also calculates the companion model for the 
charge-storage elements during the transient analysis.) The LOAD 
routine stores the equivalent currents in the current array and the 
equivalent conductances in the conductance array as prescribed by the 
element templates. 

When all the devices have been stored in the current and conduc¬ 
tance arrays, SPICE passes the arrays to the Newton-Raphson solver. 
The Newton-Raphson solver uses the matrix form of Eq. 2.18 to calcu¬ 
late a new set of node voltages. This new set of node voltages is the first 
set of iterative voltage values, and these values replace the previous 
voltages in the voltage array. Once the new node voltages are stored, 
SPICE again calls the LOAD routine to calculate new linearized circuit 
models based on the new node voltage values. The resulting conduc¬ 
tance and current values are stored in the current and conductance 
arrays, and the iterative solution process continues. 

This cycle of performing a Newton-Raphson iteration, calling the 
LOAD routine to calculate new current conductance values based on 
the previous iterative voltage values, and then performing another 
Newton-Raphson iteration, occurs over and over again. The iterative 
cycle continues until either all the node voltages and branch currents 
match the previous iterative values, or until the number of Newton- 
Raphson iterations exceeds the allowed number of iterations. 

To calculate the DC bias point, SPICE uses an initial guess for the 
node voltages in the circuit. The initial guess is used as a starting point 
in the Newton-Raphson iterations. Normally, the DC bias point solu- 
tion requires somewhere between 10 and 500 iterations. If the solution 
is found, SPICE prints the bias voltages in the output file. If the solu¬ 
tion is not found before the allowed number of iterations, SPICE prints 
a nonconvergence failure message in the output file and aborts the 
simulation. 
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DC sweep analysis 

The DC sweep analysis is simply a series of DC operating point calcu¬ 
lations. In the DC sweep analysis, a source voltage or current is 
stepped over a range of values. At each step in the progression, SPICE 
performs the same DC bias point calculation described previously. But 
in the sweep analysis, each step may require from 2 to 50 iterations to 
obtain the solution just for that point! 

Because the sweep analysis may require hundreds or thousands of 
solution points, the simulation usually takes much longer to complete 
than the DC operating point calculation. But the authors of SPICE rec¬ 
ognized something which would help speed the sweep analysis. 
Because the sweep consists of many evenly spaced steps, the node volt¬ 
age changes between any two steps is usually small. Knowing this, the 
authors of SPICE decided to use the previous solution point node volt¬ 
ages as the initial guess to the Newton-Raphson iterations for the next 
solution point. This technique works well and significantly reduces the 
number of iterations required to solve each step in the DC sweep. 

As the sweep progresses, at each new solution point SPICE saves the 
voltages and currents named on the .PRINT statement in memory. 
When the analysis is complete, the output routines in SPICE print the 
analysis results in either tabular or graphical form. 

AC frequency sweep analysis 

Like the DC sweep analysis, the AC frequency sweep starts by calcu¬ 
lating the DC bias of the circuit. Once the bias has been established, 
the nonlinear large-signal transistor and diode models are replaced by 
their linear small-signal models. The small-signal models are deter¬ 
mined from the DC bias point of the circuit. 

The small-signal models replace the large-signal models because in 
SPICE the AC frequency sweep is defined as a small-signal linear 
analysis. This means that during the AC frequency sweep, distortion, 
clipping, saturation, and other nonlinear effects are ignored after the 
bias point has been established. To SPICE users who are unaware of 
the use of small-signal models during the AC frequency sweep, the out¬ 
put results are sometimes surprising. (For example, if you perform an 
AC frequency sweep on an op-amp with an open-loop gain of le6, a sup¬ 
ply voltage of+15VDC to -15VDC, and an input voltage magnitude of 
1VAC, SPICE will predict an output voltage magnitude of 1,000,000 
volts. To see the op-amp clip at the power supply voltages, a DC sweep 
analysis or a transient times sweep analysis would be required.) 

Small-signal analysis also implies the use of complex quantities like 
voltage phase and magnitude. For SPICE to simulate complex quanti¬ 
ties, the voltage, current, and conductance arrays must have both a 
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real and an imaginary component. During AC analysis, the system 
matrices become complex quantities. 

In addition to storing the DC models, the circuit capacitor and indue- 
tor impedances are added to the complex conductance arrays. 

Once the complex currents, voltages, and impedances of the circuit 
elements are stored in the solution arrays, SPICE solves the system 
equations at each frequency point in the"analysis. At each frequency 
point, the frequency-dependent impedances are calculated and stored 
in the conductance array. The solution algorithm then determines the 
complex voltages and currents which satisfy the circuit equations. 

The invocation of the linear models simplifies the task of finding the 
solution to the circuit equations. (Remember, SPICE only uses the non¬ 
linear solver when nonlinear elements are present.) In the AC fre¬ 
quency sweep, after the bias point has been found, SPICE uses the 
computer equivalent of Gaussian elimination (LU decomposition) to 
solve the linear system equations at each frequency point. Because the 
arrays contain only linear elements, SPICE solves the small-signal 
voltages in one Gaussian decomposition (SPICE does not iterate on the 
AC small-signal solution), and because of this, the AC analysis tends to 
be much faster than any other sweep analysis. 

Transient time sweep analysis 

Of all the analysis types, the transient analysis is the most compli¬ 
cated. In many aspects, the transient analysis is similar to the DC 
sweep analysis but adds several complicating factors to the analysis. 

The transient analysis starts as with a DC bias point calculation. 
During this calculation, the circuit capacitors are modeled as open cir¬ 
cuits and the circuit inductors are modeled as short circuits. The bias 
point of the circuit describes the state of the circuit at time T =0 during 
the transient analysis. 

But the instant after the bias point has been found, the analysis 
changes. For every analysis point after the DC bias point, the time- 
dependent capacitor and inductor impedances must be added to the 
system equations. To do this, SPICE uses a numeric integration rou- 
tine to transform the circuit inductors and capacitors into simplified 
equivalent circuits. The equivalent circuit of a capacitor or an inductor 
represents the instantaneous I-V relationship of the device, and the 
equivalent circuit model changes at every timepoint in the analysis. At 
each new timepoint, the values from the capacitor and inductor equiv¬ 
alent circuits are stored in the current and conductance arrays. 

After the DC bias has been found, SPICE switches into the transient 
solution mode. The transient solver is similar to the DC sweep solver 
except with the addition of the capacitor and inductor equivalent cir¬ 
cuits. At each timepoint in the analysis, SPICE goes through a series of 
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Newton-Raphson iterations. As in the DC sweep analysis, the starting 
point for the transient iterative search is the previous set of node volt¬ 
ages from the last solution timepoint. As the transient analysis pro¬ 
ceeds, SPICE uses the previous solution voltages as an initial guess for 
the next series of iterations. The iterative process continues until the 
solution voltages for that timepoint are found, or the number of itera¬ 
tions exceeds the allowed limit. As the solution voltages for each time- 
point are found, the voltage and current values named on the .PRINT 
statement are saved in memory. After the transient analysis is com¬ 
plete, SPICE prints the stored solution points in the output in either 
tabular or graphical form. 

Summary 

On the surface, circuit simulation is a complex process of numerical 
analysis. But if you break each of the procedures into simple blocks, 
the process becomes much more understandable. Any of the analysis 
results which SPICE computes could also be predicted with simple 
hand calculations (although it might take years to do the same number 
of calculations SPICE can do in seconds). There is nothing magical 
going on in SPICE. 

The focus of this text is to show SPICE users how each of the algo¬ 
rithms works, why sometimes they don’t work, which type of circuits 
lead the simulator to highly accurate results, and which circuits lead 
the simulator to erroneous results. Maybe most important of all, in this 
text SPICE is seen as a tool, a tool with a unique set of capabilities and 
a unique set of limitations. 

Like every tool, SPICE has limitations and does not always produce 
the desired result. (Using a hammer as a screwdriver rarely leads to 
the desired result either.) Whether the tool is an automobile engine or 
a hammer drill or a circuit simulator, to use a tool properly the abilities 
and limitations of the tool must be understood. From the user’s point of 
view, the limits of software tools are often difficult to define, because 
looking inside software is nearly impossible, and in many cases illegal. 

This text examines SPICE and the internal workings of the simula¬ 
tor. Understanding how SPICE works, understanding the abilities and 
limitations of the simulator, and understanding why and where the 
simulator introduces inaccuracies are the keys to using SPICE effec¬ 
tively and productively. 
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Nonconvergence 


Understanding Nonconvergence 

Many SPICE users have, at one time or another, observed a simulation 
which fails to converge. Nonconvergence is the failure of the nonlinear 
solution algorithm; it means the simulator failed to find a set of node 
voltages and branch currents which conform to Kirchhoff’s voltage and 
current laws. Nonconvergence is probably the most persistent and 
frustrating problem facing simulation users. 

But what causes nonconvergence? Is it caused by the solution algo¬ 
rithm or by the type of analysis being performed? Do the models play a 
role in nonconvergence? And, the most important question of all, what 
can be done to minimize or eliminate nonconvergence? 

Although computer simulation has been with us for over 20 years, 
and nonconvergence has been with us since the beginning of computer 
simulation, very little has been written about the causes and cures of 
nonconvergence. Many experienced simulation users learned to over¬ 
come nonconvergence problems with years of experience cultivated 
from trial and error experimentation. But learning through trial and 
error is inefficient, ineffective, and very frustrating. 

In this chapter, the mechanisms which cause nonconvergence will be 
explored and explained. Each section of this chapter will examine a dif¬ 
ferent nonconvergence mechanism. Once the nonconvergence mecha¬ 
nisms are known, a simple, methodical, systematic technique will be 
illustrated to eliminate the nonconvergence catalists. 

Convergence and 

the Newton-Raphson Algorithm 
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SPICE converges on a set of node voltages will lead to a better under¬ 
standing of why SPICE sometimes fails to converge. Chapter 2 demon¬ 
strates how the Newton-Raphson algorithm is applied to a set of 
nonlinear circuit equations, but there is a lot to be learned by looking 
at the Newton-Raphson process graphically. 

Figure 3.1 shows the simple diode circuit used in Chap. 2 to illus¬ 
trate the Newton-Raphson iterative solution technique. This same cir¬ 
cuit will be used to illustrate the Newton-Raphson process graphically. 

Figure 3.2a illustrates the I-V characteristics of the circuit diode, 
and 3.2 b illustrates the load line imposed by the resistor and current 
source. If Figs. 3.2a and b are superimposed, the voltage at node VD 
may be determined from the intersection of the diode I-V characteris¬ 
tics and the load line as shown in Fig. 3.3. This is a graphical technique 
for finding the voltage at node VD. 

In SPICE, the voltage VD is found with the Newton-Raphson algo¬ 
rithm. Figure 3.4a illustrates the superposition of the diode I-V char¬ 
acteristics and the load line of the diode circuit. The Newton-Raphson 


V* 



l d = 1 pA * [exp(40 * V d ) -1 ] 


Figure 3.1 A simple diode circuit. 




Figure 3.2 The I-V characteristic of Fig. 3.1 circuit diode (a); and the load line of Fig. 3.1 
linear elements (6). 
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l d = 4.638A 



Vfj — .7291V 


Figure 3.3 The solution of Fig. 3.1 circuit is found at 
the intersection of the load line and I-V characteristic 
of the diode. 


algorithm begins searching for the intersection of the two curves with 
an initial voltage guess. After the initial guess (Vo) has been chosen, 
SPICE passes the V 0 voltage to the LOAD subroutine. The LOAD 
subroutine determines the linear model of the diode at the voltage 
Figure 3.46 illustrates the linear model of the diode superimposed 
on Fig. 3.4a. Once the linear model has been defined, the Newton- 
Raphson algorithm calculates the second iterative voltage value (Vi) 
by locating the intersection of the linear model and the circuit load line 
as shown in Fig. 3.4c. 

The iterative process continues as the linear model of the diode is 
recomputed for the voltage Vi. The new linear model is again extended 
to the circuit load line to predict the third iterative voltage value ( V 2 ) 
as shown in Fig. 3.4d. 




Sjj* 3-4 (a) The superposition of the I-V characteristic of the diode and the load line and the 
^“toal iterative guess V 0 ; (6) the linear approximation to the diode characteristics at V 0 - 
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Figure 3.4 (c) The linear approximation to the diode characteristic and the load line deter¬ 
mine the iterative voltage Vi; (d) the linear approximation to the diode curve at Vi and the 
load line determine V 2 . 


The Newton-Raphson iterations continue, and eventually the itera¬ 
tive voltage values approach the intersection of the diode current and 
the circuit load line. As the iterative voltage values approach the exact 
intersection of the two curves, the difference between sequential volt¬ 
age iterations decreases, as shown in Fig. 3.5. SPICE converges on the 
solution voltage and stops iterating when the voltage-to-voltage itera¬ 
tions are within a predefined error tolerance. Figure 3.6 illustrates 
the iterative node voltage values and the error between iterations. The 
error tolerance is an important part of the solution algorithm and 
defines how close SPICE converges to the exact solution before termi- 


1.8 

Iterative 1.5 
Node 

Voltage 1 2 

.9 
.6 
.3 
0 

Number of Iterations 

Figure 3.5 Iteration-to-iteration voltage values of the diode circuit 
with an initial voltage estimate of 1.5 volts. C Reprinted from Suc¬ 
cessfully Simulating Circuits with SPICE. Used with permission .) 
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jtating the iterative process. The error tolerances are defined by the 
user on the .OPTIONS statement. 

Causes of Nonconvergence 

There are several reasons the simulator may fail to converge. Some of 
the causes may be traced to the Newton-Raphson algorithm, some to 
the type of analysis being performed, and others to the device models. 
Since there are several nonconvergence mechanisms, learning to iden¬ 
tify where the nonconvergence came from is an important lesson in 
overcoming nonconvergence failures. 

Overcoming Nonconvergence Failures 

Seasoned users will testify to the problems of nonconvergence while 
newer users may only have seen problems occasionally. For those who 
have not seen nonconvergence, simulate the circuit file ch3-l.cir with 
the command: 

SIM CH3-1.CIR 

Then look at the output file. In this simulation, SPICE will not be able 
to determine the proper DC bias of the circuit. 

The ch3-l.cir is a typical nonconvergence failure, and the noncon¬ 
vergence can be corrected with the proper application of .OPTIONS 
statement parameters. In fact, somewhere between 75 to 90 percent of 


Node 

Voltage 



Error 

between 

Iterations 


J«ur.3.6 Iteration-to-iteration voltage values and the error between iterations for the 
<hode circuit. (Reprinted from Successfully Simulating Circuits with SPICE. Used with 
Permission.) 
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all nonconvergence occurrences may be eliminated with proper applica¬ 
tion of the parameters in the .OPTIONS statement or model parameters 
on the .MODEL statement. 

Unfortunately for many users, the .OPTIONS statement parameters 
have never been clearly defined, nor have the true meaning of the 
parameters been explained. Due to the lack of good documentation, 
many users simply rely on trying a variety of option settings until one 
of them will correct the nonconvergence problem. But this technique is 
neither very scientific nor very reliable, and this technique is only 
marginally successful. 

In this chapter, the causes of nonconvergence will be explained and, 
for each of the causes, simple step-by-step procedures will be suggested 
to correct the cause of the problem. This chapter focuses on resolving 
the problems of nonconvergence. 

Problem-solving approach 

In addressing the problem of nonconvergence, several different causes 
must be examined. Some nonconvergence problems are caused by the 
Newton-Raphson algorithm as it is applied in SPICE. These problems 
will be addressed in the section titled “General Newton-Raphson Con¬ 
vergence Aids,” and the techniques presented here will apply to every 
nonlinear analysis in SPICE. The General Newton-Raphson Conver¬ 
gence Aids will reduce nonconvergence problems during the DC operat¬ 
ing point analysis, the DC sweep analysis, and the transient time 
sweep. The General Newton-Raphson Convergence Aids define the 
error tolerance for convergence and set the minimum and maximum 
allowed conductance values of the circuit. Many nonconvergence prob¬ 
lems will be eliminated simply by setting the General Newton-Raphson 
Convergence Aids parameters. Because these convergence aids apply to 
all analysis types, the General Newton-Raphson Convergence Aids 
should be applied to every simulation before the analysis begins. 

But the General Newton-Raphson Convergence Aids will not elimi¬ 
nate all of the nonconvergence problems designers see. Many noncon¬ 
vergence failures can be linked to the type of analysis being performed. 
For this reason, the sections titled “Nonconvergence and the DC Oper¬ 
ating Point Solution,” “Nonconvergence and the DC Sweep Analysis,” 
and “Nonconvergence and the Transient Analysis” will look at the roof 
causes of nonconvergence which are related to specific analysis types 
Again in each of these sections, the cause of nonconvergence will be 
discussed and simple step-by-step procedures will be presented tc 
eliminate the cause of the nonconvergence. 

In essence, this chapter is a cookbook approach to resolving noncon 
vergence problems. First, set the General Newton-Raphson Conver 
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re nee Aids parameters. Then, if your simulation fails because of non- 
jonvergence, determine which analysis type SPICE was performing 
vhen the nonconvergence occurred and follow the procedures outlined 
‘or that analysis type. 

general Newton-Raphson Convergence Aids 

rhe General Newton-Raphson Convergence Aids are a set of .OPTIONS 
jtatement parameters and model parameters which should be set to 
ippropriate levels for your circuit. The General Newton-Raphson 
Convergence Aids define the accuracy SPICE must achieve to converge 
>n a solution and stop the iterative procedure. The General Newton- 
[taphson Convergence Aids also determine the minimum and maxi- 
num allowed circuit resistance values; these values will determine how 
quickly SPICE converges to the proper solution. 

The acronym SPICE stands for Simulation Program with Integrated 
Circuit Emphasis. During the late 1960s and early 1970s, breadboard- 
ng was found to be nearly impossible for many integrated circuits, and 
die need to analyze integrated circuit behavior was crucial to the 
levelopment of a successful design. Much of the funding for the devel- 
ipment of SPICE came from the integrated circuits industry, and 
therefore many of the default values in SPICE were selected for these 
ypes of circuits. 

But today, in addition to integrated circuits, SPICE is used to simu¬ 
late discrete, board-level, and even high-power circuits. Many of the 
program’s default settings are inappropriate for these types of circuits, 
rhe General Newton-Raphson Convergence Aids reset the program’s 
iefault settings to levels appropriate for the circuit being simulated. 

The General Newton-Raphson Convergence Aids consist of learning 
to set the accuracy error tolerances (RELTOL, VNTOL, and ABSTOL), 
learning to set the minimum conductance GMIN, and learning to set 
the series resistance parameters of the semiconductor models. 

Accuracy error tolerances 

h the early versions of CANCER and SPICE1, convergence was 
achieved when the program found a set of node voltages which satis¬ 
fied KirchhofFs voltage law or was sufficiently close to those* voltages. 
For diodes, bipolar transistors, and other devices which possess an 
exponential I-V characteristic, a small change in voltage may produce 
1 large change in device current. This behavior is shown in Fig. 3.7. 
Neither CANCER nor SPICE 1 checked for convergence of the nonlin¬ 
ear device branch currents. Often, simulating circuits which contained 
exponential I-V characteristics led to results which did not satisfy 
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Figure 3.7 For forward-biased 
PN junctions, a small change in 
voltage may lead to a large 
change in current. 


KirchhofFs current law; because of this, later versions of SPICE 
(including all versions of SPICE2) added a mechanism to check for cur¬ 
rent convergence in addition to the mechanism which checks for volt¬ 
age convergence. This means that not only do the voltage-to-voltage 
iterations have to converge on a solution, the nonlinear device branch 
currents must converge on a solution as well. This also means that 
error tolerances for both voltage and current must be defined. 

The error tolerances are important for two reasons. First, the error 
tolerances define the accuracy of the solution. Second, the error toler¬ 
ances define how many iterations are required to find the solution. 
Figure 3.8 illustrates how a different error tolerance will require more 
or fewer iterations to achieve a given level of accuracy. 

In SPICE, both the voltage error tolerance and the current error 
tolerance are composed of a relative limit and an absolute limit. Figure 
3.9 illustrates the procedure SPICE uses to check for convergence dur¬ 
ing the iterative process. 

During the iterative process, as SPICE gets closer to the exact set of 
node voltages which solves the circuit equations, the per-iteration volt¬ 
age change and the per-iteration current change become very small. 
Ideally, SPICE would stop iterating when the per-iteration voltage 
change and the per-iteration current change fell to zero. In this case, 
the voltage and current values between iterations would be exactly 
equal. But because of the round-off errors associated with digital com- 
puters, SPICE defines convergence when the per-iteration node volt- 
age for every node and the per-iteration current for every nonlinear 
branch is less than the error tolerance defined by Fig. 3.9. 

After every Newton-Raphson iteration, SPICE checks the error 
between the previous and the present iterative node voltages, and also 
checks the error between the previous and present iterative branch 
currents. Once SPICE has found a solution which produces nearly 
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Number of Iterations 

Figure 3.8 To reach higher levels of solution accuracy, more iterations are required. 
(Reprinted from Successfully Simulating Circuits with SPICE. Used with permission .) 


identical results between two or more iterations, the iterations are 
stopped, and the solution voltages and currents are saved. If a solution 
is not found within the allowed number of iterations, SPICE termi¬ 
nates the iterative process, prints the nonconvergence warning mes¬ 
sage, and halts the simulation. 


// NC is the number of nonconvergence circuit nodes 
// or nonlinear branch currents 

ITERATION_NUMBER = ITERATI0N_NUMBER+1 
NC = 0 

DO 1=1, # CIRCUIT NODES 

IF ( | V(n) - V(n-l)| > RELT0L*V(n) + VNTOL ) NC=NC+1 
ENDDO 

DO J=l, # NONLINEAR BRANCH CURRENTS 

IF ( 11 (n) - I (n-1) | > RELTOLM(n) + ABSTOL ) NC=NC+1 
ENDDO 

// If all the nodes and branches have converged, we have a 

// valid solution point. If not, perform one more Newton 

// iteration. 

IF (NC = 0)THEN 

GO TO SAVE_S0LUTI0N_P0INT 
ELSE 

IF (ITERATION_NUMBER < ITERATION_LIMIT) THEN 
GO TO PERFORM_ANOTHER_NEWTON_ITERATION 
ELSE 

GO TO S0LUTI0N_FAILED_T0_C0NVERGE 
ENDIF 
ENDIF 


^*9ur® 3.9 Checking for a convergent solution voltage. 
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The relative error tolerance. The error tolerance is defined by the 
.OPTIONS statement parameters RELTOL, VNTOL, and ABSTOL. 
RELTOL defines the relative error tolerance for convergence. The rel¬ 
ative error tolerance requires the iterative process to continue until the 
per-iteration voltage change and the per-iteration current change is 
less than a percentage of the final result. For example, RELTOL is set 
by default to .001 (.1 percent). For a 5-volt circuit node, the node volt¬ 
ages predicted between iterations must be 5mV or less to satisfy the 
relative error tolerance convergence criterion. When the node voltage 
is within 5mV of previous iteration, SPICE will terminate the iterative 
process. A 5mA branch current must converge to 5uA or less to termi¬ 
nate the iterative process. RELTOL may be reset by the user with the 
.OPTIONS statement. 

The absolute error tolerance. In Fig. 3.9, the relative error tolerance 
RELTOL is complemented with the absolute error tolerances VNTOL 
and ABSTOL. The absolute error tolerances are necessary because 
when a node voltage or branch current approaches or crosses zero, the 
relative error tolerance approaches zero, which means the simulation 
result must be infinitely accurate. To avoid problems under these 
conditions, the authors of SPICE added the absolute error tolerances 
VNTOL and ABSTOL to the convergence check. 

As shown in Fig. 3.9, to achieve convergence, the per-iteration volt¬ 
age change and the per-iteration current change must be less than the 
sum of the relative error tolerance and the absolute error tolerance. 
But VNTOL has a default value of only luV and ABSTOL has a default 
value of only IpA. Under normal operating conditions, the relative 
error tolerance is much larger than the absolute error tolerance, and, 
under these conditions, convergence is defined by the relative error 
tolerance. But when a node voltage or branch current falls to or 
approaches zero, the relative error tolerance falls to zero too. Under 
these conditions the absolute error tolerances define when the itera¬ 
tions stop. 

Setting the error tolerances. For many types of circuits, the default 
setting of RELTOL produces an acceptable amount of accuracy and a 
reasonable simulation run time. A few circuits will require a more 
accurate solution (smaller RELTOL), and some may require less 
accuracy (larger RELTOL). The value of RELTOL can be set on the 
.OPTION statement as shown below. 

.OPTIONS RELT0L=.0001 

Changing the value of RELTOL will change the amount of time a 
simulation requires. A smaller RELTOL increases the accuracy of the 
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result but takes longer to simulate because of the extra iterations 
required to achieve the extra accuracy. A larger RELTOL decreases 
the accuracy of the result but simulates faster. As a general guide¬ 
line, decreasing RELTOL by a factor of ten (more accurate result) 
approximately doubles the number of iterations required to solve the 
circuit. Increasing RELTOL by a factor of ten (less accurate result) 
approximately halves the number of iterations required to solve the 
circuit. 

The default values for RELTOL, VNTOL, and ABSTOL were set to 
produce an acceptable result for the type of integrated circuits being 
simulated in the late 1970s. (Remember SPICE is an acronym for Sim¬ 
ulation Program with Integrated Circuit Emphasis.) But today, SPICE 
is used to simulate many different types of circuits including discrete, 
board-level, and high-power circuits. While RELTOL produces the 
same percentage of accuracy regardless of the voltage and current lev¬ 
els in a circuit, circuits which contain voltage and current levels much 
higher or much lower than typical integrated circuits necessitate reset¬ 
ting VNTOL and ABSTOL to appropriate levels. (Do you really need to 
simulate your 20-volt switching power supply to an accuracy of luV 
and IpA?) 

The first of the General Newton-Raphson Convergence Aids is proper 
selection of the simulation error tolerances. To set the error tolerances, 
follow the procedure outlined in Fig. 3.10. Setting the error tolerances 
requires knowledge of the voltage and current levels of the circuit, but 
setting the error tolerances to levels appropriate to your circuit will 
result in faster simulations and fewer nonconvergence failures. 


Once RELTOL has been determined, 

set VNTOL = RELTOL * V small 

Where V Bma ii is the smallest voltage (magnitude) of interest in the circuit. 

(For example, V sma ii might be an input offset voltage for an op-amp circuit or a low 
logic level for a digital circuit.) 

Once VNTOL has been determined, 

set ABSTOL = RELTOL * I SIHll 

Where I 8mall is the smallest current (magnitude) of interest in the circuit. 

(For example, I sma n might be an input offset current for an op-amp circuit or the 
reverse leakage current of a diode circuit.) 

Then, enter the values on the .OPTION statement (.OPTIONS RELTOL=.0001 
VNT0L=.0001V ABSTOL=10NA) 

Figure 3.10 Setting the SPICE error tolerances. 
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Circuit conductance values 

Like the error tolerances, the conductance values of the circuit also 
determine how quickly SPICE converges on a solution. But too few 
SPICE users realize the importance of the conductance terms. When 
SPICE begins the iterative process, the solution voltages and currents 
must be found within an allowed number of iterations, and, although 
the maximum and minimum conductances have a minor impact on the 
accuracy of the circuit, the upper and lower limits on conductance will 
influence how quickly the Newton-Raphson algorithm converges to a 
solution. 


y y _ EiYA . ( 3 . 1 ) 

n + 1 " G(V n ) 

Equation 3.1 is the equation of the Newton-Raphson algorithm when 
applied to a single node circuit. In Eq. 3.1, the denominator of the 
second term is the conductance of the circuit elements. The Newton- 
Raphson algorithm uses conductance values to predict where the next 
voltage iteration will fall. If the conductance becomes very small, the 
second term of Eq. 3.1 ( F(V n )/G(V n )) becomes quite large. Worse yet, if 
the conductance value ever reaches zero, the next Newton iteration 
will cause a floating divide-by-zero error and crash the program. 

Figure 3.11 illustrates the diode characteristics under reverse bias 
conditions. In this region, SPICE models the I-V characteristics as a 
constant current equal to -IS. When the current is constant and no 
longer a function of the diode voltage, the conductance of the diode goes 
to zero. Under these conditions, the Newton-Raphson algorithm should 
fail because the diode conductance is zero. 



Figure 3.11 The diode conductance (g d ) falls to zero in the 
reverse-bias region of operation. 
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Since all of the semiconductor devices contain one or more regions of 
zero conductance (constant current output), the conductance problem 
had to be resolved. The authors of SPICE resolved this problem by 
placing a shunt resistor in parallel with every PN junction of every 
semiconductor model in SPICE. The resistor has a default value of 
1000G ohms, and the value of the resistor is set with the .OPTION 
GMIN=X. GMIN represents the conductance (reciprocal of the resis¬ 
tance) of the shunt resistor. The GMIN resistor is built into the model 
equations and always provides a small voltage-dependent current for 
the devices. Equation 3.2 is the diode equation under reverse bias with 
the GMIN resistor term included, and Eq. 3.3 is the conductance of the 
diode. 


-IS + V/GMIN 

(3.2) 

dId = 0 + GMIN 

(3.3) 

dV d 


The GMIN resistor is found across the diode, the base-emitter and 
base-collector junctions of the bipolar, the gate-drain and gate-source 
junctions of the JFET, and the drain-bulk and source-bulk junctions of 
the MOSFET. Every PN junction of every device in SPICE includes the 
GMIN resistor. Normally the resistance value is so high that the 
resulting current is significantly below the error tolerances of the 
simulation. This means the current through the GMIN resistor does 
not contribute to the value of the simulation result. For example, if 
a circuit diode is held in at -5V reverse bias with luA reverse satura¬ 
tion current, the current flowing through the GMIN resistor is 5pA 
(5V/1000G ohms = 5pA). 5pA is less than the InA relative error 
(RELTOL) tolerance of SPICE (luA * .001 = InA). The default value 
of GMIN was chosen so the simulation accuracy would not be affected 
by its presence. 

But how does GMIN influence the convergence characteristics of 
the circuit? As long as the conductance term in Eq. 3.1 is finite, the 
Newton-Raphson algorithm will continue iterating to a solution, but if 
the conductance is very small, the second term of Eq. 3.1 (F(V n )/(V n )) 
becomes quite large, and the next iterative voltage value will be far 
from the previous voltage value as illustrated in Fig. 3.12. Often, very 
small conductance values force the Newton-Raphson algorithm to sig¬ 
nificantly overshoot the correct solution voltage. When this occurs, the 
Newton-Raphson algorithm will require dozens of iterations to work 
back to the correct solution voltage. 

This effect can be demonstrated with the use of the DIODE.EXE pro¬ 
gram found in the ch2 subdirectory of the text floppy disk. Starting 
with an initial diode voltage of .5 volts, the first iteration is 9.6 volts! 
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Figure 3.12 Small conductance values (large resistances) lead to 
large iteration-to-iteration voltage changes. 


To reach the final solution voltage of .7291 requires an additional 360 
iterations. With an initial diode voltage of .55 volts, the first iteration 
is 7.8 volts and requires 289 additional iterations to reach .7291 volts. 
An initial voltage of .6 volts generates 3.6 volts on the first iteration 
and requires 118 more iterations to reach .7291 volts. 

SPICE must find the proper solution value within a fixed number of 
iterations. Because of this restriction, finding the solution with as few 
iterations as possible is an important aspect of achieving convergence. 
For simulation users, the larger the value of GMIN, the faster the 
Newton-Raphson algorithm will converge on a solution. To help avoid 
nonconvergence problems , GMIN should be set as large as possible 
without affecting the accuracy of the simulation output. Raising GMIN 
decreases the size of the shunt resistor. As long as the shunt resistor’s 
current contribution is lower than the relative error tolerance current 
resolution, there will be no difference in the accuracy of the simulation 
result. 

Figure 3.13 illustrates how to select a value of GMIN that is appro¬ 
priate for your circuit. 

GMIN is a global setting for the entire circuit, so when selecting 
the value of GMIN, consider the most current sensitive portions of 
the circuit. 

Just as the Newton-Raphson algorithm has problems with very 
small values of conductance, very large values of conductance can lead 
to problems also. Figure 3.14 illustrates the I-V characteristics of a 
diode under high forward bias conditions. For very large values of con¬ 
ductance, the second term of Eq. 3.1 becomes very small. Because of 
this, the next voltage iteration (V„ + i) will be only slightly different 
than V n , and if the device is biased far from the proper solution voltage, 
many iterations will be required to work back to the proper solution. 
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To Set GMIN 

Determine the smallest parasitic resistance value ( Rp ) which could be placed across 
any two nodes without influencing the behavior of the circuit. 

SetGMIN = 1/Rp 

Enter the value of GMIN on the .OPTIONS statement (.OPTIONS GMIN-1E-10) 
Figure 3.13 Procedure to set the GMIN value. 


High conductance values are extremely troublesome when diodes or 
PN junctions have a forward bias of more than .8 volts. Beyond this 
point, the conductance values become unrealistically large (resistance 
becomes very small) and will lead to nonconvergence problems. While 
very few circuits require diode biases this high, during the iterative 
process, SPICE may apply a high forward bias to any of the diodes in 
the circuit while iterating to a solution. If this occurs, SPICE may fail 
to converge before reaching the proper circuit solution. 

But SPICE users can guard against the extremely large values of 
conductance seen in diodes and PN junctions by always specifying the 
series resistance model parameter for all circuit diodes and bipolar 
devices. (While it is a good habit to always set the series resistance for 
the MOSFET and JFET devices, these are less prone to forward bias¬ 
ing the internal PN junctions of the device.) The default value of series 
resistance is zero ohms, and if during the iterative process SPICE 
applies a high forward bias voltage to a diode or PN junction, a non¬ 
convergence condition may occur. But if the series resistance terms 
nonzero, under high forward bias conditions, the series resistance 



Figure 3.14 Large conductance 
values (small resistance values) 
lead to small iteration-to-itera- 
tion voltage changes. 
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dominates the conductance of the device and helps reduce the occur¬ 
rence of nonconvergence. 

The series resistance terms are found on the .MODEL parameter 
statements for each of the semiconductor devices and are not global 
parameters found on the .OPTION statement. Look at the model 
parameters of your circuit. A good model should always have the series 
resistance defined. If the model does not have a resistance term in the 
parameter list, select a value of series resistance which is small enough 
to remain unobtrusive to the operation of the circuit and add it to the 
model parameters. Table 3.1 illustrates the series resistance terms for 
each of the semiconductor devices. 

General Newton-Raphson 
Convergence Aids Summary 

The General Newton-Raphson Convergence Aids include learning to 
set the error tolerances for the circuit, choosing a value of GMIN which 
reflects the largest resistances in the circuit, and setting the series 
resistance terms for diodes and bipolar transistors to a nonzero value. 
By following these guidelines, simulations will run faster and fail to 
converge less often. Table 3.2 illustrates the General Newton-Raphson 
Convergence Aids. 

Analysis-Specific Convergence Aids 

Specifying the General Newton-Raphson Convergence Aids will reduce 
the occurrence of nonconvergence. But many nonconvergence failures 
are tied to one or more analysis types. In many of these cases, the Gen¬ 
eral Newton-Raphson Convergence Aids will not be enough to eliminate 
the nonconvergence. In most cases of an analysis-specific nonconver¬ 
gence failure, the nonconvergence may be corrected by understanding 
what caused the problem and setting the appropriate .OPTION state¬ 
ments. The key to overcoming analysis-specific nonconvergence failures 
is to understand what caused the failure, and understanding the cause 
of nonconvergence is where most engineers have trouble. 

SPICE is capable of many different analysis types, but only three of 
these may fail from nonconvergence. SPICE may fail to converge dur- 


table 3.1 Series Resistance Model Parameters 


Semiconductor type 

Series resistance parameters 

Diode 

RS 

Bipolar transistor 

RE and RC 

JFET 

RD and RS 

MOSFET 

RD and RS 





Nonconvergence 


59 


table 3.2 General Newton-Raphson Convergence Aids 


Convergence aid 

Type of circuits 

Error tolerances 

Discrete circuits 

(.OPTION RELTOL=X 

Hybrid circuits 

+VNTOL=X ABSTOL=X) 

Board-level circuits 

Power circuits 
Power-integrated circuits 

GMIN 

(.OPTIONS GMIN=X) 

All circuits 

Series resistance model parameters 
(.MODEL DMOD D RS=X 
.MODEL QMOD NPN RE=X RC=X) 

Diode and bipolar circuits 


ing the DC operating point calculation, during a DC sweep analysis, or 
during a transient time sweep. In overcoming analysis-specific non¬ 
convergence, the sequence of events SPICE follows for a given analysis 
must be known. Table 3.3 illustrates the analysis flow for the four 
major SPICE analysis types. 

Every analysis type in SPICE starts with a DC bias point calculation. 
For this reason, the authors of SPICE added more DC bias point non¬ 
convergence aids than any other analysis-specific nonconvergence aids. 


DC Bias Point Convergence Aids 

The task of determining the bias point of a circuit is both the most dif¬ 
ficult and the most important. The bias point is important because it 
serves at the initial circuit condition for all other analysis types. The 
bias point calculation is difficult because often SPICE has little or no 
information on how the circuit should be biased. (Compare this with 
the DC sweep analysis or transient time sweep where the prior solu¬ 
tion voltages are used as the initial guess for the next series of Newton 
iterations.) 


TABLE 3.3 Analysis Flow in SPICE 


Analysis type 

Analysis flow 

DC operating point 

DC bias calculation 

DC sweep analysis 

DC bias calculation 

DC sweep calculation 

(Use bias point as first point in DC sweep.) 

AC frequency sweep 

DC bias calculation 

AC frequency sweep 

(Use the bias point to determine the linear small-signal 
models of all nonlinear devices.) 

Transient time sweep 

DC bias calculation 

Transient time sweep 

(Use the bias point to determine the time T- 0 value of 
capacitor voltage and inductor current.) 
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Functionally, the DC bias point calculation has already been covered 
both in the diode circuit example earlier in this chapter and in Chap. 2. 
But the specifics of how the DC bias calculation is performed and the 
impact of the .NODESET statement and the ITL1 option have yet to be 
shown. 

DC bias point calculation 

As shown in Chap. 2, SPICE constructs the system (matrix) equations 
from the elements found in the input file. Once the system equations 
are in place, SPICE needs an initial guess for the circuit node voltages 
to start the Newton-Raphson algorithm. 

SPICE makes the initial node voltage array guess by setting any 
nodes connected to a voltage source to the time zero or DC level 
described in the input file. The nodes in the current array which rep¬ 
resent connects to a circuit current source are set to the appropriate 
time zero or DC current level. All the remaining nodes in the circuit 
are set to zero, and this forms the initial guess to start the Newton- 
Raphson algorithm. 

Once the initial guess is entered in the voltage array, SPICE 
begins the iterative process. The Newton iterations continue until 
either all of the nodes and all of the nonlinear branch currents 
converge to within the specified error tolerances or the number of 
iterations exceeds ITL1. ITL1 is the .OPTIONS statement parameter 
which determines the maximum number of iterations SPICE can 
use to determine the DC operating point of a circuit. If all the nodes 
and branch currents have not converged to within the specified error 
tolerances before ITL1 iterations have been used, SPICE aborts the 
iterative process and prints the dreaded “NO CONVERGENCE IN 
DC OPERATING POINT” error message. ITL1 defaults to 100 in 
SPICE. 

If SPICE fails to converge to a DC bias solution, the cause can most 
often be traced to either lack of a good initial guess to start the 
Newton-Raphson iterations, or the lack of enough iterations to reach a 
satisfactory solution for the specified error tolerances. 

Raising ITL1 

The first DC bias point nonconvergence aid is ITL1. Many cases of 
nonconvergence can be eliminated by simply raising ITL1 on the 
.OPTIONS statement. The default setting of ITL1 was determined in 
the mid-1970s when simulating more than a few dozen nodes was rea¬ 
son for a coffee break. Today circuits with hundreds or thousands of 
nodes are common. Many circuits will require more than 100 iterations 
to reach a stable bias point. The first rule of the DC bias point noncon - 
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mergence aids is to raise ITL1. ITL1 sets the limit on the number of iter¬ 
ations SPICE uses before the search for the bias point is aborted. By 
raising ITL1, the extra iterations will only be used by SPICE if they 
are needed. The .OPTIONS statement to raise ITL1 is shown here. 

.OPTIONS ITL1=500 

Table 3.4 shows the number of iterations and the percentage of cir¬ 
cuits which converged to a stable DC bias point with no other conver¬ 
gence aids added. 1 

From Table 3.4, setting ITL1 to 500 yielded the best results. In the 
circuits tested, no measurable improvement was observed by setting 
ITL1 higher. For all circuits, set ITL1 to 500 on the .OPTIONS state¬ 
ment. Some circuits may require more, most less. ITL1 is the first of 
the DC operating point nonconvergence aids to set. 

As an example, simulate the ch3-15a.cir circuit file with the command: 

SIM CH3-15A.CIR 

After 100 iterations, the simulation will terminate with the warning 
“DC Operating Point Nonconvergence.” Then simulate ch3-15b.cir. The 
ch3-15a and ch3-15b circuits are identical except for the .OPTIONS 
ITL1=500 statement. The ch3-15b circuit allows SPICE to use the 
extra iterations required to find the correct circuit bias. 

Setting initial node voltages 

The DIODE.EXE program from Chap. 2 illustrates the importance of 
starting the Newton-Raphson algorithm with an initial guess which is 
close to the final value. Setting the initial node voltages reduces the 
number of Newton-Raphson iterations required to reach a solution. In 
SPICE, circuit nodes can be set to an initial value when starting an 
analysis with the .NODESET statement. Below is an example of using 
the .NODESET statement to set several circuit node voltages. 

.NODESET V(2)=12.5 V(5)=5.0 V(7)=16.8 


table 3.4 


Number of iterations 

Percentage of circuits 

to converge 

which converged 

100 

60 

200 

75 

500 

92 

1000 

92 
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The .NODESET statement forces the designated nodes to a specific 
voltage by applying a Norton-equivalent voltage source to the specified 
node. The Norton-equivalent source is shown in Fig. 3.15. The source 
contains a 1-ohm resistor and a current source. The current source’s 
current is set to the voltage specified on the .NODESET statement for 
that node. For most circuits, the load seen at the specified node is much 
larger than 1 ohm. Under these conditions, all of the current from the 
generator is forced through the resistor, thereby forcing the node to the 
proper voltage. 

When SPICE encounters one or more .NODESET statements, the 
simulator performs two DC bias point calculations. The first calcula¬ 
tion is done while the Norton-equivalent sources are held in the circuit. 
If the circuit converges with the Norton-equivalent sources in place, 
the bias point voltages are saved, the Norton-equivalent sources are 
removed, and a second series of Newton-Raphson iterations deter¬ 
mines the final bias solution. The second bias point calculation is done 
to ensure the Norton-equivalent sources do not impose a load on the 
circuit. 

Circuits of moderate complexity with a single, stable bias point 
rarely require the use of the .NODESET statement. Often, very large 
circuits will converge faster if several circuit nodes are “.NODESET” to 
a specified value. But the best use of the .NODESET statement is on 
circuits with more than one stable operating point. For example, flip- 
flops and latches may start in either the high or low state. Circuits 
with hysteresis may start in one of two states. For these types of cir¬ 
cuits, the .NODESET statement can be used to enhance convergence 
and ensure the circuit starts in the proper state. 

Often in cases of DC bias point nonconvergence, the number of nodes 
which fail to converge is small, usually one or two. When SPICE aborts 
the DC bias calculation, the last iterative node voltage values will be 
printed in the output file. Inspecting these node values often uncovers 


To .NODESET node 



Figure 3.15 The Norton-equivalent voltage source used in SPICE to 
model .NODESET voltages. 
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just one or two nodes which are far from the proper bias values. These 
are the nodes causing the nonconvergence problem, and setting these 
nodes with the .NODESET statement often results in successfully 
locating the bias point in the next simulation. 

A word of caution about the .NODESET statement: Do not try to 
apply the .NODESET statement to circuits which do not have a stable 
DC bias point (such as oscillators and other unstable circuits). When 
the .NODESET statement is applied to unstable circuits, although 
SPICE may converge while the Norton-equivalent source is attached to 
the circuit, during the second series of iterations, after the source has 
been removed, the instabilities of the circuit will dislocate the stable 
condition which existed and cause the simulator to fail to converge. 
SPICE does have a mechanism to start unstable circuits, the .IC state¬ 
ment. But discussion of the .IC statement will be delayed until the sec¬ 
tion on transient analysis. 

Source stepping 

In most cases, DC operating point nonconvergence can be eliminated 
by increasing ITL1 and proper use of the .NODESET statement. But 
some circuits will prove resistant to even these techniques, and some 
circuits will contain subcircuits. Nodes within subcircuits cannot be set 
with a .NODESET statement. For these types of circuits, the Source 
Stepping algorithm may be used to compute the DC bias. 

Source Stepping is a technique which steps the circuit power sources 
from zero (where the solution to all the nodes is zero) to full power. At 
each step, the previous node voltages are used as the initial guess for 
the next series of Newton-Raphson iterations. 

The Source Stepping algorithm is used to calculate the DC bias point 
when the .OPTION ITL6=X is set to a nonzero value. ITL6 is the 
parameter which determines the number of iterations allowed at each 
step in the progression. ITL6 is to Source Stepping as ITL1 is to the 
normal DC bias point calculation. For a circuit which will not converge 
after raising ITL1 and/or after adding the .NODESET statement, set 
ITL6 to 500 or more on the .OPTIONS statement. When ITL6 is set to 
a nonzero value, the Source Stepping algorithm replaces the normal 
DC bias point calculation. 

To demonstrate the Source Stepping algorithm, simulate the ch3- 
18a.cir disk file with the command: 

SIM CH3-18A.CIR 

This large circuit will fail to converge to a DC bias point. Due to the 
size and complexity of the circuit, appropriate .NODESET values can¬ 
not bp, dptprminpd miir.klv. Bv renlaciner the .OPTIONS ITL1=500 with 
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.OPTIONS ITL6=500 statement, the circuit bias point is quickly 
located with the Source Stepping algorithm. Simulate the corrected 
circuit with the command: 

SIM CH3-18B.CIR 

At this point, readers might be tempted to use the Source Stepping 
algorithm in place of the normal DC bias point algorithm for all cir¬ 
cuits. While the Source Stepping algorithm works on many nonconver- 
gent circuits, a flaw within the SPICE implementation of the 
algorithm prevents the Source Stepping method from converging on 
other circuits. 

The Source Stepping algorithm actually begins by ramping the cir¬ 
cuit sources down from full power in a binary descending order (full 
power, Yi power, H power, % power, Me power). During the ramp-down 
phase, SPICE is searching for a convergent bias point to serve as the 
starting point of the ramp-up phase. When a bias point has been found, 
SPICE ramps the sources up in the same binary order. If, during the 
ramp-up phase, a nonconvergence occurs, the algorithm repeats the 
ramp-down, ramp-up phase to work through the nonconvergence. But 
because of a flaw in the algorithm, during a second ramp-down phase, 
the power will be held at a local minimum value and will not be 
ramped below this point. Unable to ramp the sources down further, 
SPICE quickly exceeds the allowed number of iterations for that step 
and aborts the DC bias point calculation. 

For most circuits, raising ITL1 and using the .NODESET statement 
properly will produce an accurate bias point solution. When these tech¬ 
niques fail, the Source Stepping algorithm is a good alternate solution 
technique. 


Turning active elements off 

SPICE contains one final DC bias point convergence aid known as the 
OFF statement. The OFF statement, like the .NODESET statement, 
forces the simulator to compute two DC bias point calculations. During 
the first calculation, one or more active devices may be turned off. Dur¬ 
ing the second calculation, the active devices are turned on, and the 
bias point from the previous calculation is used as the initial guess for 
the second set of Newton iterations. 

The OFF statement is added to one or more semiconductor devices in 
the circuit. Figure 3.16 illustrates the use of the OFF statement for the 
four semiconductor types in SPICE. 

A note of caution about the OFF statement: For many nonlinear cir¬ 
cuits, the circuit bias point with one or more active devices turned off 
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Semiconductor 

Netlist statement 

Diode 

Dxxx na nc D_model OFF 

Bipolar transistor 

Qxxx nc nb ne Q_model OFF 

JFET transistor 

Jxxx nd ng ns J_model OFF 

MOSFET transistor 

Mxxx nd ng ns nb M_model OFF 


Figure 3.16 OFF statement syntax for SPICE semiconductor devices. 


is substantially different from the bias point when the devices are 
turned on. For this reason, the OFF statement rarely achieves conver¬ 
gence on a circuit where the ITL1, the .NODESET statement , and the 
ITL6 techniques have failed. 

DC bias point convergence aids summary 

DC bias point nonconvergence problems can be significantly reduced 
by following the procedures outlined in this section. Finding the DC 
bias point is crucial to simulation because the bias point calculation 
serves as the starting condition for every other analysis in SPICE. For 
this reason, learning to overcome DC bias point nonconvergence is 
paramount in producing accurate, high-quality simulation results. 
Table 3.5 lists the DC bias point convergence aids. Follow the conver¬ 
gence aids in the order shown in Table 3.5 if nonconvergence persists 
after setting the General Newton-Raphson Convergence Aids. 

DC Sweep Convergence Aids 

In many ways, performing the DC sweep analysis is like performing a 
series of DC bias point calculations. Tb begin the analysis, SPICE per¬ 
forms a DC bias point calculation. Once the bias point has been found, 


TABLE 3.5 DC Bias Point Convergence Aids 


Convergence aid 

Order 

Raise ITL1 

First 

(.OPTION ITL1=500 or more) 


Set .NODESET statements 

Second 

(.NODESET V(3)=2.1 V(4)=7.9) 


Use Source Stepping 

Third 

(.OPTION ITL6=500) 


Use OFF statement 

Fourth 

(Dxxx n+ n- d_model OFF) 
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the node voltages are saved and used as the initial guess for the New¬ 
ton iterations at the next point in the analysis. At each solution point 
in the analysis, SPICE uses the previous solution as the initial guess 
for the next series of iterations. This process continues for every point 
in the sweep analysis. When the analysis is complete, SPICE prints 
the output result. 

But the nonconvergence problems associated with a DC sweep anal¬ 
ysis are surprisingly different than the problems of the DC bias point. 
The two nonconvergence mechanisms of the DC sweep analysis are 
rapid voltage transitions and model discontinuities. Either of these 
might cause the simulator to fail to find a solution. 


Model discontinuities 

The semiconductor device models SPICE uses for simulation were pat¬ 
terned after the physical behavior of real devices. Figure 3.17 illus¬ 
trates the family of curves from an MOS transistor. Classical physics 
splits the transistor curves into the linear region and the saturated 
region of operation. The device equations in SPICE follow these same 
regions of operation. But, unlike the real device, SPICE uses separate 
equations for each region of operation. Because of the mathematical 
difficulty in writing an equation which describes the entire family of 
curves, two different sets of equations were written, one for the linear 
region and one for the saturated region, and joined together. Unfortu¬ 
nately, because of the way the equations were joined, there is a discon¬ 
tinuity in the conductance characteristics (the slope of the I-V curve) of 
the device. 2 In Fig. 3.18, the discontinuity exists at the intersection of 
the linear and saturation regions. 



Figure 3.17 The linear and saturated regions of operation of 
an MOS field-effect transistor. 
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Figure 3.18 Model discontinuities of the Level 3 MOS field-effect SPICE transistor 
model. 


A discontinuity in conductance may lead to problems for the Newton- 
Raphson algorithm. Figure 3.19a illustrates the conductance vs. volt¬ 
age characteristics around one of the model discontinuities. On the 
first Newton iteration close to the discontinuity, the conductance value 
leads to a new iterative voltage value on the other side of the disconti¬ 
nuity. In Fig. 3.196, the next Newton iteration falls on a conductance 
value which predicts a solution back on the original side of the discon¬ 
tinuity. In Fig. 3.19c, the third Newton iteration again predicts a solu¬ 
tion on the far side of the discontinuity. When SPICE steps close to or 
on top of a model discontinuity, the Newton-Raphson iterations may 
begin to oscillate around the discontinuity. These oscillations use up 
iterations without progressing towards a solution. 

Although model discontinuities may cause a problem when search¬ 
ing for the DC bias point solution, model discontinuity nonconvergence 
becomes more of a problem during a sweep analysis. While performing 
the DC bias point calculation, only a single solution point is sought. 
During the DC sweep analysis, many solution points are required and, 
in most analyses, the nonlinear circuit elements are being swept 
through one or more regions of operation. Model discontinuities only 
pose a problem if the solution steps align with, or fall very close to, the 
discontinuity. The more steps taken during the sweep analysis, the bet¬ 
ter the change of stepping into or close to a model discontinuity. This is 
why model discontinuity nonconvergence is more of a problem during a 
sweep analysis than during the DC bias point calculation. 






























Figure 3.19(a) and ( b ) Iterations 
around a model discontinuity 
may jump from one side of the 
discontinuity to the other. 


Rapid voltage transitions 

The second nonconvergence mechanism found during a DC sweep 
analysis is caused by rapid voltage or current transitions. Figure 3.20a 

al transfer cha 
haracterist 
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Figure 3.19(c) Model discontinu¬ 
ities may cause oscillation dur¬ 
ing the iterative process. Such 
oscillation will cause SPICE to 
use iterations without approach¬ 
ing a valid solution. 


To start the analysis, SPICE determines the DC bias point of the cir¬ 
cuit. Once the bias point has been found, the input voltage is raised to 
the first voltage step of the .DC statement, and the bias point voltages 
are used as the initial guess for the next series of Newton iterations. In 
any SPICE sweep type analysis, the previous solution voltages are used 
as the initial guess for the next set of Newton iterations . Using the bias 



figure 3.20(a) During a DC sweep analysis, SPICE uses the 
Previous voltage solution as an initial guess for the next set 
of Newton iterations. 
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Figure 3.20(d) During large voltage transitions, a new solu¬ 
tion point may be far from the previous solution point. 
Under these conditions, SPICE often fails to converge. 


voltages as the initial guess, SPICE iterates until a solution for the sec¬ 
ond sweep point has been found or until ITL2 iterations have been 
used. ITL2 is the number of iterations allowed at each step in the DC 
sweep analysis. 

This process repeats at each solution point in the sweep analysis. But 
during large voltage transitions, the voltage change between sweep 
points may be too large for the Newton iterations to find a solution. 
Nonconvergence at a transition point is a common DC sweep failure. 

When nonconvergence occurs at a voltage transition, many engi¬ 
neers reduce the step size of the .DC statement (Fig. 3.20c). While 
reducing the step size of the .DC statement will reduce the size of the 
voltage transition between steps and usually result in a successful sim¬ 
ulation, reducing the step size is not the optimum way to eliminate DC 
sweep voltage-transition nonconvergence failures. 

By reducing the .DC statement step increment, many more solution 
points will be required. Extra solution points require extra iterations 
and lead to much longer simulation runs. And as shown in Fig. 3.21a 
and b, reducing the step size increases the chance of stepping into or 
very close to a model discontinuity and failing to converge because of 
the resulting oscillation. 

Rapid voltage-transition discontinuities are best overcome by rais¬ 
ing the number of iterations allowed at each step in the analysis with 
ITL2. ITL2 limits the number of iterations allowed at each solution 
point in a DC sweep analysis. By default, ITL2 is set to 50. For many 
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Figure 3.20(c) Reducing the DC sweep step size is one 
way of resolving voltage-transition nonconvergence 
problems. But, reducing the step size increases simu¬ 
lation time and increases the chance of encountering a 
model discontinuity. Raising the DC sweep iteration 
limit (ITL2) is the optimum way of resolving voltage- 
transition nonconvergence problems. 


circuits, more than 50 iterations are required to work through large 
rapid voltage transitions. During rapid voltage transitions, 100, 200, or 
more iterations may be required to locate the solution. During a DC 
sweep analysis, for circuits which contain one or more steep transi¬ 
tions, raise ITL2 to 200 on the .OPTIONS statement (.OPTIONS 
ITL2=200). 

If, during a DC sweep analysis, raising ITL2 does not correct the 
nonconvergence problem, the nonconvergence is probably being caused 
by one or more model discontinuities. In this case, users have a choice 
of either obtaining a new model with parameter values which mini¬ 
mize the discontinuity, or increasing or offsetting the analysis step 
size. Often, model discontinuities may be overcome by increasing or 
offsetting the analysis step size so the steps are far enough away from 
the discontinuity to avoid oscillation. Figure 3.22 illustrates how users 
*&ay increase or offset the DC sweep analysis step size. 

DC sweep convergence summary 

DC sweep analysis nonconvergence failures are usually caused by one 
of two failure mechanisms, rapid voltage transitions or device model 
discontinuities. Overcoming these problems is summarized in Table 3.6. 
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Figure 3.21(a) and ( b ) Decreasing the DC sweep step size increases the chance of stepping close to 
a model discontinuity and failing to converge. The smaller the step made, the larger the chance of 
encountering a discontinuity. 


AC Frequency Sweep Convergence Aids 

The AC frequency sweep analysis is less prone to nonconvergence 
problems than any other analysis type in SPICE. The reason for this 
(as explained in Chap. 2) is that the AC frequency sweep is a linear 
small-signal analysis which does not include any nonlinear behavior 
after the bias point has been found. SPICE first determines the DC 
bias point of the circuit; during the bias-point calculation, SPICE may 
fail to converge. In these cases, the DC bias point convergence aids 
may be employed to achieve convergence. After the bias point has 
been found, the nonlinear device models are all replaced with their 
linear small-signal equivalents. With no nonlinear devices in the cir¬ 
cuit, SPICE switches to the simpler LU decomposition (Gaussian- 
elimination) solution algorithm. With the LU decomposition solver, no 
Newton iterations are required to solve the circuit and, because of 
this, nonconvergence disappears. Once the AC frequency sweep anal¬ 
ysis has determined the DC bias point of the circuit, the frequency 
sweep analysis will always converge. The only nonconvergence prob¬ 
lems associated with the AC frequency sweep are problems which 
arise in finding the DC bias point of the circuit, and these problems 
have already been discussed in the section titled “DC Bias Point Con¬ 
vergence Aids,” 
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Figure 3.22 Changing the analysis step size to overcome nonconvergence. 
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table 3.6 DC Sweep Convergence Aids 

Convergence aid 

Cause/order 

Raise ITL1 (.OPTIONS ITL2=200) 

Increase or offset analysis step size (See Fig. 3.24.) 
Develop new model parameter set 

Rapid voltage transitions/First 
Model discontinuity/First 
Model discontinuity/Second 


Transient Convergence Aids 

For many SPICE users, transient analysis is the analysis type most 
often used. In many aspects, transient nonconvergence problems are 
similar to the DC sweep analysis nonconvergence problems, but in 
some ways the two are quite different. 

To begin the transient analysis, SPICE open-circuits the circuit capac¬ 
itors, short-circuits the circuit inductors, and computes the DC bias 
point of the circuit. Immediately after the bias point has been found, the 
capacitors and inductors are restored to the circuit and immediately 
assume the voltages established by the DC bias point calculation (no 
charge-up time is required). After the charge-storage elements are 
placed in the circuit, the first timepoint is calculated, and a series of 
Newton iterations begins computing the solution voltages and currents 
at the first timepoint of the analysis. At each new solution point (time- 
point), a new set of Newton iterations is used to compute the node volt¬ 
ages and branch currents of the circuit. As the solution algorithm 
iterates, the numeric integration algorithms determine the capacitor 
currents and inductor voltages as a function of time. All of this occurs at 
each timepoint after the DC bias of the circuit has been found. 

Dynamic timestep control 

Like the DC sweep analysis, nonconvergence problems primarily arise 
from rapid voltage transitions and device model discontinuities. But, 
unlike the DC sweep analysis, SPICE uses a dynamic timestep control 
algorithm to compute the solution points rather than being specified by 
the user. This comes as a surprise to many users. The analysis state¬ 
ment 

•TRAN INS 100NS 

is an instruction which directs SPICE to perform a transient analysis 
from time T=0 to T= lOOnS. But the first parameter of the statement, 
InS, is known as the print interval and is an instruction to print the 
output results every InS of the analysis. SPICE does not solve the cir¬ 
cuit every InS of the analysis. The internal timestep control algorithm 
Las the job of selecting the actual solution points for the analysis. 
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During transient analysis, SPICE saves all the solution points in 
memory. After the analysis is complete, SPICE uses a linear interpola¬ 
tion routine to compute the evenly spaced output points. Figure 3.23 
illustrates the difference between the internal (solution) timepoints 
and the printed results. (The RSPICE option NOINTR forces RSPICE 
to print these internal timepoints rather than the interpolated time- 
points.) 

The dynamic timestep control algorithm changes the size of the 
timestep depending on the circuit activity. During periods of high cir- 
cuit activity, and large voltage and current changes, the timestep con¬ 
trol algorithm maintains a small timestep to help increase accuracy 
and reduce nonconvergence due to rapid voltage transitions. During 
periods of low activity, or little or no voltage and current change, the 
timestep control algorithm increases the step size to speed the simula¬ 
tion to a conclusion. 

The dynamic timestep control algorithm also enables the transient 
analysis to do something no other analysis can do. If SPICE ever fails 
to converge at a timepoint, the timestep control algorithm automati¬ 
cally discards the nonconvergent timepoint, reduces the timestep to % 
of the original size, and reattempts the solution with a smaller 
timestep. If SPICE fails a second time, the timestep is cut by a factor of 
eight again, and again, and again, until either the circuit converges or 
the timestep is reduced below the allowed minimum for the circuit. 
(The minimum internal timestep is fixed internally and cannot be 
adjusted by the user.) If the internal timestep is reduced below the 
allowed minimum, SPICE prints the “Internal Timestep Too Small’ 5 



Figure 3.23 The timepoints SPICE uses during transient analysis are not 
the same as the points printed in the output file. Many simulators offer an 
option so users may view the actual timepoints SPICE used to solve the 
circuit. 
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warning message (Pspice users: Pspice prints “No Convergence During 
Transient Analysis” in place of this message.) and aborts the simula¬ 
tion. Every other analysis type in SPICE aborts the simulation at the 
first nonconvergent solution point. 

The authors of SPICE recognized that rapid voltage-transition non¬ 
convergence would be a problem during transient analysis. For this rea¬ 
son, the dynamic timestep control algorithm was implemented in SPICE. 


Initial conditions 

Transient analysis is also different from every other analysis type 
because of a new initialization technique, the .IC statement. The .IC 
statement is used to set the initial conditions of the circuit during tran¬ 
sient analysis. The syntax of the .IC statement is shown below. 

.IC V(3)=2.4 V(5)=7.6 V(13)=2.8 

The .IC statement is similar to the .NODESET statement discussed 
earlier in this chapter. Both the .IC and the .NODESET are modeled 
with a Norton-equivalent voltage source. Figure 3.24 illustrates 
the Norton-equivalent circuit of the .IC statement. But this is where 
the similarities between the .NODESET statement and the .IC state¬ 
ment end. 

The .IC statement is only valid during transient analysis, whereas 
the .NODESET statement is valid throughout all analysis types. Dur¬ 
ing the DC operating point, DC sweep, and AC frequency sweep analy¬ 
ses, SPICE disregards any .IC statements in the input file. If, however, 
during a transient analysis, both a .NODESET statement and an .IC 
statement are included in the input file, only the .IC statement will be 
used; the .NODESET statement will be discarded. 


To .IC node 



figure 3.24 The Norton-equivalent voltage source used 
to model .IC statement voltages in SPICE. 
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Another difference between the .NODESET and .IC statements 
involves how SPICE determines the bias point. During a DC analysis, 
the .NODESET Norton-equivalent power supplies are held in the cir¬ 
cuit until a bias point is found. SPICE then removes the power supplies 
and performs a second DC bias calculation so the Norton-equivalent 
supplies do not load the actual circuit. Unlike the .NODESET state¬ 
ment, the .IC Norton-equivalent power supplies are held in for only one 
series of iterations. When a stable bias point has been found, SPICE 
uses the bias point as the initial transient condition of the circuit. At 
times, this may cause problems for some circuits because the bias point 
is established with the Norton-equivalent source still in the circuit. 

The .IC statement is often required for setting the initial conditions 
for charged circuit nodes or for starting unstable circuits. Figure 3.25 
is a simple RC circuit with a IV voltage initial condition across the 
capacitor. To simulate the RC decay of the output voltage, a transient 
analysis might be performed. Simulate the disk file ch3-26a.cir with 
the command: 

SIM CH3-26A.CIR 

Your results should match Fig. 3.26a. In the ch3-26a.cir circuit file, a 
.NODESET statement is used to set the initial IV across the capacitor. 
The problem with the .NODESET statement is the Norton-equivalent 
voltage source is removed after a stable bias point has been found. In 
the ch3-26a.cir, SPICE quickly establishes IV as the initial DC condi¬ 
tion. But after a stable bias point is found, SPICE removes the .NODE¬ 
SET voltage source and performs a second series of iterations. In the 
ch3-26a.cir circuit file, without the Norton-equivalent voltage source, 
the voltage across the resistor immediately falls to zero. (Remember, 
the capacitor is not reintroduced into the circuit until after the bias 
point has been determined.) Because of the need to initialize charged 



Figure 3.25 A simple circuit to demonstrate the dif¬ 
ference between the .NODESET and .IC initial con¬ 
dition statements. 
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circuit nodes, the authors of SPICE introduced the .IC statement. Sim¬ 
ulate the ch3-26b.cir circuit file with the command: 

SIM CH3-26B.CIR 

Your results should match Fig. 3.2 6b. 

The circuit file ch3-26b.eir is identical to ch3-26a.cir except with an 
.IC statement in place of the .NODESET statement in the original file. 

Nonconvergence in transient analysis 

With the differences between the DC sweep analysis and the transient 
analysis defined, the nonconvergence problems of transient analysis 
can be examined. The two dominant nonconvergence problems associ¬ 
ated with transient analysis are rapid voltage transitions and device 
model discontinuities. These are the same two failure mechanisms of 
the DC sweep analysis. 

During a DC sweep analysis, rapid voltage-transition nonconver¬ 
gence was corrected by increasing the number of iterations allowed at 
each step in the analysis or by reducing the step size of the analysis. 
In transient analysis, the dynamic timestep control algorithm auto¬ 
matically reduces the step size during large circuit transitions 
automatically. In a DC sweep analysis, device model discontinuity 
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Figure 3.26(a) The output response of Fig. 3.25 circuit with a .NODESET statement used 
to set the initial capacitor voltage. ( Reprinted from Successfully Simulating Circuits 
With SPICE. Used with permission. ) 
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Figure 3.26 (b) The output response of Fig. 3.25 circuit with a .IC statement used to set 
the initial voltage. {Reprinted from Successfully Simulating Circuits with SPICE. Used 

with ru>rmissinn 'I 


nonconvergence was overcome by increasing the step size but, during 
transient analysis, the dynamic timestep control algorithm deter¬ 
mines the step size. 

During transient analysis, as the circuit approaches a voltage tran¬ 
sition, two potentially conflicting events occur. First, the timestep 
control algorithm reduces the size of the timestep, and second, during 
voltage transitions, the semiconductor devices of the circuit change 
state and transit from one region of operation into another. As the 
semiconductor devices change between regions of operation, one or 
more model discontinuities may be exposed. A voltage transition 
during transient analysis causes the timestep to be reduced during 
the same period that device model discontinuities are most likely to 
be uncovered. Worse yet, if a discontinuity causes a nonconvergence, 
the timestep control algorithm automatically reduces the size ol 
the timestep, and the reduced timestep is almost guaranteed to 
stumble into the same discontinuity and fail to converge again. (In 
the section on DC sweep nonconvergence, model discontinuities were 
best overcome by increasing the step size to help step over the dis¬ 
continuity.) 

Device model capacitance 

Transient nonconvergence is primarily caused by a combination ol 
model discontinuities and a significantly reduced step size brought on 
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by voltage transitions within the circuit. This failure mechanism can 
be corrected by first recognizing that the model discontinuities are 
found in the DC I-V characteristics of the device models. All of the 
device models have a built-in capacitance. The built-in capacitances 
are either junction capacitance associated with the PN junction of the 
device or overlap capacitances associated with the insulated gate of a 
MOSFET. Often these capacitive device currents may help bridge the 
DC discontinuity. Unfortunately, the built-in capacitances all have a 
default value of zero, which could lead to convergence problems if a DC 
discontinuity exists. All real semiconductor devices have some capaci¬ 
tive component; therefore, all simulation models should have their 
associated capacitance terms set to a nonzero value. Adding the capaci¬ 
tance terms to a model helps to improve both the accuracy of the tran¬ 
sient analysis and the convergence characteristics of the simulator. 
Table 3.7 illustrates the capacitance model parameters for each of the 
semiconductor devices. 

If a model parameter set does not include the capacitive terms, either 
find an accurate value for the capacitance and add the value to the 
model set, or set the capacitance parameter to a value small enough to 
be inconsequential to the normal operation of the circuit. Even a small 
capacitive component is often enough to overcome nonconvergence 
problems. Table 3.8 illustrates general minimum values for each capac¬ 
itance term for both discrete and integrated circuit components. For 
some circuits these will be too large. Use Table 3.7 with caution. 

A simple SR flip-flop is shown in Fig. 3.27. The flip-flop is con¬ 
structed from CMOS equivalent gates. The circuit file also contains 
two pulse voltage sources, one on the set pin and one on the reset pin. 
The simulation file performs a transient simulation on the flip-flop, 
driving the latch through several state changes. Simulate the ch3- 
27.cir file with the command: 

SIM CH3-27.CIR 

The circuit will fail to converge at the first voltage transition. 

Now simulate the ch3-28a.cir circuit file. This file is identical to the 
ch3-27.cir file except for the addition of the MOSFET capacitance 


TABLE 3.7 Capacitance Model Parameters 


Device type 

Capacitive model parameters 

Diode 

CJO 

Bipolar transistor 

CJE, CJC, CJS 

JFET transistor 

CGD, CGS 

MOSFET transistor* 

CGDO, CGSO, CGDO, CBD, CBS, CJ, CJSW 


* For MOSFET models, use either CBD and CBS or CJ and CJSW. 
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terms CGBO, CGSO, CGDO, CBD, and CBS. The capacitance terms 
were determined from Table 3.7. 

SIM CH3-28A.CIR 

Your result should match Fig. 3.28a. In this circuit, adding a small 
capacitance term eliminates the nonconvergence failure. 

Raising the iteration limit 

In the section on DC sweep analysis, raising the iteration limit allowed 
the simulator more iterations before aborting the simulation. For tran¬ 
sient analysis, the iteration limit is ITL4. If a timepoint fails to con¬ 
verge within ITL4 iterations, SPICE discards the timepoint, cuts the 
timestep by a factor of eight, then reattempts the solution with the 


TABLE 3.8 Minimum Capacitance Parameter Settings 


Device/capacitance parameter 

ICs 

Discrete 

Diode/C JO 

.05pF 

.lpF 

Bipolar/CJE 

.lpF 

.2pF 

Bipolar/C JC 

.lpF 

.2pf 

Bipolar/C JS 

.lpF 

.2pf 

JFET/CGD 

.2pF 

.2pf 

JFET/CGS 

.2pF 

.2pf 

MOSFET/CGDO t 

20pF 

•5 P f 

MOSFET/CGSO f 

20pF 

•5pf 

MOSFET/CGBO 1 

lOpF 

.025pf 

MOSFET/CBD* 

.lpF 

.5pf 

MOSFET/CBS* 

.lpF 

.5pf 

MOSFET/CJ* 

20pF 

N/A 

MOSFET/CJSW* 

lOpF 

N/A 


* For MOSFET models, use either CBD and CBS or CJ and CJSW. 
f The MOSFET capacitances CGDO, CGSO, and CGBO are per meter of 
gate length or per meter of gate width. That is why these values appear so 
much larger than the other capacitive terms. 
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Figure 3.28(a) SR flip-flop transient response with capacitance model parameters 
included in circuit file. C Reprinted from Successfully Simulating Circuits with SPICE. 
Used with permission.) 


smaller timestep. By raising ITL4, SPICE has more iterations to con¬ 
verge on a solution. This reduces the number of times SPICE is forced 
to shrink the timestep; the larger the timestep, the less likely SPICE is 
to run into or close to the model discontinuity and fail to converge. 

By default, ITL4 is set to only 10 iterations. For all transient 
circuits, raise ITL4 to 40 iterations on the .OPTIONS statement 
(.OPTIONS ITL4=40). 

Simulate the ch3-28b.cir circuit file with the command: 

SIM CH3-28B.CIR 

Your results should match Fig. 3.28b. The ch3-28b.cir circuit file con¬ 
tains the original SR flip-flop without the capacitance terms but with 
the ITL4 parameter raised to 40. Notice the result obtained with this 
circuit is identical to the previous result hut the run time for the ch3- 
28b.cir is less than half of the ch3-28a.cir simulation! 

A side benefit to raising ITL4 is an increase in simulation speed. 
Chapter 5 will show that the timestep control algorithm in SPICE 
always cuts the timestep by a factor of 8, yet the timestep control algo¬ 
rithm will increase the timestep by no more than a factor of 2. If SPICE 
fails to converge and cuts the timestep, at least three more timepoints 
must be solved before the simulator returns to the point in the simula¬ 
tion where the timestep was originally cut. For many circuits, raising 
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Figure 3.28 (b) SR flip-flop transient response with transient iteration limit (ITL4) set to 
40. (Reprinted from Successfully Simulating Circuits with SPICE. Used with permission.) 


ITL4 will increase the speed of the simulator and reduce the occur¬ 
rence of nonconvergence. 

Transient convergence aids summary 

Overcoming transient nonconvergence failures requires an under¬ 
standing of DC bias point and DC sweep analysis nonconvergence 
failures. The dynamic timestep control algorithm helps reduce non¬ 
convergence failures due to rapid voltage transitions but may aggra¬ 
vate failure due to model discontinuities. Restoring the device model 
capacitance terms from their default of zero to a more realistic value 
helps minimize model discontinuities. Raising ITL4 gives the timestep 
algorithm more time to converge on a solution before the timestep is 
reduced. Maintaining a larger timestep decreases the chance of step¬ 
ping into or close to a model discontinuity. Used together, these two 
techniques will eliminate most transient nonconvergence failures. 
Table 3.9 outlines the transient convergence aids and the order in 
which they should be applied. 

Nonconvergence Summary 

Most of the common nonconvergence problems users face can easily be 
overcome with the controls and options available in SPICE. Many 
problems occur simply because the error tolerances and program 
defaults are not appropriate for a given circuit (Table 3.2). Once these 
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TABLE 3.9 Transient Convergence Aids 


Convergence aid 

Cause/order 

Add capacitance model parameters 

First 

(.MODEL DIO D (CJ0=.lpf) 


Raise ITL4 

Second 

(.OPTIONS ITL4=40) 



are set, specific analysis types may cause problems. Table 3.5 illus¬ 
trates the convergence aids for computing the DC bias point of the cir¬ 
cuit. Table 3.6 shows the convergence aids for a DC sweep analysis. 
Finally, the transient convergence aids are shown in Table 3.9. 

Set the general convergence aids first. When nonconvergence occurs, 
determine which analysis caused the nonconvergence, then go to the 
table for that analysis type. Do not confuse convergence aids. Setting 
the DC bias point aids will not help resolve transient nonconvergence 
problems. Setting transient aids will not help during DC sweep. Set¬ 
ting a DC sweep aid will not help during an AC analysis. Identify 
which analysis caused the nonconvergence, and learn which conver¬ 
gence aids apply to each analysis type. 

Convergence Aids 

in the SPICE-Like Simulators 

Hspice 

Hspice has many of the same convergence aids as standard SPICE, as 
well as several enhancements. Convergence aid options include GMIN, 
RELTOL, ABSTOL, VNTOL, ITL1, ITL2, and ITL4. The enhancements 
include Source-Stepping, GMIN-Ramping, and Pseudo-Transient bias 
determination. 

Hspice uses a modified Source-Stepping algorithm when the option 
CONVERGE=3 is set. Hspice uses a variable step size algorithm in 
place of the binary stepping algorithm found in SPICE. 

When the Hspice option GRAMP=X is specified, Hspice calls the 
GMIN-Ramping algorithm. GMIN-Ramping involves setting the GMIN 
conductance (resistance) to a large (small) value. With a large conduc¬ 
tance (small resistance) across every nonlinear device in the circuit, 
much of the nonlinear behavior of the circuit is suppressed. Under 
these conditions, HSPICE converges very quickly. After converging, 
Hspice drops the value of GMIN by a factor of 10 and reattempts a 
solution. This process continues until the circuit fails to converge or 
GMIN is reduced to the default setting. If Hspice fails to converge, 
GMIN is increased by a factor of 10 and held at that value for the rest 
of the simulation. 
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The GRAMP=X parameter sets the exponential value of the GMIN 
conductance to start the ramping method. For example, if the .OPTION 
GRAMP=6 was used^Hspice would ramp the GMIN conductance (resis¬ 
tance) from le-6 down to le-12 mhos (le6 ohms up to lel2 ohms). 

To find the DC bias point of the circuit, Hspice added a technique 
known as Pseudo-Transient analysis. With the Pseudo-Transient anal¬ 
ysis, a small capacitor is attached to every node in the circuit. The 
analysis continues by ramping the DC power supplies from zero up 
to full power. Once the power supplies are at full power, the capacitors 
are removed from the circuit, and the result is the proper DC bias 
point. The Hspice option CONVERGENT invokes the Pseudo-Transient 
analysis. 

IS_Spice 

IS_Spice has the same convergence aid as standard SPICE. Conver¬ 
gence aid options include GMIN, RELTOL, ABSTOL, VNTOL, ITL1, 
ITL2, ITL3, ITL4, and ITL6. 

Micro-Cap IV 

Micro-Cap IV has many of the same convergence aids as standard 
SPICE, with the addition of the SAVE BIAS command. Convergence 
aid options include GMIN, RELTOL, ABSTOL, VNTOL, ITL1, ITL2, 
ITL3, and ITL4. 

While not directly a convergence aid, the SAVE BIAS command 
allows users to reuse the bias point voltages from a previous simula¬ 
tion. For large circuits, reusing the bias point may increase the speed 
of subsequent runs by eliminating the time and iterations consumed in 
the bias point calculation. 

Pspice 

Pspice has many of the same convergence aids as standard SPICE, 
with three additional enhancements. Convergence aid options include 
GMIN, RELTOL, ABSTOL, VNTOL, ITL1, ITL2, and ITL4. The con¬ 
vergence enhancements include an additional GMIN resistor across 
the independent current sources and dependent voltage and current 
sources, automatic source-stepping for the DC bias point and DC 
sweep calculations, a reduced resistor for the .NODESET and .IC 
statements, and the SAVE BIAS command. 

A GMIN resistor is placed across each independent current source 
and every dependent source in the circuit. Changing the value of gmin 
on the .OPTIONS statement alters this value. 

If Pspice fails to converge within ITL1 iterations during a DC bias 
point calculation, or ITL2 iterations during a DC sweep analysis, the 
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program automatically switches into the Source-Stepping algorithm. 
The Source-Stepping algorithm in Pspice uses a variable step size 
rather than the binary stepping method of standard SPICE. 

In SPICE, .NODESET and .IC statements are modeled with Norton- 
equivalent source. The Norton-equivalent source contains a current 
source and a 1-ohm resistor. For many low-impedance circuits, the 
1-ohm resistor adds significant loading to the circuit. In Pspice, 
the 1-ohm resistor is replaced with a .002-ohm resistor, and the value 
of the current source is increased to match the reduced resistance. 

While not directly a convergence aid, the SAVEBIAS command 
allows users to reuse the bias point voltages from a previous simula¬ 
tion. For large circuits, reusing the bias point may increase the speed 
of subsequent runs by eliminating the time and iterations consumed in 
the bias point calculation. 

Summary 

Nonconvergence is one of the most common and frustrating problems 
facing simulation users. Resolving nonconvergence problems involves 
identifying and eliminating the cause of the nonconvergence. Set the 
general convergence aids on every circuit you simulate. If nonconver¬ 
gence occurs, identify which analysis type caused the nonconvergence 
and set the analysis-specific convergence aids. Following the guide¬ 
lines outlined here will resolve between 80 and 90 percent of all non¬ 
convergence problems. 

Some nonconvergence problems may be more stubborn than others. 
For tough problems, follow the guidelines exactly as shown. Some¬ 
times, the first convergence aid you try will resolve the problem. Other 
problems will require all of the convergence aids. Do not worry about 
adding too many convergence aids. You will not introduce any appre¬ 
ciable error in your simulation if you follow the guidelines offered in 
this chapter. The guidelines presented here represent a proven, sys¬ 
tematic approach to resolving nonconvergence problems. 

Overcoming convergence problems is possible and practical. Before 
the accuracy of a simulation can be determined, before the speed of a 
simulation can be improved, before the analysis results can be exam¬ 
ined, the circuit must converge. 
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Numeric Integration 


What is numeric integration? Why is numeric integration used in 
SPICE? These are two important questions which any knowledgeable 
simulation user should be able to answer. If you don’t know the 
answers to these questions, read on. Too few users realize the signifi¬ 
cant role that numeric integration plays in producing an accurate tran¬ 
sient simulation. 

Numeric integration is used in SPICE to calculate the current flow¬ 
ing through the circuit capacitors as a function of time; in a dual role, 
numeric integration is used to calculate the voltage across the circuit 
inductors as a function of time. 

As to what numeric integration is, numeric integration is analytic 
integration applied in a piecewise linear manner. Whereas analytic 
integration produces a continuous function which represents the solu¬ 
tion, numeric integration produces discrete solution points. If all the 
solution points are plotted on a graph, the resulting curve will match 
the function generated from the analytic integration. 

Most often, numeric integration is used to solve differential equations 
which have an unknown or difficult-to-calculate analytic solution. In 
the case of circuit simulation, the differential equations of interest are: 


and 


C*dV 

dt 


(4.1) 


L*dl 

dt 


(4.2) 


Numeric integration determines the capacitor currents and inductor 
voltages for a transient simulation by solving the differential equa¬ 
tions relating voltage and current. 
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During transient analysis, SPICE uses one of three user-selectabl 
integration methods—trapezoidal, backward-Euler, or Gear—to detei 
mine the capacitor currents and inductor voltages. In addition to th 
standard Gear integration formula, four additional variations of th 
Gear equation may be selected. This gives the SPICE user a wide var 
ety of integration methods from which to choose. 

Transient Simulation Warning 

A variety of integration methods are provided in SPICE simp] 
because no one integration method is always the best for all transier 
simulations. When a numeric integration formula is applied to a di 
ferential equation, the piecewise linear solution generated by the into 
gration is only an approximation of the exact function. At each solutio 
point, the integration algorithm may introduce a small amount < 
error. The accuracy of the numeric integration method will depend o 
the size of the simulation time step, the shape of the voltage and cu 
rent waveforms, and the integration method being used. Because of tl 
many different waveforms which arise during transient analysis, tl 
numeric integration algorithms must be extremely robust to maintai 
accuracy throughout the variety of circuit responses. A given integr; 
tion method may work well on one function or waveform and produ< 
completely inaccurate results on another. Because no one integratic 
method is best suited for all transient simulations, SPICE users mui 
be able to decide which integration method produces the most accural 
result for their circuit. 

When used without caution, numeric integration inaccuracies can 1 
an Achilles’ heel for the simulation user. Unlike nonconvergence fai 
ures, SPICE will not generate a warning message when the numer 
integration method is producing inaccurate results. The SPICE us< 
must learn to look for signs which indicate something has gone wror 
during the simulation. Often, when the integration method does fa: 
simulation results will display well-defined, well-characterized anor 
alous behavior. A numeric integration failure may appear as pr 
nounced as a catastrophic failure or as subtle as a voltage or currei 
ripple. But to the trained eye, these same signs are the signals 
switch to a different integration method. Learning the characterist 
signs of numeric integration failure, and what to do to correct the fa: 
ure, are important lessons in simulation craftsmanship. 

Warning summary 

Simulation users should not presume the results obtained from the 
simulations are always correct. There are two reasons for this: firs 
model inaccuracies will lead the simulator to inaccurate results; se 
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ond, even when the models are completely accurate, several failure 
mechanisms of the simulator can lead to the wrong solution (numeric 
integration failure is just one of the failure mechanisms) without 
warning the user. Never simply assume the simulation output is cor¬ 
rect . If questionable results appear, use good engineering judgment to 
determine whether the anomaly was simulator-related or design- 
related. Never simulate a circuit without a reasonably good idea of 
what the simulation output should be. A good designer would never 
breadboard a collection of resistors, capacitors, inductors, and transis¬ 
tors without some reasonable expectation of the behavior of the circuit. 
Simulation users should treat computer simulation the same way. 

Numeric Integration 

Before going further, a simple example of numeric integration might be 
useful here. A time function is shown in Eq. 4.3. 

Fit) = 2 1 (4.3) 

The derivative of Eq. 4.3 may easily be computed and is shown in Eq. 
4.4. The derivative F'(t) is a constant and is by definition a differential 
equation. Equation 4.4 will serve as the starting point for our example. 

FXt) = 2 (4.4) 

Equation 4.4 is a differential equation. A differential equation is 
nothing more than the derivative of a function. Often the original func¬ 
tion is unknown, but in this case both the derivative and the original 
function are known. Solving a differential equation involves trans¬ 
forming the differential equation back into the original function. 

Equation 4.4 may be solved with an analytic technique known as 
Separation of Variables. Begin by rewriting Eq. 4.4 in the form of a dif¬ 
ferential equation. 


2 (4.5) 

dt 

Multiplying both sides of Eq. 4.5 by dt results in Eq. 4.6. 

dY(t) = 2 *dt (4.6) 

Equation 4.6 may be integrated over time, yielding Eq. 4.7. 

Y(t) = 2 t + C (4.7) 

Not surprisingly, Eqs. 4.7 and 4.3 are identical except for the con¬ 
stant C. 
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The technique of separation of variables is an analytic method of 
solving a differential equation. Equation 4.5 can also be solved with 
numeric integration. 

Numeric integration produces results similar to the analytic solu¬ 
tion but in a piecewise linear form. Notice the word “similar” was used, 
not “identical.” Numeric integration is an approximation to the solu¬ 
tion. At each solution point, a finite amount of error may be introduced. 
Often, the error will be insignificant and the numeric integration solu¬ 
tion will be very close to the analytic solution, but at times the inte¬ 
gration error can be significant. 

During numeric integration, the original function is reconstructed as 
a series of discrete data points at each step in the integration. To apply 
numeric integration to Eq. 4.5, construct an X-Y coordinate grid as 
shown in Fig. 4.1. Since Eq. 4.4 is a function of time , discrete units of 
time will be the interval of integration. 

To apply numeric integration to a differential equation, an integra¬ 
tion formula must be chosen, and a starting point or initial value of the 
solution must be known. For this example, the backward-Euler 
numeric integration formula will be used, and Y(t = 0) will be set to 
zero. The backward-Euler formula is shown in Eq. 4.8. 

Y(x + 1) = Y(x) + step size * dF( * - + 1} (4.8) 

dx 


Y(t) 


J-1-1-1_I_I_L 


Figure 4.1 X-Y coordinate grid. 


Time 
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Equation 4.9 results when the backward-Euler formula is applied to 
Eq. 4.5. 


Y(t + 1) = Y(t) + timestep * dY{t + 1} (4.9) 

dt 

In Eq. 4.9, Y(t + 1) represents the solution to the differential equa- 
tion at the timepoint T+ 1. Y(t) is the solution at the past timepoint T. 
For our example, the timestep will be set to a constant .5 seconds. The 
integration may be started by setting Y{t = 0) to zero and evaluating 
Eq. 4.9 at discrete .5-second timepoints. Table 4.1 displays the calcula¬ 
tion and the computed value of Y(t) at each timepoint. Figure 4.2 recon¬ 
structs the graph of Y(t). 

Many practical engineering problems must be solved with differen¬ 
tial equations. The goal of numeric integration is to reconstruct the 
original function from the known differential equation. At each step 
in the integration, the integration formula produces a discrete solu¬ 
tion value. The original function is represented by the set of solution 
values. 

In this example, the numeric integration solution was exactly equal 
to the analytic solution, but only because this was a simple function 
where the differential equation is constant. In more complex functions, 
such as those with a nonconstant first or nonzero second derivative, 
many of the numeric integration methods will produce small amounts 
of error at each step in the integration. 

SPICE and Numeric Integration 

In SPICE, the numeric integration routines solve the differential equa¬ 
tions describing the current-voltage relationships for capacitors and 
inductors. Equation 4.10a is the backward-Euler expression for the I-V 


TABLE 4.1 Calculated Timepoint Solution Values 


Y(0) 

= 0 



Y(.5) 

= 0 

+ .5*2 

= 1 

Y(1.0) 

= 1 

+ .5*2 

= 2 

Y(1.5) 

= 2 

+ .5*2 

= 3 

Y(2.0) 

= 3 

+ .5*2 

= 4 

Y(2.5) 

= 4 

+ .5*2 

= 5 

Y(3.0) 

= 5 

+ .5*2 

= 6 

Y(3.5) 

= 6 

+ .5*2 

= 7 

Y(4.0) 

= 7 

+ .5*2 

= 8 

Y(4.5) 

= 8 

+ .5*2 

= 9 

Y(5.0) 

= 9 

+ .5*2 

= 10 
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relationship of a capacitor. A dual expression can be written for the 
inductor. 


VAt + 1) = V c (t) + timestep *JAt + 1) 


(4.10a) 


To determine the capacitor current from Eq. 4.10a, the voltage and 
current terms must appear on opposite sides of the equality. Using sim¬ 
ple algebra to rearrange the terms of Eq. 4.10a results in a more famil¬ 
iar expression for capacitor current. 

[V c (t + D-Vcffl] *c=I c {t + 1) (4.106) 

timestep 

At each timepoint in the transient simulation, the node voltages 
from the present solution point V c (t + 1) and the previous timepoint 
V c (t) are used to determine the change in capacitor voltage. 

[v c (t) - v c (t -1)] = dVc_ (4 n) 

timestep dt 

Once the change in voltage is known, the capacitor current can be 
determined by multiplying Eq. 4.11 by the capacitance value. 

c * mt +1) - vm _ c*dv c - IJf + v 

timestep dt 


( 4 . 12 ) 
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When SPICE begins the transient simulation, an equation similar to 
Eq. 4.12 is generated for every capacitor in the circuit, and a dual equa¬ 
tion is generated for every inductor in the circuit. 

In searching for a solution, SPICE iterates on node voltage values over 
and over again until a set of voltages is found that satisfies KirchhofTs 
voltage and current laws and satisfies Eq. 4.12 for every capacitor/ 
inductor in the circuit. Since V c {t) is the node voltage from the previous 
timepoint (not from the previous iteration) and remains fixed, only 
V c (t + 1) changes with each new iterative voltage value. As V c (t + 1) 
changes, I c (t + 1) changes to match the new capacitor voltage. When the 
exact solution is found, Eq. 4.12 becomes the relationship between the 
capacitor voltage and current at the transient timepoint. 

Types of Numeric Integration 

SPICE offers users three different numeric integration methods for use 
during transient simulation. The three methods are backward-Euler, 
trapezoidal, and Gear integration. Of these, both the backward-Euler 
and trapezoidal methods were developed in the mid-1700s. 1 Only the 
Gear method was developed during this century. 2 The different inte¬ 
gration methods were written into SPICE because Nagel 3 found that 
no one method was reliable under all transient simulations. For some 
circuits the trapezoidal method produces the most accurate result; for 
other circuits, the Gear method works best; and yet other circuits will 
produce the best result with the backward-Euler method. 

Although each integration method must perform the same task, the 
way each method generates the solution will differ. An understanding 
of how each method accomplishes the task of solving the capacitor cur¬ 
rents and inductor voltages, and knowing the difference between 
methods, will help in choosing an appropriate method for a given cir¬ 
cuit type. 

Although SPICE has only three different integration types, four will 
be presented here to illustrate the differences and similarities between 
methods. The four of interest are forward-Euler, backward-Euler, 
trapezoidal, and Gear integration. To illustrate the differences between 
methods, a graphical technique utilizing two X-Y coordinate grids will 
be employed. The Y-axis of the first graph represents the differential 
equation dY(t)/dt. The second graph represents both the exact analytic 
solution and the numeric integration solution. The analytic solution is 
shown in the second graph as a solid line. The numeric integration solu¬ 
tion is shown in the graph as a series of data points. 

Figure 4.3 represents the forward-Euler integration method. The 
forward-Euler estimates the value of the function Y(t + 1) by adding 
the previous solution value Y(t) to the product of the timestep and the 
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derivative of the timepoint T, dY(t)fdt. The name forward-Euler comes 
from the use of the previous timepoint T in establishing the derivative 
for a solution at T + 1. In a sense, the forward-Euler method is esti¬ 
mating the solution based on the derivative of the past solution point. 
The forward-Euler integration formula is shown in Eq. 4.13. 

Y(t + 1) = Y(t) + timestep*<iY(£) (4.13) 




Time 

y(7+1) = Y(T) + h » 

at 

Figure 4.3 Forward-Euler numeric integration formula and graphical interpretation of 

method 
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Figure 4.4 illustrates the backward-Euler integration method. The 
forward- and backward-Euler methods are very similar. The difference 
between the two methods is the selection of the timepoint where the 
derivative of the function is being evaluated. The backward-Euler 
method estimates the value of the function Y{t + 1) by adding the pre¬ 
vious solution value Y(t) to the product of the timestep and the deriva¬ 
tive of the timepoint T + 1, dY(t + 1 )/dt. Whereas the forward-Euler 
method uses the derivative at the past solution point T, the backward- 




Time 

y(r+1) = Y(T) + h * dYi ^' ) - 

Figure 4.4 Backward-Euler numeric integration formula and graphical interpretation of 
method. 
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Euler method uses the derivative of the current solution point T + 1. 
The backward-Euler integration formula is shown in Eq. 4.14. 

Y(t + 1) = Y{t ) + timestep*c?Y(£ + 1) (4.14) 

Figure 4.5 represents the trapezoidal integration method. The trape- 
zoidal is similar to the Euler methods but again differs in the selection 
of where the derivative of the function will be evaluated. Like the 
Euler methods, the trapezoidal method estimates the value of Y{t + 1) 
by adding the previous solution value Y(t) to the product of the 
timestep and an averaged derivative. The trapezoidal algorithm uses 
the derivative of the past solution point T and the derivative at the cur¬ 
rent solution point T + 1. These two derivative values are averaged, 
and this average derivative is used in the calculation. Rather than 
evaluating the derivative at the past timepoint (forward-Euler) or at 
the present timepoint (backward-Euler), the trapezoidal method aver¬ 
ages the two derivatives as shown in Fig. 4.5. Because of this, the 
trapezoidal method tends to be more accurate than either the forward- 
or backward-Euler methods. The trapezoidal integration formula is 
shown in Eq. 4.15. 

Y(t + 1) = Y(t) + timeste P*fo^) + dY (t + (4 15) 

2 

The Gear integration method 2 is a relative newcomer. Figure 4.6 rep¬ 
resents the Gear integration method. The Gear integration method is 
very different from traditional integration methods and estimates the 
value of the function at time T + 1 with information from the present 
Solution point and two past solution points! The Gear-2 integration for¬ 
mula is shown in Eq. 4.16. 

y(( + D = (|)*Ytt) 

+ (| j*timestep*cm« + 1) (4.16) 

In addition to the Gear-2 formula, there are Gear-3, Gear-4, Gear-5, 
and Gear-6 formulas. Where the Gear-2 formula uses two previous 
timepoints, Y(t ) and Y{t - 1), the Gear-3 formula uses three, Y(t ), Y(t - 1), 
and Y{t - 2). This progression continues up to the Gear-6 formula which 
uses six previous timepoints to help estimate where the function will be 
at time T + 1. 
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Time 

y(7+1) = Y{T) + h* dY ( T+1 )J d t 

Figure 4.5 Trapezoidal numeric integration formula and graphical interpretation of 
method. 


Accuracy and Stability 
of Integration Methods 

Due to the difference in the integration formulas, each method will 
produce a different result when applied to a given function. How well a 
given integration formula will predict the correct solution is deter¬ 
mined by the accuracy and stability of the method. 

An integration method is approximating the solution of a differential 
equation at each timepoint in the transient simulation. Because the 



























































98 Chapter Four 



Time 

V(7+1) = f * Y(T) - ~ » Y(T- 1) + | * h * dYi ™ } 

Figure 4.6 Gear method numeric integration formula and graphical interpretation of 
method. 

numeric integration solution is only an approximation to the exact 
solution value, a finite amount of error may be introduced at each time- 
point. The amount of error at a given timepoint is a measure of the 
accuracy of the integration method. Of the four integration methods 
presented (backward-Euler, forward-Euler, trapezoidal, and Gear), 
each will have a different amount of error introduced at a given time- 
point. The error introduced at each timepoint in a numeric integration 
is known as the local truncation error. Figure 4.7 illustrates the accu¬ 
racy of an integration method and the local truncation error for an 
arbitrary function. 
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Figure 4.7 The accuracy of a numeric integration method refers to the 
amount of error at each timepoint of the integration. The error at each 
timepoint is known as Local Truncation Error. 


A second criterion for a numeric integration method is stability. How 
the local truncation error accumulates over a large number of time- 
points is a measure of the stability of an integration method. A stable 
integration method may generate a result which overestimates the 
actual solution at some timepoints and underestimates the actual solu¬ 
tion at others. But over a large number of timepoints (a long transient 
simulation), a stable integration method produces a result which 
closely approximates the actual solution. Conversely, an unstable 
method tends to accumulate the local truncation error at each time- 
point. Over a large number of timepoints, an unstable integration 
method will diverge from the exact solution. Figure 4.8 shows the 
result of a stable and an unstable integration method when applied to 
a general function. 

Both the accuracy and stability of a given integration method will be 
determined by the first- and second-order derivatives of the function 
being integrated and by the timestep used during integration. Because 
of this, no one integration method will always be the most accurate or 
the most stable. An integration method may generate large amounts of 
local truncation error when applied to one function, and generate little 
or no error when applied to a different function. A given integration 
method may be stable on some functions and unstable on others. 

Although this text will set some general guidelines on which inte¬ 
gration method you should use on your circuits, any integration 
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Figure 4.8 The stability of an integration method refers to how the numeric integration 
solution behaves over a large number of timepoints. Accurate numeric integration meth¬ 
ods may not always be stable and stable methods may not be accurate. 


method is prone to accuracy and stability problems. A knowledgeable 
simulation user will learn to recognize numeric integration inaccura¬ 
cies and will take corrective actions. The following sections will show 
you how to look for anomalous behavior and what to do to correct your 
simulation. 

Comparison of Integration Methods 

To illustrate the difference between integration methods, three sample 
time functions were chosen. All three had smooth and continuous 
derivatives. The derivative of each function was calculated and then 
the forward-Euler, backward-Euler, trapezoidal, and Gear-2 integra¬ 
tion methods were applied to the computed derivatives. The resulting 
curves represent the four numeric integration approximations of the 
original function. 

Equations 4.17, 4.18, and 4.19 show the three sample functions used 
for this illustration. 


fit) = t 2 - 5t 

(4.17) 

fit) = 5 - 5 *expi-t) 

(4.18) 

fit) = sin(£) 

(4.19) 
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Figures 4.9a, b, and c represent Eq. 4.17 for a step size of Is, .5s, and 
.25s respectively. By decreasing the step size of the integration, the 
local truncation error decreases. While no one integration method is 
always the most accurate or the most stable, decreasing the step size of 
any integration method improves the accuracy of the solution. Figures 
4.10a, b, and c represent Eq. 4.18, and Figs. 4.11a, b, and c represent 
Eq. 4.19. In each illustration, decreasing the step size produces a pro¬ 
nounced improvement in the accuracy of solution. 

The previous statement deserves careful consideration. By decreas¬ 
ing the step size of the integration, or by decreasing the timestep dur¬ 
ing transient analysis, the resulting numeric integration solution for 
capacitor currents and inductor voltages improves. Many SPICE users 
know that decreasing the timestep cures many transient anomalies. (A 
more accurate numeric integration solution is only one of several 
improvements brought on by a smaller transient timestep.) But 
decreasing the timestep may accentuate two undesirable characteris¬ 
tics of the simulator. Decreasing the timestep forces the simulator to 
solve more solution points and results in longer simulation runs, and 
decreasing the timestep tends to increase the chance of stepping into or 
close to a model discontinuity and failing to converge. In any simulator, 
a delicate balance exists between speed, accuracy, and the ability to com 
verge. This natural tension may be seen when attempting to select a 



TIME 

(a) 

Figure 4.9(a) Analytic solution vs. the numeric integration solution for the func¬ 
tion F{t) = t 2 - 5 1. (.Reprinted from Successfully Simulating Circuits With SPICE. 
Used with permission.) 
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nfl COMPARISON OF INTEGRATION METHODS FOR .5 SEC. TIMESTEP 
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Figure 4.9(0) Analytic solution vs. the numeric integration solution for the func¬ 
tion F(t) -t 2 - 5 t. (.Reprinted from Successfully Simulating Circuits With SPICE. 
Used with permission .) 



COMPARISON OF INTEGRATION METHODS FOR .25 SEC. TIMESTEP 
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Figure 4.9(c) Analytic solution vs. the numeric integration solution for the func¬ 
tion F(t ) = t 2 - 5 1. (.Reprinted from Successfully Simulating Circuits With SPICE. 
Used with permission.) 
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COMPARISON OF INTEGRATION METHODS FOR 1 SEC. TIMESTEP 
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Figure 4.10(a) Analytic solution vs. the numeric integration solution for the func¬ 
tion Fit) = 5 - 5e~*. (Reprinted from Successfully Simulating Circuits With SPICE. 
Used with permission.) 
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COMPARISON OF INTEGRATION METHODS FOR .5 SEC. TIMESTEP 
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Figure 4.10(h) Analytic solution vs. the numeric integration solution for the func¬ 
tion Fit) = 5 - 5e^. C Reprinted from Successfully Simulating Circuits With SPICE. 
Used with permission.) 
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COMPARISON OF INTEGRATION METHODS FOR .55 SEC. TIMESTEP 
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Figure 4.10(c) Analytic solution vs. the numeric integration solution for the func¬ 
tion Fit) = 5 - be*. (Reprinted from Successfully Simulating Circuits With SPICE. 
Used with permission.) 



TIME 

(a) 

Figure 4.11(a) Analytic solution vs. the numeric integration solution for the func¬ 
tion Fit) = sin(£). (Reprinted from Successfully Simulating Circuits With SPICE. 
Used with permission.) 
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TIME 

(b) 

Figure 4.11(d) Analytic solution vs. the numeric integration solution for the func¬ 
tion Fit) = sinU). (. Reprinted from Successfully Simulating Circuits With SPICE. 
Used with permission.) 
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(c) 

Figure 4.11(c) Analytic solution vs. the numeric integration solution for the func¬ 
tion Fit) = sinit). iReprinted from Successfully Simulating Circuits With SPICE. 
Used with permission .) 
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timestep for transient analysis which maintains accuracy, reasonable 
convergence characteristics, and reasonable simulation time. 

Numeric Integration of Electronic Circuits 

If we narrow our focus to the general class of functions most often seen 
in electronic circuits (linear, piecewise linear, exponential, and sinu¬ 
soidal), some general observations may be made about numeric inte¬ 
gration. 

First, the forward-Euler method is generally inaccurate and unsta¬ 
ble on these types of functions, and therefore is not used in SPICE. 
Second, remaining methods—backward-Euler, trapezoidal, and Gear— 
may be represented on a continuum with accuracy on the left-most side 
and stability on the right-most side as shown in Fig. 4.12a. The trape¬ 
zoidal is generally the most accurate; the Gear methods, because of 
their smoothing or averaging function, tend to be the most stable. (The 
backward-Euler method is sometimes referred to as a Gear-1 integra¬ 
tion method although no formal Gear-1 method exists.) But too much 
stability is not always good. On highly nonlinear circuits, past solution 
values may yield extremely poor estimates of future solution values. In 
general, the Gear-3, Gear-4, Gear-5, and Gear-6 integration methods 
show little or no improvement over the Gear-2 integration method, 
and, since these methods are more complex and involve more compu¬ 
tational steps, resulting in longer simulation runs, the Gear-3, Gear-4, 
Gear-5, and Gear-6 methods are not recommended for use in simula¬ 
tion. For general electronic circuit simulation, backward-Euler, trape¬ 
zoidal, and Gear-2 integration methods will produce the best results. 
These three methods are shown in Fig. 4.126. 

Before we go further, it is important to take note of the difference 
between Eqs. 4.10a and 4.106. Although these two equations are 
essentially the same, a is the equation of an integrator, while 6 is the 
equation of a differentiator. Integration tends to be a smoothing or 
averaging function, while differentiation tends to be a noisy and some¬ 
times erratic function. Although known as numeric integration, SPICE 
reverses the equality in the integration formula and uses numeric dif- 
ferentiation to determine the capacitor currents and inductor voltages. 

The step response is one measure of the stability of a system. By 
applying a step function to the numeric integration routines, several 
very important observations may be made. In Fig. 4.13, a ramp volt¬ 
age is applied to a 1-farad capacitor. The voltage increases at a con¬ 
stant 1 volt per second. At time t = 1, the voltage is held constant. 
Using analytic integration, the capacitor current can be calculated 
directly from Eq. 4.1. Both the capacitor voltage and current are illus¬ 
trated in Fig. 4.13. 
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Figure 4.12(a) Continuous exchange between accuracy and stability. 



Trapezoidal B.E. Gear-2 

(b) 


Figure 4.12(b) Suggested integration methods for SPICE sim¬ 
ulation. 


By replacing the analytic integration with a numerical integration, 
we can compare how each integration method performs when used as 
a numeric differentiator. 

Figures 4.14a and b illustrate the capacitor current response gener¬ 
ated with the application of the backward-Euler integration method. 
In Fig. 4.14a the backward-Euler integration method was applied with 
a fixed step size of .025S. Figure 4.146 is the same calculation with a 
fixed step size of ,0025s. The smaller step size improves the integration 
result but takes longer to calculate. 
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Figure 4.13 Capacitor voltage and calculated current response. C Reprinted from 
Successfully Simulating Circuits With SPICE. Used with permission.) 
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Figure 4.14(a) Capacitor voltage and current response calculated from backward- 
Euler integration and a step size of .025S. (Reprinted from Successfully Simulating 
Circuits With SPICE. Used with permission.) 



TIME 

(b) 

Figure 4.14(h) Capacitor voltage and current response calculated from backward- 
Euler integration and a step size of .0025S. (Reprinted from Successfully Simulat¬ 
ing Circuits With SPICE. Used with permission.) 
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Figures 4.15a and 6 illustrate the capacitor current response gener¬ 
ated with the use of the trapezoidal integration method. In Fig. 4.15a the 
capacitor current after time t- 1 appears to oscillate around the correct 
value of zero. While this result may surprise many readers, the effect of 
ringing is a well-known trapezoidal failure mechanism. In many func¬ 
tions which contain an abrupt change or discontinuous derivative, 
dV/dt, using the trapezoidal integration method with a large timestep 
results in the phenomenon known as trapezoidal oscillation. The oscilla¬ 
tion occurs because of the eigenvalues 5 of the circuit and the size of the 
timestep being used. The oscillation does not occur often but will occur 
when the eigenvalues of the circuit are small and the timestep is large. 
Trapezoidal oscillation is not a result of components in the circuit; the 
oscillation is simply an erroneous result of the numeric integration algo¬ 
rithm. Trapezoidal oscillation is a readily observable clue which indi¬ 
cates the numeric integration method is failing to predict the correct 
response, and corrective actions should be taken. 

Of the three integration methods available in SPICE, the trape¬ 
zoidal method is the only one which exhibits this oscillatory behavior. 
For most functions, the oscillation can be suppressed or eliminated by 
reducing the step size of the integration. By reducing the step size to 
.0025S, the oscillation seen in Fig. 4.156 is reduced significantly from 
the previous illustration. Reducing the step size further will eliminate 
the oscillation. 



TIME 

(a) 

figure 4.15(a) Capacitor voltage and current response calculated from trapezoidal 
integration and a step size of .025S. (Reprinted from Successfully Simulating Cir¬ 
cuits With SPICE. Used with permission .) 
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, „„„ CAPACITOR CURRENT USING TRAPEZOIDAL INTEGRATION 
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Figure 4.15(6) Capacitor voltage and current response calculated from trapezoidal 
integration and a step size of .0025S. C Reprinted from Successfully Simulating Cir¬ 
cuits With SPICE. Used with permission.) 



SPICE employs a dynamic timestep control algorithm and source- 
breakpoint adjustment (these will be discussed in detail in Chap. 5) to 
help minimize trapezoidal oscillation. For many simulations these 
techniques work well at controlling the oscillation, but, for some cir¬ 
cuits, trapezoidal oscillation will be readily observable. 

Figures 4.16a and b illustrate the capacitor current response gener¬ 
ated with the application of the Gear-2 integration method. In Fig. 
4.16a the capacitor current after time t - 1 appears to overshoot the 
correct value of zero. This overshoot is a well-known characteristic of 
the Gear integration method. Because the Gear integration method 
uses past timepoints to help determine the next solution point, circuit 
functions with an abrupt change often display a result which appears 
to overshoot the correct response. As in the case of trapezoidal oscilla¬ 
tion, the Gear-induced overshoot is not the result of components in the 
circuit; the overshoot is simply an error introduced by the integration 
algorithm. The Gear-2 integration method is the most stable integra¬ 
tion method in SPICE. However, the same mechanism which adds to 
the stability of the Gear method tends to resist abrupt changes. This 
resistance to change produces the overshoot seen in Fig. 4.16a. 

One way to reduce or eliminate the overshoot is to again reduce the 
size of the timestep being used during the integration. Figure 4.166 
illustrates how a much smaller step size reduces the overshoot at time 
t = 1. 
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CAPACITOR CURRENT USING GEARS INTEGRATION 
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Figure 4.16(a) Capacitor voltage and current response calculated from Gear-2 inte¬ 
gration and a step size of .025S. (.Reprinted from Successfully Simulating Circuits 
With SPICE, Used with permission .) 
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figure 4.16(6) Capacitor voltage and current response calculated from Gear-2 inte¬ 
gration and a step size of .0025S. (.Reprinted from Successfully Simulating Circuits 
With SPICE. Used with permission.) 
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Trapezoidal oscillation and Gear overshoot seen in the previous illus¬ 
trations are two basic failure mechanisms which may contribute inac¬ 
curacies to your transient simulation. These same characteristics are 
the kind of clues which users need to watch for during transient simu¬ 
lation. SPICE will not flag these as anomalous behaviors. As a SPICE 
user, you must observe these behaviors and make the appropriate 
adjustments to your simulation. 

Selecting an integration method 

Because SPICE uses only one integration method for the entire circuit, 
simulation users need to choose the best method for their circuit. The 
authors of SPICE originally chose the trapezoidal integration as the 
default integration method. The trapezoidal method is a good choice as 
the default method because, on the type of functions most often seen in 
electronic circuits, the trapezoidal tends to produce the most accurate 
result. In addition to the algorithm’s accuracy, the trapezoidal method 
is also faster to evaluate than the Gear integration methods. For many 
types of circuits, the trapezoidal integration method is the best choice 
for simulation. 

But the trapezoidal integration method is not always the best choice. 
In addition to the trapezoidal integration method, the authors of 
SPICE also implemented the backward-Euler and the Gear integra¬ 
tion methods and gave simulation users the option of selecting which 
integration method to use. 

Changing integration methods 

SPICE users may select the trapezoidal, backward-Euler, or Gear meth¬ 
ods of integration. The switch to change integration methods is appro¬ 
priately found in the .OPTIONS statement. To select either Gear-2 or 
trapezoidal integration methods, use the following option: 

.OPTIONS METHOD=TRAP (default) 

.OPTIONS METHOD=GEAR 

If the METHOD=GEAR option is chosen, SPICE will use the Gear-2 
integration method. Although this text does not recommend using the 
higher-order Gear methods (orders higher than 2), the higher-order 
Gear methods may be selected with the MAXORD option as shown 
below: 

.OPTIONS METHOD=GEAR MAX0RD=3 
.OPTIONS METHOD=GEAR MAX0RD=4 
.OPTIONS METHOD=GEAR MAX0RD=5 
.OPTIONS METHOD=GEAR MAX0RD=6 
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To select the backward-Euler integration method, use the MU 
option: 

.OPTIONS MU=0 

The backward-Euler integration method was added to SPICE in a 
later release, after Nagel 3 had finished his work on the program. For 
this reason, the MU option was never documented in the original 
SPICE2 User’s Guide. 4 For this reason, many simulation users are not 
aware of the MU option. 

In SPICE, the original trapezoidal algorithm has been modified to 
include the scaler MU in the formulation. Because of similarities in the 
formulation of the two algorithms, if the value of MU is 0.5, the formula 
reduces to the simple trapezoidal integration method. If MU is set to 0, 
the formula reduces to the backward-Euler integration method. And if 
MU is set between 0 and 0.5, a weight-averaged value of the two meth¬ 
ods is used. In SPICE, the value of MU is set by default to 0.5. Setting 
the option MU = 0 changes to backward-Euler integration, and setting 
MU between 0 and 0.5 combines the two integration methods. 


Detecting and Correcting Integration Failures 

None of the three integration methods in SPICE will always produce 
an accurate result on all transient circuits. All three integration meth¬ 
ods contain different failure mechanisms. Occasionally, during tran¬ 
sient analysis, certain voltage and current waveforms will excite one or 
more of the integration failure mechanisms. Exciting one of these fail¬ 
ure mechanisms will introduce error in the simulation output. The 
amount of the error may range from millivolts and picoamps to tens of 
volts or amps. It is up to the simulation user to detect numeric inte¬ 
gration errors and take corrective actions. SPICE does not have an 
internal mechanism to detect numeric integration failures and there¬ 
fore does not flag the erroneous output. Numeric integration failures 
must be detected in the simulation output by the user Only then can 
corrective measures be applied. 


Trapezoidal integration failures 

The trapezoidal integration method is the default method used in 
SPICE. The trapezoidal method combines good accuracy with low com¬ 
putational requirements. For many circuits, the trapezoidal integra¬ 
tion method is the best choice for simulation. 

When is the trapezoidal method not a good choice? When the trape¬ 
zoidal method starts to contribute significant error to the simulation 
result. There are two well-known failure mechanisms of the trape- 
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zoidal integration method that simulation users should learn to recog¬ 
nize. One failure mechanism is trapezoidal oscillation; you have 
already seen an example of this in Fig. 4.15a. The other failure mech¬ 
anism is accumulated error. Once the existence of these failure mecha¬ 
nisms is known, both mechanisms can be identified from characteristic 
features in the simulation output. The characteristic features of these 
failure mechanisms are easy to identify if you know what to look for in 
the output. By learning the characteristic features of the failure mech¬ 
anisms, SPICE users can learn to detect the presence of trapezoidal 
integration failure and swiftly switch to a more appropriate integra¬ 
tion method. 

Trapezoidal oscillation. Trapezoidal oscillation is a well-known, well- 
documented characteristic of the trapezoidal integration method. 5 
Trapezoidal oscillation occurs when the integration step size is too 
large to follow the curvature of a given function. The result of this fail¬ 
ure mechanism is a predicted solution which appears to oscillate 
around the correct solution. The disk files ch4-17a.cir and ch4-17b.cir 
contain netlists for the transient simulation of an MOS capacitor. Fig¬ 
ure 4.17a shows the simulation result from ch4-17a.cir. The simulation 
result of Fig. 4.17a appears to oscillate around the correct result pre¬ 
dicted by Fig. 4.176. In simulating the circuit ch4-17a.cir, trapezoidal 
oscillation introduces large amounts of error in the output. Simulate 
this circuit with the command: 

SIM CH4-17A.CIR 

Your result should match Fig. 4.17a. 

To eliminate trapezoidal oscillation from the simulation output, the 
SPICE user may switch to either Gear or backward-Euler integration. 
Neither Gear nor backward-Euler integration suffers from the oscilla¬ 
tion anomaly. Figure 4.176 shows the simulation result from the 
circuit ch4-17b.cir. The ch4-17b.cir circuit netlist is identical to 
ch4-17a.cir except for the addition of the option statement which calls 
the backward-Euler integration method (.OPTIONS MU=0). The 
backward-Euler method completely eliminates the oscillation. 

Simulate this circuit with the command: 

SIM CH4-17B.CIR 

Your result should match Fig. 4.176. 

The ch4-17c.cir circuit netlist is identical to ch4-17a.cir except for 
the addition of the option statement which calls the Gear-2 integration 
method (.OPTIONS METHOD=GEAR). The simulation results of ch4- 
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Figure 4.17(a) Transient simulation of an MOS capacitor. {Reprinted from Successfully 
Simulating Circuits With SPICE. Used with permission.) 
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Figure 4.17(0) Corrected transient simulation of an MOS capacitor. {Reprinted from 
Successfully Simulating Circuits With SPICE. Used with permission.) 
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17c.cir are shown in Fig. 4.17c. Figure 4.17c shows the simulation 
results of ch4-17c.cir and illustrates a distinct overshoot at both of the 
transition points. During the simulation of ch4-17c.cir, and on many 
switching circuits, the Gear integration methods, because of their 
averaging effect, tend to produce overshoot where circuit voltages and 
currents change abruptly. Like the oscillation caused by trapezoidal 
integration, the overshoot is not present in the actual circuit. The over¬ 
shoot is merely an error introduced by the Gear integration method. 

While either the Gear or backward-Euler integration methods 
could be chosen to eliminate trapezoidal oscillation, for this circuit 
the backward-Euler method is the best choice. In comparing the sim¬ 
ulation results from all three simulations, and by knowing what the 
simulated output should be, the backward-Euler integration method 
is the best choice for this simulation. 

Several SPICE-like simulators do not allow users to switch to Gear or 
backward-Euler integration. What can the users of these tools do to cor¬ 
rect trapezoidal oscillation problems? The answer is to reduce the max¬ 
imum transient timestep. (Note: Chap. 5 will present a detailed look at 
how to limit the maximum transient timestep.) Figures 4.146, 4.156, 
and 4.166 demonstrate how reducing the transient step size eliminates, 
or reduces, both trapezoidal oscillation and Gear overshoot. Decreasing 
the maximum transient step size is another way of reducing the error 
and improving the accuracy of all the integration methods in SPICE. 



TIME 

(c) 

Figure 4 . 17 (c) Transient simulation using Gear-2 integration. (Reprinted from Success¬ 
fully Simulating Circuits With SPICE. Used with permission.) 
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But decreasing the timestep is not always desirable. In Chap. 3, it 
was shown that decreasing the timestep made the simulator run 
longer and made simulations more prone to nonconvergence problems. 
Whenever possible, changing integration methods is the preferred way 
to correct numeric integration problems. 

Trapezoidal oscillation is caused by a failure of the numeric integra¬ 
tion method, not by an improper circuit description or device model. At 
times, as in Fig. 4.17a, trapezoidal oscillation may be large and radi¬ 
cally unrealistic; at other times, the oscillation may be very subtle. Not 
all trapezoidal oscillation is as pronounced as the previous example. 
Figure 4.18a shows the simulation result of the circuit file ch4-18a.cir. 
The subtle ripple in the simulation output is the result of trapezoidal 
oscillation. Figure 4.186 shows the result of the same simulation with 
the Gear integration method. 

Simulate these circuits with the command: 

SIM CH4-18A.CIR 

then 

SIM CH4-18B.CIR 

Your results should match Figs. 4.18a and b. 
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Figure 4.18(a) Transient simulation with trapezoidal oscillation. C Reprinted from Suc¬ 
cessfully Simulating Circuits With SPICE. Used with permission.) 
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Figure 4.18(d) Corrected transient simulation. (Reprinted from Successfully Simulating 
Circuits With SPICE. Used with permission.) 


Because SPICE does not have an algorithm to detect trapezoidal 
oscillation, the program will not warn users of the integration method 
failure. Trapezoidal oscillation is a phenomenon which must be 
observed by the user. The clues to look for are unexplained or unex¬ 
pected oscillation, especially oscillation around what should be the cor¬ 
rect result. This naturally leads to the question of how you know 
whether the oscillation seen in a simulation output is real or merely 
failure of the trapezoidal algorithm. The best way to differentiate 
between the two is to rerun the simulation with either the Gear or 
backward-Euler integration method. With the new integration method 
selected, the oscillation will disappear if it was caused by a failure of 
the trapezoidal algorithm. If the circuit response is the cause of the 
oscillation, all three integration methods—trapezoidal, Gear, and 
backward-Euler—will reproduce the oscillitory behavior. 

Accumulated error. Accumulated error is a second failure mechanism 
of the trapezoidal integration method. Accumulated error usually 
occurs in periodic circuits and long transient simulations. The result of 
this failure mechanism is usually an unexplainable circuit behavior. 

A typical transient simulation may require thousands of timepoint 
calculations. At each timepoint, the numeric integration algorithm 
may introduce a finite amount of error in the solution. At some time- 
points, the error will overestimate the exact solution, while at other 
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timepoints the error will underestimate the exact solution. Since the 
error from the previous solution point is carried forward to the next 
solution point, the overestimation/underestimation errors usually can¬ 
cel each other over a large number of timepoints. But sometimes the 
accumulated error tends to increase with each new timepoint. When 
this occurs, the simulation result deviates from the correct solution. 
Often the simulation output appears to be completely inaccurate. In 
many cases, the erroneous output will appear to defy the laws of 
physics (or electronics). 6 

The circuit file ch4-19a.cir is a listing for the transient behavior of a 
555 timer. In the ch4-19a.cir netlist, an external capacitor and several 
external resistors are connected to the timer. In this configuration, the 
timer behaves like an asynchronous oscillator. The simulation results 
are shown in Fig. 4.19a. In this simulation, accumulated error intro¬ 
duces a mathematically stable state which is not present in the real 
555 timer. Because the circuit is an asynchronous oscillator, once the 
integration error accumulates to a point where the circuit is stable, the 
output remains in that state for the rest of the simulation. Simulate 
this circuit with the command: 

SIM CH4-19A.CIR 

Your result should match Fig. 4.19a. 

The circuit file ch4-19b.cir is the same 555 timer circuit, with the 
addition of the option statement which calls the Gear-2 integration 
method. The result of this simulation is shown in Fig. 4.195. Simulate 
this circuit with the command: 

SIM CH4-19B.CIR 

Your result should match Fig. 4.195. 

Accumulated error in trapezoidal integration is usually character¬ 
ized by an unexplainable simulation behavior. In these cases, switch¬ 
ing to a more stable integration method (Gear-2) or reducing the 
maximum transient timestep minimizes the accumulated error and 
restores the normal circuit behavior. (Author’s note: Experience has 
shown that accumulated numeric integration error is not the only 
cause of unexplainable circuit behavior!) 


Backward-Euler integration failures 

Although the backward-Euler method does not suffer from the numer¬ 
ical oscillation of the trapezoidal method, it may suffer from larger 
amounts of local truncation error when applied to nonlinear wave- 
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Figure 4.19(a) Accumulated error simulation failure. (.Reprinted from Successfully Sim 
ulating Circuits With SPICE. Used with permission.) 
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Figure 4.19(b) Corrected transient simulation. (Reprinted from Successfully Simulating 
Circuits With SPICE. Used with permission .) 
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forms. For sinusoidal and exponential waveforms, the backward-Euler 
integration method may generate large amounts of error which lead to 
accumulated-error failures and unexpected simulation results. For 
these types of circuits, the backward-Euler integration method is not 
suggested. The backward-Euler is best used on linear and piecewise 
linear waveforms. The backward-Euler method is a good choice for dig¬ 
ital circuits. 

The circuit file ch4-20a.cir is a listing for the transient behavior of an 
ideal resonator circuit. The circuit is composed of a single capacitor 
and inductor. The transient simulation is started with an initial cur¬ 
rent flowing through the inductor. Because both of the elements are 
ideal (no parasitic leakages), the circuit should continue to oscillate for 
the entire transient simulation. Figure 4.20a shows the simulation 
result when using the backward-Euler integration method. For sinu¬ 
soidal and nonlinear waveforms, the backward-Euler method may 
introduce large amounts of local truncation error in the simulation. On 
these circuits, the local truncation error quickly accumulates and 
results in an accumulated-error failure. In the simulation of ch4- 
20a.cir, the error quickly dampens the ideal oscillation. Simulate this 
circuit with the command: 

SIM CH4-20A.CIR 

Your result should match Fig. 4.20a. 

The circuit file ch4-20b.cir is the same oscillator, except for the 
option which resets the simulator to the trapezoidal integration 
method. Simulate this circuit with the command: 

SIM CH4-20B.CIR 

Your result should match Fig. 4.206. Figure 4.206 illustrates proper 
simulation output. 

Like trapezoidal accumulated-error failures, accumulated error in 
backward-Euler integration is characterized by an unexpected or 
unexplained simulation result. Backward-Euler numeric integration is 
especially prone to accumulated-error failures when simulating non¬ 
linear or sinusoidal waveforms, and is best used on linear or piecewise 
linear waveforms. 

Gear integration failures 

Like backward-Euler integration, Gear integration does not suffer 
from trapezoidal oscillation, but Gear integration may suffer from 
large amounts of local truncation error, especially on circuits which 
contain highly nonlinear or switching waveforms. 



Figure 4.20(a) Accumulated error simulation failure. {Reprinted from Successfully Sim 
ulating Circuits With SPICE. Used with permission.) 
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Figure 4.20(d) Corrected transient simulation. (.Reprinted from Successfully Simulatinj 
Circuits With SPICE. Used with permission.) 
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The Gear integration method uses an averaging technique to predict 
the next solution point. In theory, this averaging technique should 
make the algorithm more stable (less accumulated error over a large 
number of timepoints). But, because the Gear method uses past time- 
point values to help predict where the next solution point should be 
found, for switching waveforms and piecewise linear waveforms, the 
Gear method may overshoot the correct solution value. (You have 
already seen an example of this in Figs. 4.16a and 4.17c.) Pronounced 
overshoot in switching circuits is a characteristic failure of Gear 
method integration in SPICE. 

Many switching circuits exhibit overshoot at the transition points. 
Once again, simulation users face the dilemma of whether the over¬ 
shoot seen in their simulation is real or simply a Gear integration 
error. The best way to determine the true cause of the overshoot is to 
switch integration methods. Neither trapezoidal nor backward-Euler 
integration methods generate the overshoot of the Gear integration 
method. Rerun the simulation with one of the other integration meth¬ 
ods. If the overshoot was caused by the mechanics of the circuit, the 
simulation result should be identical to the Gear integration results. If 
the overshoot was caused by a Gear-introduced error, the overshoot 
will disappear. 

Suggested Integration Settings 

When simulating transient circuits, users must learn to look for 
numeric integration anomalies. Table 4.2 is a list of suggested integra¬ 
tion methods for different types of circuits. Table 4.3 is a summary of 
numeric integration failure mechanisms and corresponding corrective 
actions. In Table 4.3, the list of corrective actions includes reducing the 
maximum transient timestep. Setting the maximum timestep is the 
subject of the next chapter. Read Chap. 5 before using this technique to 
correct numeric integration problems. Use these tables to minimize 
numeric integration failures and correct integration failures if they do 
occur. 


TABLE 4.2 Suggested integration Methods 


Circuit waveforms 

Suggested integration method 

Sinusoidal (sine wave circuits, oscillators, general 
amplifiers, power supplies, transmission lines) 

Trapezoidal 

Linear and piecewise linear (pulse circuits, digital 
and logic) 

Trapezoidal or backward-Euler 

Exponential (bipolar and diode amps, RC tree 
networks) 

Trapezoidal or Gear 
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TABLE 4.3 Common Numeric Integration Failures 


Failure mechanism 

Corrective action 

Trapezoidal oscillation 

Switch to backward-Euler or Gear method or 
reduce maximum timestep 

Trapezoidal accumulated error 

Switch to Gear method or reduce maximum 
timestep 

Gear overshoot 

Switch to backward-Euler or trapezoidal 
integration or reduce maximum timestep 

Gear Accumulated error 

Reduce maximum timestep 

Backward-Euler accumulated error 

Switch to trapezoidal or Gear method or 
reduce maximum timestep 


Numeric Integration in Other Simulators 

Hspice 

Meta-Software’s Hspice contains all three integration methods of stan¬ 
dard SPICE, but when Hspice users simulate ch4-17a.cir and ch4- 
18a.cir they will not see the anomalous behavior seen in the standard 
SPICE output. The explanation for this has to do more with how 
Hspice performs transient analysis than how the integration methods 
are applied. 7 

Standard SPICE and all the other simulators reviewed here, by 
default, use a timestep control algorithm known as the Local Trunca¬ 
tion Timestep Control algorithm. Meta-Software has developed a pro¬ 
prietary timestep control algorithm for transient analysis known as 
DVDT. The DVDT algorithm improves on several of the deficiencies of 
the standard timestep control algorithm and as a result tends to suffer 
from fewer numeric integration problems. (Remember the relationship 
between step size and numeric integration accuracy.) 

A second reason Hspice does not suffer from the same problems is 
that many of the I-V and C-V device equations have been rewritten and 
many of the sharp edges of the equations (discontinuities) have been 
smoothed. Since these same sharp edges initiate trapezoidal oscilla¬ 
tion, a reduction in discontinuities reduces the occurrence of trape¬ 
zoidal oscillation. (Hspice users will notice how the Hspice simulation 
of circuit ch4-17a.cir produces a very different looking output.) 

IS_Spice 

Intusoft’s IS_Spice is based directly on SPICE 2G.6. IS_Spice offers 
users the standard choices of trapezoidal, backward-Euler, and Gear 
integration methods. Trapezoidal integration is the default method in 
IS_Spice. 
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Micro-Cap IV 

Spectrum Software’s Micro-Cap IV is also based on SPICE2G.6 but has 
been completely rewritten in C. Micro-Cap IV does not support either 
backward-Euler or the Gear integration methods. 9 If (trapezoidal) 
integration problems do occur, users must reduce the maximum time- 
step of the program. Micro-Cap IV users should simulate ch4-17a.cir 
and ch4-18a.cir to demonstrate the effects of numeric integration fail¬ 
ure. Then resimulate the circuits while reducing the maximum time- 
step until the errors disappear. 

Pspice 

MicroSim’s Pspice uses either trapezoidal or Gear-2 numerical integra¬ 
tion. But unlike standard SPICE, Pspice does not offer users an option 
to manually select the integration method. By default, Pspice starts a 
transient simulation with trapezoidal numeric integration. As the sim¬ 
ulation progresses, the integration routines in Pspice look for the 
numerical conditions which might induce trapezoidal oscillation. If the 
conditions for trapezoidal oscillation occur, Pspice automatically changes 
to Gear-2 integration before the oscillation begins. 10 

Summary 

Any of the integration methods in SPICE will produce a finite amount 
of error at each timepoint in the transient analysis. In many simula¬ 
tions, the error is small. In these cases, the overall accuracy of the sim¬ 
ulation result will be within the error tolerances of the program. But 
when the integration inaccuracies become large, the simulation result 
will deviate from the correct solution. The amount of deviation may be 
infinitesimally small or orders of magnitude larger than the original 
signal! When analyzing simulation results, SPICE users must learn to 
recognize numeric integration errors. SPICE does not have an algo¬ 
rithm to detect and warn users of numeric integration failure. Be cau¬ 
tious and question your output. Does the simulation generate the 
result you expect? If not, scrutinize the simulation. Question unex¬ 
pected behavior. Being aware of integration-caused failures and know¬ 
ing how to resolve these failures are the watermarks of a true 
simulation craftsperson. 

The previous sections describe the type of errors that SPICE 
numeric integration methods are most likely to produce. These are the 
types of erroneous behavior SPICE users must learn to look for. SPICE 
users should never run a simulation without some reasonable expecta¬ 
tion of what the output should be. This statement is important, 
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because without a reasonable expectation of the output, detecting 
numeric integration failures may be impossible. 

Because SPICE and almost all of the SPICE-like simulators use 
trapezoidal integration as the default integration method, you should 
begin looking for trapezoidal integration failures first. The two inte¬ 
gration failures associated with trapezoidal integration are trapezoidal 
oscillation and accumulated error. Of these two failure mechanisms, 
trapezoidal oscillation is the most common. Figures 4.17a and 4.18a 
are illustrations of trapezoidal oscillation. Users should also look for 
signs of accumulated integration error. Figures 4.19a and 4.20a are 
illustrations of accumulated-error failures. 

When you switch to the Gear integration method, watch for unex¬ 
pected overshoot in switching circuits, and watch for accumulated 
error. When you switch to the backward-Euler integration method, 
watch for accumulated error failures. 

In SPICE, the numeric integration algorithms routinely produce 
stunningly accurate, highly reliable results during transient simula¬ 
tion. But occasionally, the numeric integration routines will generate 
erroneous results. Since SPICE does not automatically detect numeric 
integration errors, the responsibility of detecting and correcting these 
errors falls on the user. Most numeric integration failures can be read¬ 
ily identified if users know to look for them, and all numeric integra¬ 
tion failures can be fixed. Correcting the integration failure may 
involve switching algorithms or reducing the maximum internal 
timestep. Becoming a knowledgeable SPICE user involves much more 
than simply learning to create input netlists. It also involves knowing 
the strengths and weaknesses of SPICE, and knowing how to compen¬ 
sate for the weaker facets of the program. 
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Chapter 

5 

Timestep Control 


The timestep control algorithms determine the timepoints where 
SPICE solves the circuit equations. The selection of timepoints relates 
directly to the accuracy of the numeric integration routines and indi¬ 
rectly to the convergence properties of the simulator. For these rea¬ 
sons, the timestep control algorithms of SPICE are the most important 
routines of the transient simulation. But like so many other routines in 
SPICE, the timestep control algorithms are not foolproof. Accurate, 
high-quality transient simulations require a precise tailoring of the 
timestep control routines. 

This chapter will present a look at the two timestep control algo¬ 
rithms in SPICE. Like numeric integration, each algorithm has advan¬ 
tages and disadvantages, each algorithm will work well with certain 
types of circuits, and each algorithm will suffer from a unique set of 
problems. Knowing which timestep control algorithm to use on a spe¬ 
cific circuit and knowing the limitations of each algorithm is the next 
lesson in becoming a true simulation craftsperson. 

Timestep Control in the Early Years 

In the late 1960s and early 1970s when CANCER and SPICE 1 were 
being used to simulate circuits, the transient command: 

.TRAN INS 100NS 

forced the simulator to solve the circuit equations at each InS time 
interval. This is known as constant timestep control. With this com¬ 
mand, both CANCER and SPICE 1 solved the circuit equations and 
printed the output results every InS during the transient simulation. 
Figure 5.1a illustrates a constant timestep control simulation. For 
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these early simulators, the print interval and the internal timepoints 
were identical. 

But in sweeping the timestep at a constant rate, two problems were 
soon observed. First, circuits with large, rapid voltage and current 
transitions often failed to converge at the switch point. (This is due to 
the same nonconvergence mechanism which causes the DC sweep 
analysis to fail around voltage transitions. See Chap. 3 for a thorough 
discussion of why this occurs.) Only by significantly reducing the step 
size would the simulation successfully cross the voltage transition. But 
small timesteps generate extremely long simulation times and may 
aggravate nonconvergence problems. The second problem with con¬ 
stant timestep control relates to the numeric integration algorithms. 
With a constant step size, numeric integration algorithms produce the 
most accurate result when the circuit is in a stable state, and the inte¬ 
gration routines produce the least accurate solution when the circuit is 
in a transition state. With a constant step size, the error introduced by 
the numeric integration algorithms is proportional to the transition 
rate of circuit voltages and currents. For circuit simulation, constant 
timestep control forces users to select between a slow (long run time), 
accurate simulation and a fast, inaccurate one. 

Improving the Timestep Control Algorithm 

The introduction of SPICE2 brought two different forms of dynamic 
timestep control to circuit simulation. Both of the timestep control 
algorithms decrease the step size when the circuit approaches a volt- 



Figure 5.1(a) Simulation with constant timestep control. 
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age or current transition and then increase the step size when the 
circuit begins to stabilize again. These dynamic timestep control algo¬ 
rithms are substantially better than simple constant timestep control 
and significantly increase both the accuracy and speed of a transient 
simulation. 

But, like the numeric integration algorithms, neither of the timestep 
control algorithms are completely foolproof. Both timestep control rou¬ 
tines have problems on certain types of circuits, and these problems 
may introduce significant error in the simulation output. Even with 
completely accurate device models, failure mechanisms within the 
timestep control algorithms may introduce significant error into the 
simulation result. A knowledgeable simulation user is aware of these 
possible problems and, when questionable simulation results arise, 
knows the corrective action to employ. 

An Overview of Dynamic Timestep Control 

The two timestep control algorithms in SPICE both have basically the 
same purpose. The timestep control algorithms must alter the step size 
so, during rapid voltage transitions, the timestep is decreased to avoid 
nonconvergence and maintain accuracy through the transition region, 
and, during periods of little or no circuit activity, the timestep is 
increased to speed the simulation to a conclusion. But there are subtle 
differences between the two timestep control algorithms which may 
amount to significantly different simulation results. 

The dynamic timestep control algorithms significantly changed the 
way SPICE computes transient timepoints, but the transient com¬ 
mand did not change when the timestep control routines changed. This 
has lead to a slightly confusing situation. The transient command 

.TRAN INS 100NS 

still seems to imply that SPICE solves the transient simulation at 
evenly spaced InS intervals. This command is the same command the 
earlier version of CANCER and SPICE 1 used with their constant 
timestep control routines. The dynamic timestep control routines do 
not solve the circuit equations at evenly spaced intervals. The timestep 
control routines perform the task of deciding where the circuit equa¬ 
tions should be solved. As the transient analysis begins, SPICE uses 
Eq. 5.1 to compute the first transient timepoint. After the program 
computes the first timepoint, the timestep control algorithms assume 
the responsibility for increasing or decreasing the size of the timestep. 

STEP(l) - transient duration ^ ^ 

50 
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As SPICE continues the transient simulation, the timestep control 
algorithms monitor and adjust the size of the timestep. At each com¬ 
puted timepoint, SPICE saves the solution values in memory. At the 
and of the transient simulation, SPICE uses a simple linear interpola¬ 
tion routine and the data saved from the transient simulation to gen¬ 
erate the evenly spaced InS output result. Figure 5.16 illustrates a 
dynamic timestep control simulation. In most simulations, the number 
3 f timepoints SPICE solves during the transient analysis is much 
higher than the number of output points. The InS print interval has 
very little influence on the selection of timepoint values. 

Both timestep control algorithms found in SPICE use three sensors 
to control the size of the timestep. The sensors monitor the rate of 
change of circuit voltages and currents, nonconvergence timepoints, 
and source breakpoints. Of these, the first, the rate-of-change sensor, is 
the predominant mechanism which influences the selection of time- 
steps. The others modify the timestep predicted by the rate-of-change 
sensor. In addition to the control sensors, the timestep control algo¬ 
rithms also have upper and lower limits on the size of the timestep. All 
these mechanisms work together to predict the best timestep for the 
simulation. 
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Figure 5.1(b) Simulation with dynamic timestep control. 

Dynamic timestep control increases both simulation speed, by sampling fewer 
timepoints, and simulation accuracy, by sampling more points where they are 
needed. SPICE offers users two types of dynamic timestep control. 

When using dynamic timestep control, the values printed in the output file are 
interpolated from the two closest solution points. The print interval has little influ¬ 
ence on the selection of solution timepoints. 
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Timestep control and circuit activity 

Both timestep control algorithms have a method of measuring the rate 
of change of the circuit voltages and currents, and this measurement 
will be used to increase or decrease the current timestep. The timestep 
is decreased during voltage and current transitions and increased dur¬ 
ing periods of lesser circuit activity. Each of the timestep control algo¬ 
rithms uses different measurements to determine the circuit’s rate of 
change. One of the algorithms monitors the number of iterations at 
each new solution point. During rapid circuit transitions, the change in 
voltage and current levels requires more Newton-Raphson iterations 
to achieve convergence. As the number of iterations increases, the 
timestep is decreased. During periods of latency, far fewer iterations 
are required because the circuit voltages are not changing. Under 
these conditions, the timestep can be increased. As the analysis pro¬ 
gresses, the step size is adjusted, based on the number of iterations 
used at the past timepoint. The other algorithm uses a routine which 
monitors the amount of error the numeric integration routines are pro¬ 
ducing. During periods of little or no circuit activity, the numeric inte¬ 
gration routines are extremely accurate and produce very little error. 
Because the error is low, the timestep control algorithm increases 
the step size. As the circuit activity increases, the error increases, and 
this forces the timestep control algorithm to decrease the timestep. 
Throughout the entire transient analysis, the rate-of-change sensors 
monitor the circuit activity and make appropriate adjustments to the 
timestep. 

Timestep control 

and nonconvergent timepoints 

Both timestep control algorithms measure the number of iterations 
required to converge on a solution timepoint. If a given timepoint ever 
fails to converge to a solution within ITL4 iterations, both timestep con¬ 
trol algorithms will discard the failed timepoint, reduce the previous 
timestep to Vs of the original value and then reattempt the solution at a 
timepoint closer to the last convergent timepoint. Nagel et al. 1 deter¬ 
mined that when the transient simulation failed to converge within a 
set number of iterations, the failure could usually be attributed to a 
fast circuit transition. Because of this, Nagel reduced the timestep to % 
of the original value and reattempted a solution at the new timepoint. 
Figure 5.2 illustrates how SPICE recalculates a new timepoint after 
failing to converge. 

The ability to discard nonconvergent timepoints and reattempt a 
new solution is a significant improvement over both the DC operating 
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Figure 5.2 Attempting a new solution point after a nonconvergent timepoint. 


point and the DC sweep solution algorithms. Only during transient 
analysis will the simulator keep searching for a solution after a non¬ 
convergent timepoint has been found. SPICE will continue to search 
for a new solution point by reducing the timestep over and over until 
either the simulator converges or the size of the timestep is reduced 
below an allowed minimum. Only when the latter occurs does SPICE 
halt the solution algorithm and print the “internal timestep too small” 
warning message. 


Timestep control and source 
breakpoint adjustment 

Both pulse and PWL (piecewise linear) circuit sources generate pulse 
waveforms with sharp leading and trailing edges. Often, a change in 
the state of the source, such as a leading edge or trailing edge, is a pre¬ 
cursor to circuit transition. Because of this, SPICE generates a break¬ 
point table before the start of the transient simulation. The table 
contains the timepoints corresponding to the beginning and end of 
each rising or falling edge for pulse or PWL sources. During transient 
analysis, SPICE forces the timestep control algorithm to solve the cir¬ 
cuit at each breakpoint, and automatically reduces the timestep fol¬ 
lowing a breakpoint calculation. Forcing the simulator to solve the 
circuit at the source breakpoints greatly increases the accuracy of the 
transient output. Figure 5.3 illustrates the breakpoints of a pulse volt¬ 
age source. 
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Minimum and maximum step sizes 

Each of the timestep control algorithms has an upper and lower limit 
on the size of the timestep. During long periods of latency, the timestep 
control algorithm continues to try to increase the timestep. The 
timestep will be increased until it reaches the maximum allowed step 
size. Even with no circuit activity, the timestep will not exceed this 
maximum limit. If a solution point fails to converge, SPICE will reduce 
the size of the timestep by 8 and reattempt the solution. If the new 
timepoint fails to converge, the timestep will be reduced by 8 again. 
This process repeats until a convergence solution is found or until the 
step size reaches or becomes less than the minimum. If the step size is 
reduced below the allowed minimum, SPICE prints the “internal 
timestep too small” error message and aborts the simulation. 

For both timestep control algorithms, the mechanisms just described 
determine the size of the timestep and where SPICE tries to solve the 
circuit equations. The rate-of-change sensors select a timestep which is 
appropriate for the circuit. If the new timepoint fails to converge, the 
timestep is reduced by a factor of 8 and the solution reattempted. If the 
timestep is close to a source breakpoint, the step size is adjusted to 
align with the breakpoint. An inactive period of the simulation may 
increase the timestep to the maximum allowed limit, or a nonconver¬ 
gence problem may reduce the timestep below the minimum allowed 
limit and abort the simulation. These are the mechanisms which deter¬ 
mine the step size during transient analysis. 
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Iteration-Count Timestep Control 

The first of the two timestep control algorithms is the Iteration-Count 
method. The name comes from the way SPICE monitors the circuit 
activity. 

The Iteration-Count timestep control algorithm monitors the num¬ 
ber of iterations at each timepoint during the transient simulation. 
With the circuit in a stable condition, the number of iterations required 
to solve the circuit equations at each new timepoint is usually small. 
When the number of iterations becomes smaller than the lower thresh¬ 
old (ITL3), the timestep control algorithm doubles the size of the 
timestep before calculating the next timepoint. A low number of itera¬ 
tions indicates the circuit is relatively inactive and therefore a larger 
timestep can be used. 

When the circuit begins to enter a transition or the circuit activity 
starts to increase, the number of iterations at each timepoint increases. 
When the number of iterations exceeds the upper limit (ITL4), the 
timestep control algorithm discards the current timepoint, cuts the 
timestep to % of the present value, and reattempts the solution with 
the smaller timestep. 

If the circuit enters a latent phase and the number of iterations is 
lower than ITL3, the timestep control routine will try to increase the 
timestep. If the circuit remains in a latent stage, the timestep control 
will continue to increase the step size until the size of the timestep 
reaches its maximum limit. With the Iteration-Count method, the 
timestep is limited to become no larger than the print interval speci¬ 
fied on the transient statement. The print interval determines the 
maximum timestep for the Iteration-Count method. 

If the current timepoint fails to converge, the timestep will be 
reduced to Ys of the previous value and the solution reattempted at the 
new timepoint. If the new timepoint fails, the timestep will again be 
cut. This process will continue until either a convergent timepoint is 
found, or the step size is reduced below the allowed minimum. For the 
Iteration-Count timestep control algorithm, the minimum allowed step 
size is TSTOP/50e9, where TSTOP is the duration of the transient 
analysis. 

The Iteration-Count timestep control algorithm adjusts the timestep 
by using an upper and lower threshold. When the number of iterations 
is greater than the upper threshold, ITL4, the timestep is reduced by a 
factor of 8. When the number of iterations is less than the lower thresh¬ 
old, ITL3, the timestep is increased by a factor of 2. The Iteration-Count 
algorithm adjusts the timestep to keep the number of iterations at each 
timepoint between these two thresholds. 

Regardless of the circuit activity (or inactivity), the timestep is never 
allowed to become larger than the print interval. If the timestep is 



Timestep Control 135 


reduced below the allowed minimum, SPICE prints the “internal 
timestep too small” error message and aborts the simulation. 

Both ITL3 and ITL4 are SPICE options and may be adjusted by the 
user. ITL3 is set by default to 4 iterations, and ITL4 is set by default 
to 10 iterations. The Iteration-Count timestep control algorithm is 
selected by setting the option LVLTIM to 1. 


Local Truncatioh Error Timestep Control 

The second type of timestep control is the Local Truncation Error 
timestep control algorithm. The Local Truncation Error method is the 
default method used by SPICE and most of the SPICE-like simulators. 
The Local Truncation Error routine uses a formula which predicts the 
magnitude of the error computed in the numeric integration calcula¬ 
tions of the previous timepoint. The Local Truncation Error timestep 
control algorithm uses the amount of error being generated in the 
numeric integration routines to adjust the timestep. Equation 5.2 illus¬ 
trates the equation relating the calculated local truncation error and 
the timestep. The expression LTE in Eq. 5.2 is the computed local trun¬ 
cation error. As the truncation error increases, the timestep decreases; 
as the truncation error becomes smaller, the timestep grows. This algo¬ 
rithm maintains the timestep so the amount of error in the numeric 
integration routine is maintained within an allowed limit. 


STEP(T + 1) = / tRTOL*(RELTOL* I I e (T) 1 +ABSTOL) (5 2) 
V MAX(ABSTOL, LTE(T)) 

In Eq. 5.2, the computed local truncation error is scaled by the 
options ABSTOL, RELTOL, and TRTOL. ABSTOL and RELTOL are 
the convergence error tolerance options presented in Chap. 3; TRTOL 
will be discussed later in this chapter. Equation 5.2 re-illustrates the 
importance of setting the error tolerance parameters to appropriate 
levels for the voltages and currents in your circuit. Inappropriate set¬ 
tings for ABSTOL and RELTOL may force SPICE to use an unusually 
small timestep during transient simulation. While a small timestep 
will help ensure accuracy of the simulation, an excessively small 
timestep will cause long run times and may aggravate nonconvergence 
problems. 

In addition to the local truncation error rate-of-change sensor, the 
Local Truncation Error timestep control algorithm also uses ITL4 as a 
limit on the number of iterations allowed at a given timepoint. If the 
number of iterations at a timepoint is greater than ITL4, the timestep 
control algorithm cuts the timestep to % of the original value and reat¬ 
tempts the solution at the new timepoint. 
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As the circuit activity decreases, the local truncation error intro¬ 
duced by the numeric integration routines decreases. As the local trun¬ 
cation error decreases, the timestep control routine tries to increase 
the timestep. If the circuit remains stable, the timestep control routine 
continues to increase the step size until the step size reaches the max¬ 
imum allowed limit. For the Local Truncation Error timestep control 
routine, the maximum allowed timestep is TSTOP/50, where TSTOP is 
the transient duration. 

If the current timepoint fails to converge, the timestep will be 
reduced to % of the previous value and the solution reattempted at the 
new timepoint. If the new timepoint fails, the timestep will again be 
cut. This process will continue until either a convergent timepoint is 
found, or the step size is reduced below the allowed minimum. For the 
Local Truncation Error timestep control algorithm, the minimum 
allowed step size is TSTOP/50e9, where TSTOP is the total transient 
duration. Both of the timestep control algorithms in SPICE use the 
same minimum timestep limit. 

The Local Truncation Error timestep control algorithm is the default 
algorithm in SPICE. The Local Truncation Error method is selected 
with the option LVLTIM=2. 

A Comparison of Timestep Control Methods 

While both timestep control algorithms are good, the Local Truncation 
Error algorithm should produce a more accurate transient analysis 
because of the close coupling between the numeric integration trunca¬ 
tion error and the size of the timestep. The Iteration-Count method has 
no measure of the numeric integration error and therefore no means to 
adjust the timestep with this error. 

But the Local Truncation Error algorithm will only provide a supe¬ 
rior result if the estimate for local truncation error is accurate. The for¬ 
mula which SPICE uses to calculate the local truncation error is 
known as the Third Divided Difference formula. 1 Experimentally, 
Nagel et al. noticed that the local truncation error predicted from the 
Third Divided Difference equation was generally about a factor of 7 
over what the actual error should be. For this reason, Nagel added the 
option TRTOL to the Local Truncation Error timestep formula (Eq. 
5.2). TRTOL is set by default to 7. Even with this correction factor, for 
certain types of circuits the estimate for local truncation error is 
extremely poor. For sinusoidal circuits and inductive circuits, the 
Third Divided Difference formula and Eq. 5.2 lead to an extremely poor 
estimate of timestep. Because of this, both of these types of circuits 
may produce better results when simulated with the Iteration-Count 
timestep control algorithm. 
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Even when the estimate of local truncation error is accurate, a 
second failure mechanism exists in the Local Truncation Error time- 
step control algorithm. The Local Truncation Error timestep control 
algorithm limits the maximum timestep to TSTOP/50, whereas the 
Iteration-Count timestep control algorithm limits the maximum time- 
step to TSTEP. 

TSTEP is the print interval and TSTOP is the transient duration. 
For circuits with more than 50 solution points, the Local Truncation 
Error timestep control algorithm allows the timestep to become much 
larger than the print interval. For simulations requiring hundreds or 
thousands of solution points, the Local Truncation Error timestep con¬ 
trol algorithm may use a step size much larger than the print interval. 
When this occurs, the simulation result is undersampled or aliased 
and the output appears choppy and distorted. 

Timestep Control Failures 

Figure 5.4 is the schematic of a 200-MHz sinusoidal oscillator. The disk 
file ch5-5a.cir contains the listing of the circuit. Simulate this circuit in 
RSPICE with the command: 

SIM CH5-5A.CIR 

Your results should match Fig. 5.5a. 

In this simulation, the default timestep control algorithm (Local 
Truncation Error) allows the step size to become too large to accurately 
simulate the circuit. The local truncation error measurement for sinu¬ 
soidal waveforms is often poor and results in a poor estimate of the 
step size. The Local Truncation Error timestep control algorithm 
allows the timestep to grow up to one-fiftieth of the total transient 
duration (in this case 4nS). For a 200-MHz signal, 4nS is much larger 
than the worst-case Nyquist criterion of 2.5nS. The 200-MHz signal is 
undersampled! 


200 MHz 





200pF 



Figure 5.4 Circuit schematic for 
circuit file ch5-5-Hr 
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Figure 5.5(a) (Reprinted from Successfully Simulating Circuits With SPICE. Used with 
permission.) 


Now simulate the circuit file ch5-5b.cir with the command: 

SIM CH5-5B.CIR 

Your results should match Fig. 5.5 b. In this simulation, the .OPTION 
LVLTIM=1 directs SPICE to use the Iteration-Count timestep control 
algorithm. Notice the improvement in the simulation results. The Iter¬ 
ation-Count timestep control algorithm limits the timestep to no more 
than the print interval TSTEP (in this case .2nS) and, for this circuit, 
eliminates the distortion caused by undersampling. 

Limiting the maximum step size 

Although the default limit for the maximum step size of the Local 
Truncation Error timestep control algorithm is TSTOP/50, the maxi¬ 
mum step size may be manually limited by the user. 

.TRAN TSTEP TSTOP <TSTART> <HMAX> 

The maximum step size of either timestep control algorithm may be 
limited manually with the HMAX parameter of the transient state¬ 
ment. 2 Both timestep control algorithms use either the default maxi - 
mum timestep or HMAX (whichever is smaller) as the maximum 
transient step size. 
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Figure 5.5(b) The Local Truncation Error timestep control algorithm allows the timestep 
to become much larger than the print interval. Switching to the Iteration-Count method 
or manually limiting the HMAX parameter corrects the problem. (.Reprinted from Suc¬ 
cessfully Simulating Circuits With SPICE. Used with permission .) 



Now simulate the circuit file ch5-5c.cir with the command: 

SIM CH5-5C.CIR 

Your results should again match Fig. 5.5 b. In this simulation, the 
(default) Local Truncation Error timestep control algorithm’s maxi¬ 
mum step size is limited to the print interval (HMAX is set to match 
TSTEP) with the command: 

.TRAN .2NS 200NS 0 .2NS . 

Again notice the improvement over the initial simulation result. 
Also notice that if HMAX is to be set, a numeric value must fill the 
TSTART position! 

For those who are not familiar with the TSTART parameter, 
TSTART specifies the transient simulation time where SPICE begins 
printing the output results. TSTART is not where SPICE begins the 
simulation. (Transient analysis always begins at time T=0.) 

Simulate the disk file ch5-5.cir with the command: 


RSPICE CH5-5.CIR 
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Then graph the result with the command: 

RGRAPH CH5- 5.OUT CH5- 6.MAC /MULTI /MARKERS 

Your result should match Fig. 5.6. In this simulation, the oscillator 
circuit is simulated with and without the HMAX limit. The two results 
are superimposed for comparison. Notice how the simulation without 
the step size limit undersamples the corrected output. 

For many simulations, using either the Local Truncation Error 
timestep control with HMAX set to the print interval or using Iteration- 
Count will produce nearly identical results. Some simulations will 
require that HMAX be set to something smaller than the print interval; 
some simulations will produce an accurate result even with HMAX set 
larger than the print interval. A good rule of thumb for most simula¬ 
tions is to set HMAX equal to the print interval. As an error-prevention 
measure, always set the HMAX parameter on the transient statement 
to help minimize timestep control problems. 

Limitations of the breakpoint adjustment 

In addition to adjusting the step size as the circuit activity changes, and 
limiting the step size to a maximum value, SPICE also aligns the solu¬ 
tion timepoints with the breakpoints of the sources. As the transient 
simulation begins to approach a breakpoint, the timestep control algo- 



TIME 

Figure 5.6 Timestep control failure and corrected simulation. (Reprinted from Success¬ 
fully Simulating Circuits With SPICE. Used with permission .) 
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rithms adjust the timestep so the next solution point falls on the next 
breakpoint. To further enhance the accuracy just after the breakpoint, 
the timestep control algorithms automatically reduce the timesteps just 
after the breakpoint in anticipation of a voltage or current transition. 
For many circuits, timestep control with breakpoint adjustment signif¬ 
icantly improves the accuracy of the simulation result. 

While the breakpoint adjustment works well on many circuits, when 
the step size becomes too large the timestep control algorithm may 
miss one or more source breakpoints. Figure 5.7 is the schematic of an 
SR flip-flop and test generators. The disk circuit ch5-8a.cir contains 
the listing of the circuit. The pulse generator uses short pulses to set 
and reset the flip-flop. Simulate the circuit with the command: 

SIM CH5-8A.CIR 

Your results should match Fig. 5.8a. The pulse being used to set and 
reset the latch is illustrated in Fig. 5.9. 

In this simulation, the breakpoints of the short pulses are skipped by 
the timestep control algorithm! Figure 5.9 shows that the set and reset 
pulses last only 750nS from start of the rising edge to the end of the 
falling edge. Because of the short duration of the set and reset pulse 
and the large timestep, the timestep control algorithm misses the 
breakpoints and completely skips over the set and reset pulses. As a 
result, the flip-flop never receives the set or reset signal. 

This is a fairly common occurrence in simulations which contain 
short pulses. Short is a relative term. A more general description is: 
Any circuit with a voltage or current pulse with a total duration 
of less than Koo of the simulation’s maximum timestep may be 
subject to problems with the breakpoint adjustment during transient 
simulation. 



Figure 5.7 Circuit schematic for circuit file ch5-8.cir. 
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Figure 5.8(a) In this simulation, the timestep control algorithm skipped several source 
breakpoints. (Reprinted from Successfully Simulating Circuits With SPICE. Used with 
permission .) 
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Figure 5.8(6) The breakpoint detection algorithm skipped the short-duration pulses. 
For these types of circuits, an accurate simulation may require limiting the HMAX 
parameter. (Reprinted from Successfully Simulating Circuits With SPICE. Used with 
permission .) 

























Timestep Control 143 


1.250 

R 1.125 

A 

N 1.000 

^ BIS. 0m 

E 150.0m 

N 

T B25 .0m 

5 500.0m 

£ 315.0m 

E 250.0m 

P 

125. 0m 
0.000 

-125.0m 

-625 .0u 0.000 625 .0u 1.250m 1.015m 2.500m 3.125m 3.150m 4.315m 5.000m 5.625m 

TIME 

Figure 5.9 The sync pulse is 750nS from the start of the leading edge to the end of the 
falling edge. With a transient duration in the microsecond range, this small sync pulse 
could easily be missed by the breakpoint detection algorithm. For these types of circuits, 
setting the maximum timestep is crucial. 


SYNC PULSE OF CH5 -0 CIRCUIT 


. V(10 ) 



..... 



For circuits which use short pulses, set the maximum timestep to no 
more than 50 times the minimum pulse duration. This will help ensure 
that the breakpoint adjustment algorithm will detect the presence of 
the source breakpoint and adjust the timestep accordingly. 

In the case of the ch5-8 circuits, the pulse duration is 750nS (250nS 
rise + 250nS fall + 250nS duration). The above rule states the maxi¬ 
mum timestep for this circuit should be: 

50 * 750nS = 37.5uS 

The circuit file ch5-8b.cir is identical to the original circuit file except 
for the addition of the HMAX parameter on the transient statement. In 
the transient statement, HMAX is set to 37.5uS. Simulate the ch5- 
8b.cir circuit file with the command: 

SIM CH5-8B.CIR 

Your output should match Figure 5.86. 

An astute reader may also recognize that long pulse waveforms with 
short rise and fall times might also be subject to breakpoint problems. 
While this is true, pulses which have short rise and fall times but a 
long pulse duration create less problems than short pulse durations. 
Imagine a long pulse with a very short rise time. Although SPICE may 
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miss the rising edge source breakpoint, if the pulse duration is long the 
next timepoint will occur with source at its pulsed value. So while the 
timing might be off slightly, the simulation does not entirely skip the 
pulse. For those circuits where the timing of the rising and falling edge 
is critical, set the maximum timestep to no more than 50 * (rise time) 
or 50 * (fall time) to ensure the breakpoints of the rising and falling 
edges are found. 

Timestep Control Settings 
for Typical Circuits 

The previous circuit examples illustrate the common problems with the 
SPICE timestep control algorithms. Table 5.1 summarizes the differ¬ 
ences between the timestep control algorithms in SPICE. Table 5.2 pre¬ 
scribes a general set of guidelines which should be used to select the 
timestep control algorithm and the maximum step size. Using Table 5.2 
to select a timestep control algorithm and maximum timestep will elim¬ 
inate most of the problems caused by the timestep control algorithms. 


Timestep Control and Oscillator Circuits 

Oscillator circuits are a special class of circuits that present unique 
problems for both timestep control algorithms. This is especially true 
in the simulation of the oscillator during the start-up period. 

Figure 5.10 is a simple Colpits oscillator. The disk file ch5-lla.cir is 
a listing of the oscillator. Simulate this circuit with the command: 

SIM CH5-11A.CIR 

Your results should match Fig. 5.11a. 


TABLE 5.1 SPICE Timestep Control Algorithm Summary 


Characteristic 

Iteration count 

Local truncation error 

Initial timestep 

TSTOP/50 

TSTOP/50 

Minimum timestep 

TSTOP/50e9 

TSTOP/50e9 

Maximum timestep 

TSTEP 

TSTOP/50* 

(without HMAX) 

Maximum timestep 

MIN(HMAX,TSTEP) 

MIN(HMAX,TSTOP/50) 51 

(with HMAX) 

.OPTION ITL3=X 

Iteration limit to increase 

Not used 

(default=4) 

timestep 


.OPTION ITL4=X 

Iteration limit to reduce 

Iteration limit to reduce 

(default=10) 

timestep 

timestep 


*If circuit contains no energy storage elements (capacitors or inductors), the maximun 
timestep for LTE without HMAX is TSTEP, and the maximum timestep with HMAX is 
MIN(HMAX,TSTEP). 
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table 5.2 General Timestep Control Algorithm Settings with Maximum Timestep 
Recommendation 


Nature of circuit 

Timestep control algorithm 

General circuits 

Sinusoidal with frequency = FREQ 

High inductive content 

Pulse-driven circuit 

Local Truncation Error w/HMAX = TSTEP 
Local Truncation Error w/HMAX = TSTEP or 
w/HMAX = 1/(8*FREQ) 

Iteration-Count or Local Truncation Error 
w/HMAX = TSTEP 

Iteration-Count w/HMAX = 50 * Shortest 
Pulse Duration or Local Truncation Error 
w/HMAX = 50 * Shortest Pulse Duration 


In this circuit, SPICE fails to predict the start-up of the oscillator. 
The reason SPICE fails to predict the oscillation relates to the opera¬ 
tion of the timestep control algorithm. 

Most real oscillators start oscillating because of minute voltage 
and current surges (noise) within the circuit. SPICE also generates 
minute amounts of voltage and current noise in the form of round-off 
error. The error tolerance options of the program (RELTOL, VNTOL, 
and ABSTOL) determine the magnitude of the round-off error noise. So 
why doesn’t the oscillator start? The oscillator fails to start because the 
voltage surges for the size of the timestep being used were not large 
enough to unsettle the circuit. 

Assume the voltage error tolerance of the program is set to luV 
(.OPTION VNTOL=luV). The minimum amount of numerical noise at 
each timepoint is ±luV. The voltage surge at each timepoint is 
delta__V/delta_T or luV/timestep. As the transient simulation begins, 
SPICE computes an output voltage very close to zero volts. Oscillation 



Figure 5.10 Circuit schematic for circuit file ch5-ll.cir. 
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Figure 5.11(a) (Reprinted from Successfully Simulating Circuits With SPICE. Used with 
permission .) 


has not yet begun. As the transient simulation continues, the step size 
grows quickly because the timestep control algorithm senses no volt¬ 
age change at the output. As the step size continues to grow, the value 
of the numerical noise ( dV/dt ) decreases because the luV voltage is 
now being divided by larger and larger timesteps. The cycle continues. 
The timestep control algorithm increases the timestep because oscilla¬ 
tion has not started, and oscillation does not start because the amount 
of numerical noise decreases with each larger timestep. 

Often, after just a few transient timepoints, the timestep is too large 
to produce a sufficient dV/dt to start the oscillator and, because of this, 
the oscillator never oscillates. 

To correct this situation, when simulating oscillators during start¬ 
up, the maximum timestep must be set to a small value, small enough 
so the numerical noise can produce a sufficient dV/dt to start the oscil¬ 
lator. A good rule of thumb for oscillators is to set the maximum time- 
step to % of the period of oscillation. 

The disk file ch5-llb.cir is the same Colpits oscillator with the max¬ 
imum timestep (HMAX) set at 1/(8 * FREQ), where FREQ is the natu¬ 
ral frequency of oscillation. Simulate this circuit with the command: 

SIM CH5-11B.CIR 


Your results should match Fig. 5.116. 
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Figure 5.11(b) Simulating oscillators may cause timestep control problems. Often, accu¬ 
rate simulation of these circuits requires limiting the HMAX parameter. C Reprinted from 
Successfully Simulating Circuits With SPICE. Used with permission.) 


CORRECTED OSCILLATOR SIMULATION 



Setting the maximum timestep allows this simulation to correctly 
predict the behavior of the oscillator. The maximum timestep required 
to start oscillation will vary from circuit to circuit, but the Vsth rule is 
generally a good starting value. 

Even after setting the maximum timestep, some oscillator circuits 
will still fail to start. For these circuits, ramping the power supply from 
zero up to full power often results in successful start-up. Figure 5.12 
illustrates a CMOS ring oscillator. Simulate the circuit disk circuit 
ch5-13a.cir with the command: 

SIM CH5-13A.CIR 

Your results should match Fig. 5.13a. 




-OVout 


Figure 5.12 Circuit schematic for circuit file ch5-13.cir. 
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Figure 5.13(a) (Reprinted from Successfully Simulating Circuits With SPICE. Used with 
permission .) 


OSCILLATOR SIMULATION FAILURE 

V (2) 


The ch5-13a.cir circuit file is an inverter oscillator (ring oscillator). 
The circuit uses a series combination of seven CMOS inverters. Notice 
the maximum timestep (HMAX) has manually been set to Vs the period 
of oscillation in the circuit file and yet the oscillator still fails to oscillate. 

In simulating ch5-13a.cir, the circuit power supply is set to 5VDC. In 
the circuit file ch5-13b.cir, the power supply will be ramped with the 
aid of a PWL source description. In this simulation, the power supply 
transits from zero to 5VDC much like applying power to a dead circuit. 
Simulate the ch5-13b.cir circuit with the command: 

SIM CH5-13B.CIR 

Your results should match Fig. 5.136. 

On occasion, limiting the maximum timestep and ramping the power 
supply will still fail to start some stubborn oscillators. One sure way of 
starting the most stubborn oscillator circuit is to use the .IC statement 
to initialize one of the nodes within the oscillator’s loop. But this is an 
artificial means of initiating oscillation and will not represent the real 
start-up process of the oscillator. 

Oscillators by their very nature can lead to timestep control prob¬ 
lems and simulation failure. Setting the maximum timestep to an 
appropriate value, ramping the power supplies, or initializing nodes 
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Figure 5.13 (b) In addition to limiting the HMAX parameter, some oscillators will require 
ramping the power supplies to simulate the onset of oscillation. (Reprinted from Suc¬ 
cessfully Simulating Circuits With SPICE. Used with permission.) 


CORRECTED OSCILLATOR SIMULATION 



with the .IC statements are all important techniques for transient sim¬ 
ulation of oscillator circuits. 

Interpolated vs. Noninterpolated 
Output Results 

In Berkeley SPICE2G.6, the transient simulation output is interpo¬ 
lated from the actual timepoints SPICE computes and the print inter¬ 
val found on the transient statement (Fig. 5.16). Berkeley SPICE2G.6 
does not give users access to the actual timepoints where SPICE solves 
the circuit equations. Because of this, SPICE users must either set the 
print interval to a small value (which results in large output files) or 
set the print interval to a large value and risk the chance of missing 
parts of the transient response. 

In most of the vendor-offered, SPICE-like simulators, users are 
given the choice of printing either the actual (noninterpolated) time- 
points where SPICE solves the circuit or the interpolated output 
results. It is highly recommended that SPICE users view the noninter¬ 
polated transient results whenever possible. Simulate the disk circuit 
ch5-14a.cir with the command: 


SIM CH5-14A.CIR 
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Your results should match Fig. 5.14a. Now simulate the disk circuit 
ch5-14b.cir with the command: 

SIM CH5-14A.CIR 

Your results should match Fig. 5.146. 

The disk file ch5-14b.cir contains a listing identical to the ch5-14a.cir 
circuit with the addition of the RSPICE specific option NOINTR 
(.OPTION NOINTR). The NOINTR (NO INTeRpolate) option sup¬ 
presses the interpolation routine and prints the actual timepoints at 
which RSPICE solves the circuit equations. Users should note that the 
NOINTR is not a standard SPICE option and each of the SPICE-like 
simulators may use a different option or command to suppress the 
interpolation routine. 

Timestep Control in Other Simulators 
Hspice 

Hspice uses a proprietary timestep control algorithm named DVDT. 
This algorithm senses the rate of change of voltage levels in the circuit 
and adjusts the timestep accordingly. Hspice also offers the standard 
Local Truncation Error and Iteration-Count timestep control algo¬ 
rithms. 
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Figure 5.14(a) (Reprinted from Successfully Simulating Circuits With SPICE. Used with 
permission .) 
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Figure 5.14(d) Interpolation often hides the finer details of a simulation. Use the actual 
solution points (noninterpolated data) whenever possible. (.Reprinted from Successfully 
Simulating Circuits With SPICE. Used with permission .) 


All three algorithms benefit from the following enhancements. First, 
Hspice automatically limits the maximum timestep to RMAX times the 
print interval. RMAX is a user-selectable option which has a default 
value of 2. Second, an improved breakpoint adjustment decreases the 
chance that Hspice will skip a source breakpoint. 

In Hspice, the POST option forces the simulator to store the nonin¬ 
terpolated results in a raw data file for later viewing. If the option 
INTRPL is set on the option statement, Hspice saves the interpolated 
results in the raw data file. The transient results printed in the output 
file are the interpolated values of the output. 

IS_Spice 

IS_Spice uses the standard SPICE2G.6 timestep control algorithms. 
Both the Iteration-Count and Local Truncation Error methods are 
available. The Local Truncation Error algorithm is the default method 
in IS_Spice. 

IS_Spice does not offer a mechanism to generate the noninterpolated 
output results during transient analysis. 

Micro-Cap IV 

Micro-Cap IV uses a modification of the Local Truncation Error time- 
step control algorithm. This algorithm is fundamentally identical to 
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SPICE’s algorithm, with the addition of an enhancement which 
improves the algorithm’s ability to detect source breakpoints. This 
improvement helps ensure that Micro-Cap IV will not skip source 
breakpoints. 

Because of this improvement, Micro-Cap IV does not offer the 
Iteration-Count timestep control algorithm as a user-selectable 
option. 

In Micro-Cap IV, the simulator saves and graphs the noninter- 
polated transient results. Micro-Cap IV does not generate a standard 
SPICE output file. 

Pspice 

Pspice uses a modification of the Local Truncation Error timestep 
control algorithm. This algorithm is fundamentally identical to 
SPICE’s algorithm, with the addition of two enhancements which 
improve the algorithm’s performance. The first improvement is an 
automatic limit to the maximum timestep whenever a sine wave 
source is in use. In Pspice, the maximum timestep is automatically 
limited to 1/(8*FREQ), where FREQ is the highest frequency of any 
sinusoidal generator in the circuit. The second improvement is a 
change in the breakpoint calculation, which ensures that Pspice will 
not skip source breakpoints. 

Because of these improvements, Pspice does not offer the Iteration- 
Count timestep control algorithm as a user-selectable option. 

In Pspice, the .PROBE statement forces the simulator to store the 
noninterpolated transient results in a raw data file for later viewing. 
The transient results printed in the output file are the interpolated 
results of standard SPICE. 

Summary 

The dynamic timestep control algorithms greatly improve the speed 
and accuracy of the transient simulation. But the timestep control 
algorithms are not foolproof, especially the (default) Local Truncation 
Error algorithm. Too often, the timestep becomes too large for a given 
circuit response. The resulting simulation often contains a large per¬ 
centage of error. Sometimes, the timestep will become large enough to 
completely skip important behaviors of the circuit. 

To guarantee high-quality simulation results, SPICE users need to 
know which timestep control algorithm to use, how to manually limit 
the maximum timestep with the HMAX parameter on the transient 
statement, and how to view the noninterpolated transient output 
results when available. 



Timestep Control 153 


References 

1. L. W. Nagel, SPICE2: A Computer Program to Simulate Semiconductor Circuits, Elec¬ 
tronics Research Laboratory Rep. No. ERL-M520, University of California, Berkeley, 
1975. 

2. SPICE Version 2G.6 User's Guide. 



Chapter 

6 

SPICE Options 


One of the most important commands in the SPICE input file is the 
.OPTIONS statement. The .OPTIONS statement is used to set param¬ 
eter values which control the simulation. The parameters alter the 
numeric algorithms of SPICE, change the way the simulation results 
appear in the output file, generate additional information about the 
simulation or the netlist, as well as perform a variety of miscellaneous 
functions. 

Understanding the .OPTIONS statement parameters is crucial to 
producing fast, accurate, convergent SPICE simulations. While the 
default parameter settings produce accurate results on some circuits, 
many others will require adjustments to one or more parameters. 
Remember SPICE was originally optimized for integrated circuit-level 
voltages and currents. Board- and system-level designers and power 
supply designers have already seen several parameter values which 
must be changed simply because of the magnitude of voltages and cur¬ 
rents flowing through their circuits. Circuits which exhibit noncon- 
vergent behavior can be forced to a convergent state by altering the 
parameters which control the numeric algorithms of SPICE. Tran¬ 
sient simulations which display timestep control anomalies or 
numeric integration instabilities can be corrected by applying the 
proper parameter settings. Chapters 3, 4, and 5 all focused on produc¬ 
ing faster, more accurate, convergent simulations. In these chapters, 
specific problems were corrected with the application of one or more 
•OPTIONS parameters. The .OPTIONS statement serves as an 
avenue to both correcting circuit simulation problems and optimizing 
simulation run time. For fast, accurate, convergent simulations, 
understanding the .OPTIONS statement is not just desirable, it is a 
requirement. 
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SPICE has 34 user-selectable options. Most of the vendor-offered 
SPICE-based simulators use the same options and often add many new 
ones. Of the 34 available, 33 are documented in the Berkeley SPICE2G.6 
User’s Guide. 1 One was never added to the original documentation. Of 
these 34,6 operate as flags and signal an event when the option is named 
on the .OPTIONS line. As an example, to print accounting information 
(run time, number of iterations used, etc.) for your simulation run, a user 
would specify .OPTIONS ACCT. The other 29 operate as parameters and 
require that a numeric value be assigned on the .OPTIONS line. For 
example, to set the number of iterations allowed during the DC bias- 
point calculation, a user would specify .OPTIONS ITL1=X, where X is the 
number of iterations allowed. 

The proper command syntax for the .OPTIONS statement is: 

.OPTIONS <optionl option2 option3 ...> 

However, because of the way SPICE processes the lines of the input 
file, only the first three characters in the .OPTIONS command are 
required. All of the following are acceptable forms of the .OPTIONS 
statement. 

.OPTIONS <optionl option2 option3 ...> 

.OPTION <optionl option2 option3 ...> 

.OPTIO <optionl option2 option3 ...> 

.OPT I <optionl option2 option3 ...> 

.OPT <optionl option2 option3 ...> 

Users may specify more than one .OPTIONS line in an input file. If 
the same option is declared more than once, only the last value of the 
option will be used. This is important in simulations which are com¬ 
posed of different blocks or subsections, each of which may contain dif¬ 
ferent option statements. SPICE does not allow different sections of a 
netlist to operate with specific option settings. Once set, an option 
applies globally to the circuit until another .OPTIONS statement 
resets the value. 

All of the options have names composed of two to six alphanumeric 
characters. Some of the names were well chosen and clearly express 
the meaning of the option. For example, the option which selects the 
numeric integration method used in the simulation, METHOD= 
‘GEAR’, is easy to understand and use. But often the names are less 
than self-explanatory. For example, the option which selects which 
timestep control algorithm is used in the simulation is LVLTIM=X. 
Lastly, there are just a few options which have names which appear to 
be explicit and meaningful but do not operate as the name implies. For 
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example, RELTOL, ABSTOL, and VNTOL all determine the accuracy 
of the simulation results, as might be implied from the TOLerance por¬ 
tion of the name. The options CHGTOL and TRTOL are both correction 
factors for a shortcoming in the Local Truncation Error timestep con¬ 
trol algorithm. They do not have a direct effect on the accuracy of the 
simulation results even though they both have TOL in their names. A 
clear understanding of each of the options helps build a strong founda¬ 
tion for high-quality simulation. 

As a fledgling SPICE user, the author recalls trying to decode the 
mysterious one-line description of each option in the SPICE User’s 
Guide. 1 Since this young engineer had little idea of what many of the 
options did, many misbehaving simulations were often liberally 
doused with a variety of different options. With each simulation, the 
young engineer anxiously awaited the results, each time hoping that 
these option settings would be the key to success. Often, hours of sim¬ 
ulation time were consumed before stumbling on the right combina¬ 
tion. Too many times, the right combination was never found. 

SPICE Option Definitions 

In this section, each of the 34 standard SPICE options will be defined 
in clear detail. Where applicable, formulas showing how each of the 
option values should be set for specific circuits will be given. Some 
options should always be set regardless of the type of simulation. Oth¬ 
ers may only be needed occasionally. With each option, recommenda¬ 
tions will be made as to whether the option should always be set and/or 
when it should be used. 

For designers using one of the SPICE-like simulators, most of the 
options detailed here have a one-to-one correspondence with the same 
named option in your simulator. Appendix A lists the .OPTIONS 
parameters for several popular circuit simulators and, where possible, 
shows the correspondence between the SPICE2G.6 option and the 
vendor-offered option. 

The first 8 options may be classified as options which control the out¬ 
put file format, the next 18 control the numeric algorithms of SPICE, 
4 determine the default MOS gate geometries, and finally, 4 more can 
only be classified as miscellaneous options. The options will be pre¬ 
sented in this order. 

Output File Format Options 

The output file format options include ACCT, LIST, NOMOD, 
NOPAGE, NODE, OPTS, LIMPTS=X, and NUMDGT=X. 



158 Chapter Six 


.OPTIONS ACCT 

The accounting option is a flag and directs SPICE to print 47 run¬ 
time statistics about the job just simulated. The results appear in the 
output file. The accounting section of the output file contains an 
abundance of information about the simulation run. Within the 
accounting section, users will find information relating to the speed 
and efficiency of the simulation, the time and number of iterations 
each analysis type consumes, the size and number of elements in the 
circuit, the size of the solution matrix, the number of analysis points, 
and the maximum number of nodes or circuit elements which may be 
simulated. 

Figure 6.1 shows an example of the accounting information as it 
appears in the output file. The information is segmented into two dif¬ 
ferent forms, numeric values pertaining to elements of the simulation 
and run-time statistics. 

The first row of values in the accounting section relates to the ele¬ 
ments in the circuit. The second row displays information regarding 
the analysis. The third reveals details of the system equations, and the 
final row contains data about the transient analysis and the memory 
requirements of the simulation. 

The headings which describe the circuit element, analysis, system 
equation, and transient and memory entries are shown in Tables 6.1a, 
b, c, and d. The items marked in italics will be the most useful to sim¬ 
ulation users. 

Following the four rows of data are four columns of data displaying 
the run-time statistics and number of iterations each program segment 
used. The program segment headings are shown in Table 6.2. Again, 
the headings shown in italics will be of interest to most users. 

Two of the best uses for the information found in the accounting sec¬ 
tion are to improve the speed or efficiency of a simulation and to deter¬ 
mine the maximum circuit size or maximum analysis time supported 
by your simulator. 

Using the accounting information to improve simulation. To improve the 
speed of simulation, observe the number of iterations each of the major 
analysis types consumes and use this as a gauge as you alter one or 
more options or commands of the simulation. For example, Table 6.3 
shows the number of iterations required and the amount of time 
required to compute the dc bias point of disk circuit ch6-3.cir. For this 
circuit, adding the .NODESET statement reduced the number of itera¬ 
tions from 53 to 13 and reduced the time used for the bias point analy¬ 
sis by 56 percent. Try adding one or more .NODESETs to your circuit 
and watch the number of iterations listed under the DCAN program 
segment. With the proper settings, you can often reduce the number of 
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Figure 6.1 ACCT option output. 
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table 6.1(a) Circuit Element Information 


NUNODS 

NCNODS 

NUMNOD 

NUMEL 

DIODES 

BJTS 

JFETS 

MFETS 


This is the number of nodes in the main circuit. 

This is the number of circuit nodes after subcircuit expansion. 

This is the number of nodes after subcircuit expansion and includes 
the internal parasitic nodes for semiconductor devices. 

This is the number of elements (components) in the circuit. 

This is the number of diodes in the circuit. 

This is the number of BJTs in the circuit. 

This is the number of JFETs in the circuit. 

This is the number of MOSFETs in the circuit. 


TABLE 6.1(b) 

NUMTEM 

ICVFLG 

JTRFLG 

JACFLG 

INOISE 

IDIST 

NOGO 


Analysis Information 

This is the number of temperatures at which the circuit will be analyzed. 

This is the number of output points in a DC sweep analysis. 

This is the number of output points in a transient sweep analysis. 

This is the number of output points in an AC frequency sweep. 

This flag indicates whether a noise analysis was performed. 

This flag indicates whether a distortion analysis was performed. 

This flag indicates whether an error condition was encountered during 
simulation. A number of different error conditions may set this flag. 


TABLE 6.1(c) Circuit Element Information 


NSTOP 

NTTBR 

NTTAR 

IFILL 

IOPS 

PERSPA 


This is the number of nodes in the main circuit matrix (N x N). 

This is the number of matrix entries including subcircuit expansion 
and internal nodes. 

This is the number of matrix entries including subcircuit expansion 
and internal nodes after “fill-in” is complete. 

This is the number of “fill-ins” generated in the conductance matrix. 

This is the number of floating-point operations required to solve the 
system equations. 

This is the percent of the conductance array which was empty. Elec¬ 
tronic circuits generally lead to extremely sparse system equations. 
While PERSPA may not be useful directly, it is interesting to observe 
this number for a variety of circuit types. 
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TABLE 6.1(d) 

Transient Analysis and Memory Information 

NUMTTP 

This is the number of internal timepoints SPICE used during the 
transient analysis. This is not the same as JTRFLG, which is the 
number of interpolated output points. 

NUMRTP 

This is the number of times SPICE failed to converge, cut the 
timestep, and reattempted the circuit solution during transient 
analysis. For an efficient run, NUMRTP should be less than 10 per¬ 
cent of NUMTTP. In most simulations, NUMRTP may be reduced 
(transient efficiency increased) by raising ITL4. 

NUMNIT 

This is the total number of iterations used during the transient 
simulation. 

MAXMEM 

This is the maximum amount of memory available to SPICE. The 
units of memory are 4-byte words. 

MEMUSE 

This is the amount of memory used during the simulation. The units 
of memory are 4-byte words. Comparing this number to MAXMEM 
enables users to determine the maximum circuit size or the maxi¬ 
mum run time a particular version of SPICE can support. 

COPYKNT 

This is the number of times the copy routine was called in SPICE. 
The copy routine is used to move numeric values throughout the 
program. 


TABLE 6.2 Program Segment Headings 


READIN 

SETUP 

TRCURV 

DCAN 

DCDCMP 

DCSOL 

ACAN 

TRANAN 

OUTPUT 

LOAD 

CODGEN 

CODEXC 

MACINS 

OVERHEAD 


This segment reads the input file listing, parses the input file syntax, 
and stores the circuit elements in memory. 

This segment defines the system equations and initializes the circuit 
sources. 

This segment performs the DC sweep analysis. 

This segment performs the DC operating point (bias point) calculation. 

This segment performs the numerical pivoting and LU factorization of 
the system equations. 

This segment performs the back substitution of the system equations. 

This segment performs the AC small-signal frequency sweep analysis. 

This segment performs the transient time sweep analysis. 

This segment prints the simulation results to the output file. 

This segment loads the device conductance and current values into the 
solution array during the iterative solution procedure. 

This segment generates machine code to solve the solution array. This 
routine is called with the LVLCOD option. Most versions of SPICE 
have disabled this option. 

This segment solves the solution array when the LVLCOD option is 
selected. Most versions of SPICE have disabled this option. 

This segment builds the machine instructions to solve the solution 
array when the LVLCOD option is selected. Most versions of SPICE 
have disabled this option. 

This is time consumed by SPICE for general housekeeping chores. 
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TABLE 6.3 Bias Point Run-Time Statistics 



Number of iterations 

Time 

Without NODESET 

53 

4.29 sec. 

With NODESET 

17 

1.87 sec. 


iterations by 50 percent or more and eliminate a considerable portion 
of the time required for the operating point calculation. 

For transient analysis, watch the values of NUMTTP, NUMRTP, and 
NUMNIT. NUMTTP is the number of timepoints SPICE uses during 
the transient simulation; this is not the number of points printed in the 
output file. NUMRTP is the number of timepoint reversals or the num¬ 
ber of times SPICE cuts the timestep and steps back in time to reat¬ 
tempt the solution point after failing to converge on a given timepoint. 
NUMNIT is the total number of iterations required to complete the 
transient analysis. By raising the value of ITL4 on the .OPTIONS 
statement, SPICE will cut the timestep less often during transient 
analysis. The less timepoint reversals SPICE uses, the faster the sim¬ 
ulation will complete. For an efficient simulation, the values of 
NUMTTP, NUMRTP, and NUMNIT should be as low as possible. In 
general, the values of NUMTTP, NUMRTP, and NUMNIT may be 
decreased by raising ITL4. As you raise ITL4 in your simulation, watch 
the values of NUMTTP, NUMRTP, NUMNIT, and the amount of time 
used for the transient simulation. Continue raising ITL4 until no fur¬ 
ther improvement in the NUMTTP, NUMRTP, and NUMNIT can be 
achieved, the simulation time ceases to be reduced, or accuracy is 
degraded. 

The disk file ch6-4.cir simulates the transient response of an SR flip- 
flop. In the ch6-4.cir file, the original circuit description has been dupli¬ 
cated several times with different values of ITL4. Table 6.4 shows the 
ITL4 value, the number of transient timepoints, the number of time- 
point reversals, the number of transient iterations, and the overall 
simulation time. Notice that raising ITL4 beyond 40 produces no sig- 


TABLE 6.4 Transient Run-Time Statistics 


ITL4 

NUMTTP 

NUMRTP 

NUMNIT 

Job time 

10 (default) 

257 

38 

1012 

17.47 sec. 

15 

136 

7 

561 

9.77 sec. 

20 

112 

1 

471 

8.29 sec. 

40 

108 

0 

417 

7.53 sec. 

100 

108 

0 

417 

7.53 sec. 


Results obtained from a 33-MHz 386/387 computer. 
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nificant improvement in simulation speed. Simulate this circuit with 
RSPICE to verify these results. 

To determine the maximum number of nodes or elements allowed in 
your SPICE program, watch the MAXMEM and MEMUSE values. 
SPICE uses a large data array to store the circuit description, the ele¬ 
ment values, and all of the model parameters. The same data array 
stores all the analysis results. The size of the array determines the 
maximum number of elements or the maximum number of analysis 
points allowed for a given program. If the storage requirements of your 
simulation ever exceed the amount of memory available, SPICE will 
print a most unfriendly message in the output file and abort the simu¬ 
lation. MAXMEM is the maximum amount of memory your version of 
SPICE can address. MEMUSE is the amount of memory used during a 
simulation. The data array must hold both the elements of a circuit 
and all of the analysis data. Because of this, SPICE can simulate the 
largest number of elements when performing a very short analysis, say 
a dc operating point calculation. Conversely, SPICE can simulate the 
longest analysis (most analysis points) when a very small circuit is 
being simulated. 

To determine the maximum number of elements SPICE can simu¬ 
late, create a circuit netlist with several elements. Include the ACCT 
option in the circuit netlist. Simulate the circuit, record the MAXMEM 
and MEMUSE values, and add several additional elements to the orig¬ 
inal netlist. Again simulate the circuit and record the new MEMUSE 
value. The difference between the two MEMUSE values represents the 
amount of memory required for the additional elements. The amount 
of memory per element and MAXMEM can be used to determine the 
theoretical maximum number of elements SPICE can simulate. The 
disk file named ch6-5.cir contains a 7-stage and a 14-stage ring oscilla¬ 
tor. When simulated, the MEMUSE value can be used to calculate the 
amount of memory required for the seven additional stages of the cir¬ 
cuit. Table 6.5 shows the amount of memory required to simulate each 
circuit. The addition of seven extra stages required 4376 memory units 
(or 625 memory units per inverter stage). By extrapolating the memory 
units per stage to the 200,000 memory units available, the user may 
deduce that up to approximately 320 stages (320 inverters) could be 
used in this analysis! 


table 6.5 Run Statistics for Maximum Element Calculation 


Oscillator stages 

MEMUSE 

MAXMEM 

7 

4720 

200,000 

14 

9096 

200,000 
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Obviously, the maximum number of elements will be reduced when a 
more sophisticated analysis is performed. 

A similar calculation can be made for the maximum simulation 
length of any job. For example, to determine if a long transient simula¬ 
tion may be performed, run a short transient simulation on the circuit. 
Then increase the transient duration by 20 percent over the previous 
run. The difference in MEMUSE values represents the amount of mem¬ 
ory required for the 20 percent increase in transient duration. The max¬ 
imum transient duration can be calculated from the MAXMEM value 
and the amount of memory used for the longer transient run. 

The disk file named ch6-6.cir contains the same 7-stage ring oscilla¬ 
tor of the previous problem. The first circuit file simulates the tran¬ 
sient behavior for 4nS (with a lOOpS print resolution). The second 
simulates the behavior for 8nS (again with a lOOpS print resolution). 
Table 6.6 shows the amount of memory the two runs use. 

From the amount of memory per nanosecond of transient duration, 
the user may deduce that the same lOOpS print resolution for this cir¬ 
cuit will yield a maximum transient duration of 18uS. 

For the purist, both MAXMEM and MEMUSE have units of four- 
byte words; so, to determine the actual number of bytes of memory 
used for either of these, multiply the MAXMEM and MEMUSE values 
by 4. The result is the number of bytes SPICE is using. 


.OPTIONS LIST 

The LIST option is a flag and directs SPICE to print a list of all the ele¬ 
ment types of a simulation in the output. The list is ordered by element 
type. All resistors are listed first, then all capacitors, inductors, etc., 
until all the elements in the circuit have been listed. Included in the 
list are the nodal connections for each element, the element value, any 
optional element parameters, and subcircuit connectivity information. 
Figure 6.2 is an example of the LIST option. 

The LIST option is best used in diagnosing topology-related errors. 
Erroneous node connections can quickly be uncovered and resolved. 
The list is ordered by element type, which makes node connections 
much easier to read. 


TABLE 6.6 Run Statistics for Maximum Run-Time Calculation 


Transient duration 

MEMUSE 

MAXMEM 

4nS 

4720 

200,000 

8nS 

4764 

200,000 




THIS IS A SIMPLE X STAGE RING OSCILLATOR 



**** CIRCUIT ELEMENT SUMMARY 

temperature - 

27.000 DEG C 


**** RESISTORS 





NAME 

NOOES 

VALUE 

TCI 

TC2 


RUP.X1 

1000 2 

1.43E+03 

O.OOE+OO 

O.OOE+OO 


RUP.X2 

1000 3 

1.43E+03 

O.OOE+OO 

O.OOE+OO 


RUP.X3 

1000 4 

1.43E+03 

O.OOE+OO 

O.OOE+OO 


RUP.X4 

1000 5 

1.43E+03 

O.OOE+OO 

O.OOE+OO 


RUP.X5 

1000 6 

1.43E+03 

O.OOE+OO 

O.OOE+OO 


RUP.X6 

1000 7 

1.43E+03 

O.OOE+OO 

O.OOE+OO 


RUP.X7 

1000 1 

1.43E+03 

O.OOE+OO 

O.OOE+OO 


**** CAPACITORS AND INDUCTORS 




NAME 

NOOES 

IN COKO 

VALUE 



COUT.X1 

2 0 

O.OOE+OO 

2.00E-14 



COUT.X2 

3 0 

O.OOE+OO 

2.00E-14 



COUT.X3 

4 0 

O.OOE+OO 

2.00E-14 



C0UT.X4 

5 0 

O.OOE+OO 

2.00E-14 



COUT.X5 

6 0 

O.OOE+OO 

2.00E-14 



C0UT.X6 

7 0 

O.OOE+OO 

2.00E-14 



C0UT.X7 

1 0 

O.OOE+OO 

2.00E-14 



**** INOEPENOENT SOURCES 





NAME 

NOOES 

DC VALUE 

AC VALUE 

AC PHASE 

TRANSIENT 

VDD 

1000 0 

5.00E+OG 

O.OOE+OO 

O.OOE+OO 


**** MOSFETS 






NAME 0 

G S 3 

MODEL 

W AD 

PO 

RDS 




L AS 

PS 

RSS 

Hi .XI 2 

1 0 0 NCKAN 

1.3E-04 O.OE+OO O.OE+OO 

1.0E+00 




3.0E-06 O.OE+OO O.OE+OO 

1.0E+00 

M1.X2 3 

200 NCHAN 

1.3E-04 O.OE+OO O.OE+OO 

1.0E+00 




3.0E-06 O.OE+OO O.OE+OO 

1.0E+00 

M1.X3 4 

300 NCHAN 

1.3E-04 O.OE+OO O.OE+OO 

1.0E+00 




3.0E-06 O.OE+OO O.OE+OO 

1.0E+00 


Figure 6.2 LIST option output. 
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LIST generates a list of components in the circuit netlist, NODE gen¬ 
erates a list of the nodes in the circuit and the elements connected to 
each node. Figure 6.4 shows an example of the node connection infor¬ 
mation printed in the output file. 

Although the NODE option can be used to debug topology-related 
problems, NODE does not describe which terminal of a given compo¬ 
nent is connected to a specific node. For this reason, the LIST option is 
a much better topology diagnostic tool. A better use of the NODE 
option is assisting in the resolution of a dc operating point nonconver¬ 
gence. If a circuit fails to converge during a dc operating point analysis, 
SPICE prints a list of the last iterative node voltage values in the out¬ 
put. Often, one or more of these voltage values are extremely, unrealis¬ 
tically high and obviously wrong. When one or more of these nodes are 
encountered, use the NODE option to generate a list of all the elements 
connected to that node. One or more of those elements is generating 
the nonconvergence problem. Once the problem elements have been 
identified, use the .NODESET or OFF statements presented in Chap. 
3 to correct the nonconvergence problem. 


THIS 1$ A SIMPLE X STAGE RING OSCILLATOR 

**** ELEMENT NOOE TABLE TEMPERATURE ■ 27.000 DEG C 


0 

C0UT.X1 

VDD 

M1.X4 

M1.X7 

C0UT.X2 
Ml .XI 
N1.X4 

C0UT.X3 
Ml.XI 
M1.X5 

COUT.X4 

M1.X2 

M1.X5 

C0UT.X5 

M1.X2 

M1.X6 

C0UT.X6 

M1.X3 

M1.X6 

C0UT.X7 

M1.X3 

M1.X7 

1 

RUP.X7 

C0UT.X7 

Ml .XI 

M1.X7 




2 

RUP.X1 

CQUT.X1 

Ml .XI 

M1.X2 




3 

RUP.X2 

C0UT.X2 

M1.X2 

M1.X3 




A 

RUP.X3 

C0UT.X3 

M1.X3 

M1.X4 




5 

RUP.X4 

C0UT.X4 

M1.X4 

M1.XS 




6 

RUP.X5 

C0UT.X5 

M1.XS 

M1.X6 




7 

RUP.X6 

CQUT.X6 

M1.X6 

M1.X7 




1000 

RUP.X1 

VOO 

RUP.X2 

RUP.X3 

RUP.X4 

RUP.X5 

RUP.X6 

RUP.X7 


Figure 6.4 NODE option output. 
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.OPTIONS OPTS 

The OPTS option is one of the most useful options in SPICE. The OPTS 
option is a flag that directs SPICE to print the value of each of the 
optional parameters during a SPICE simulation. At first, this seems 
quite redundant, because users must specify options on the .OPTIONS 
line. But the OPTS option prints all user-selectable options, not just 
those specified on the .OPTIONS statement. 

One of the best uses of the OPTS option is to determine the default 
setting of all the SPICE options. Determining the default settings of 
any SPICE simulator is extremely important because different simula¬ 
tors will assign different default values that will have a dramatic 
impact on the speed, accuracy, and convergence properties of the simu¬ 
lator. 

To determine the default settings in SPICE, generate a simple 
netlist with a .OPTIONS line which has nothing but the OPTS option 
declared. After simulation, the output file will contain a section similar 
to the list shown in Fig. 6.5. These are the default option settings of 
SPICE. 

If you use the OPTS during a simulation with one or more option 
parameters set on the .OPTIONS line (in addition to the OPTS state¬ 
ment), the declared value, rather than the default value, will be 
printed in the output file. 

.OPTIONS LIMPTS=X 

The LIMPTS=X is the parameter that sets an upper limit on the num¬ 
ber of solution points which can be simulated or printed in the output 
file. LIMPTS is one of several safety nets which were installed in 
SPICE to prevent users from accidentally consuming excessive com¬ 
puter resources. For example, if the following line 

.TRAN .IMS 20MS 

was accidentally entered as 

.TRAN .INS 20MS 

SPICE would try to simulate 200 million data points! 

Excessive run times or output files are not a significant problem for 
today’s single-user computers. Workstation or personal computer users 
can easily abort a long simulation run without significant repercus¬ 
sions. But when SPICE was written, a central computing facility, batch 
jobs, and line printer outputs were used. Often, users knew nothing 
about how their job was running until the printed output was carried 
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****** n/24/90 ******* RSPICE V3.2.0 04/15/92 ****** 6:13 pm ******** 

OPTIONS EXAMPLE 

**** OPTION SUMMARY TEMPERATURE * 27.000 DEG C 

0C ANALYSIS « 


GMIN * 

1.0006-12 

RELT0L * 

1.000E-03 

ASSTOL > 

1.0006-12 

VNT0L » 

1.0006-06 

L VI COO - 

1 

ITL1 « 

100 

ITL2 * 

50 

PIVT0L * 

1.0006*13 

PIVREL - 

1.0006-03 

TRANSIENT ANALYSIS - 

METH00 « 

TRAP 

MAXORD • 

2 

CHGTOL » 

1.000E-14 

TRT0L ■ 

7.0006*00 

LVLTIM » 

2 

MU 

0.500 

ITL3 • 

4 

ITL4 * 

10 

ITL5 ■ 

5000 

MISCELLANEOUS 

- 

LIMPTS * 

201 

LIMTIM ■ 

2 

CPTIME * 

100000000 

NUMOGT « 

4 

TNOM ■ 

27.000 

DEFL » 

1.0006-04 

DEFW > 

1.0006-04 

DEFAO • 

0.0006*00 

OEFAS * 

0.0006*00 


Figure 6.5 OPTS option output. 


to the output room. The second .TRAN statement generates 200 mil¬ 
lion output points in the output file. Assuming 66 lines per page and 
100 pages in a one-inch stack of fanfold printer paper, 200 million 
points translates into a stack of printer paper just short of a half mile 
high! Just imagine the look on the computer center person that tried to 
take that stack off the printer! 
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By default, LIMPTS is set to 201 in SPICE. Many vendor-offered 
versions of SPICE reset this to a higher value. During the reading of 
the input file, SPICE determines the number of points the user 
requested for the simulation. If the number of points is greater than 
LIMPTS, SPICE prints an error message in the output file and termi¬ 
nates without performing the simulation. The user can reset the 
parameter to either the exact number of points he or she wishes to sim¬ 
ulate or simply to a larger default value. (The author usually sets 
LIMPTS to 10,000 in most simulations.) 

.OPTIONS NUMDGT=X 

The last of the output options is the NUMDGT=X parameter. 
NUMDGT is the number of significant digits for voltages and currents 
printed in the output file. Changing this option alters the print format 
used by SPICE, not the accuracy of the simulator. By default, 
NUMDGT is set to 4. Allowable values for this parameter are 1 to 7. 


Output File Format Options Summary 

The output file format options determine the type and form of infor¬ 
mation printed in the output file. The information obtained may be 
used to debug netlist-related problems, minimize output file size, 
determine the default SPICE option settings, and determine the max¬ 
imum circuit size or run time. 


Numeric Control Options 

The numeric control options consist of four subgroups. In the first 
group are the convergence parameters GMIN, RELTOL, ABSTOL, and 
VNTOL. The second group contains the iteration limits ITL1, ITL2, 
ITL3, ITL4, ITL5, and ITL6. The transient analysis options group fol¬ 
lows with LVLTIM, METHOD, MAXORD, MU, TRTOL, and CHGTOL. 
Finally, two options define the matrix pivoting group, PIVTOL and 
PIVREL. All of the numeric control options are parameterized options 
rather than flags. Parameterized simply means that the value of the 
option must follow the option name. 


Convergence options 

The convergence options include GMIN, RELTOL, ABSTOL, and 
VNTOL. Chapter 3 describes each of these and discusses how they 
impact the accuracy and convergence characteristics of a circuit. 
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.OPTIONS GMIN=X. The parameter GMIN=X represents the minimum 
allowed conductance (maximum resistance) of any element in the cir¬ 
cuit. GMIN is also the value of small conductance (large resistance) 
placed in parallel with every semiconductor pn junction in the circuit. 
GMIN can be equated with the parasitic leakages associated with a 
given circuit. The default value of GMIN = le - 12 mhos. 

The best use for GMIN is in aiding the convergence properties of the 
simulator. For many circuits, raising GMIN from the default will result 
in better overall convergence characteristics without sacrificing the 
accuracy of the result. To select a value for GMIN, determine the small¬ 
est parasitic resistance which could be placed across any two nodes of 
a circuit without affecting the normal circuit operation. Set GMIN=X to 
the inverse of the resistance value. 

GMIN =-1- 

(smallest parasitic resistance) 

•OPTIONS RELTOL=X. The RELTOL=X parameter is the relative error 
tolerance required for convergence. To solve the circuit equations, 
SPICE starts the Newton-Raphson iterative solution algorithm. The 
Newton-Raphson algorithm iterates over and over, searching for a set 
of node voltages and branch currents which will make the remaining 
circuit equations conform to Kirchhoff’s voltage and current laws. 
Since the Newton-Raphson algorithm never really “knows” the proper 
circuit voltages and currents, the signal to stop iterating must come 
from outside the algorithm. One of the useful properties of the Newton- 
Raphson algorithm is that when the algorithm is close to the exact 
solution, the change in iteration-to-iteration voltage and current 
values approaches zero. 

Because digital computers may have a difficult time predicting when 
the difference in voltage and current is exactly zero (because of round¬ 
off error), SPICE presumes the Newton-Raphson algorithm has found 
the exact solution when the difference between iterative voltages and 
currents is less than a given error tolerance. The first part of the error 
tolerance is a percentage change between iterative voltages and cur¬ 
rent. SPICE stops the iterative procedure when the percentage change 
between iterative voltage values for each node in the circuit is less 
than RELTOL and when the percentage change between iterative 
semiconductor branch currents is less than RELTOL. 

RELTOL=X sets the relative error tolerance allowed in SPICE. The 
default value of RELTOL = .001, or one-tenth of one percent. For most 
circuits this is a good compromise between speed and accuracy. More 
accuracy requires more iterations at a given solution point. More iter- 
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ations translate into more simulation time. Less accuracy requires less 
iterations, resulting in faster simulation execution. Nagel demon¬ 
strated that for each additional significant digit of accuracy required, 
the number of iterations required to achieve the solution doubles. 2 For 
most analog circuits, the default value of RELTOL yields satisfactory 
simulation accuracy and speed. For many digital circuits, RELTOL 
may be raised to .01 or .05 for faster simulation execution without a 
significant change in accuracy, Experiment with RELTOL for different 
circuits. Use a value which results in the best compromise between 
speed and accuracy. The value of RELTOL must be determined before 
the VNTOL=X and ABSTOL=X parameters may be set. 

.OPTIONS VNTOL=X ABSTOL=X. The VNTOL=X and ABSTOL=X 
parameters work with the RELTOL option to determine the error tol¬ 
erance for the iterative solution algorithm in SPICE. The percentage 
change criterion works well for most voltage and current values used 
in circuit simulation, except one. When the node voltage or branch cur¬ 
rent approaches zero, the percentage change error tolerance criterion 
fails to predict a good termination point for the iterative algorithm. 
Because of this, in addition to the percentage change criterion, a lower 
absolute error tolerance must also be determined. The absolute error 
tolerances will define convergence when the node voltage or branch 
current falls to zero. For node voltages in the circuit, VNTOL=X is 
the absolute error tolerance between iterations, and ABSTOL=X is the 
absolute error tolerance for semiconductor branch currents. These two 
option parameters determine the lower resolution on voltage and 
current for the simulator. 

The default value for VNTOL is luV. The default value for ABSTOL 
is IpA. The authors of SPICE set these to integrated circuit level volt¬ 
ages and currents. For any circuit which is not an integrated circuit 
(i.e., board level, discrete, and especially power circuits), the VNTOL 
and ABSTOL parameters should be reset to align with the voltage and 
current levels seen in the circuit. 

To select values for VNTOL and ABSTOL with respect to a given cir¬ 
cuit, study the circuit. Determine the lowest voltage magnitude in the 
circuit and multiply RELTOL by this value. The result is the setting 
for VNTOL. For ABSTOL, determine the lowest current magnitude in 
the circuit and multiply RELTOL by this value. The result is the set¬ 
ting for ABSTOL. These two equations align the relative and absolute 
error tolerances for SPICE. Users should notice that the values for 
VNTOL and ABSTOL are dependent on the value of RELTOL. Any 
modification to RELTOL should be followed with an alteration to 
ABSTOL and VNTOL. 
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Readers are reminded that all of the convergence parameters are 
discussed in detail in Chap. 3. 

Iteration limit options 

The next set of options are the iteration limits. Six parameters define 
the iteration limits for different analysis types. The iteration limits all 
begin with the three characters, ITL. 

.OPTIONS ITL1=X. The ITL1=X parameter defines the upper limit on 
the number of iterations allowed to compute the dc operating point. 
The default value of ITL1 is 100. Empirically it can be shown that 
about 60 percent of all circuits converge within 100 iterations, 75 per¬ 
cent converge within 200 iterations, and 92 percent converge within 
500 iterations. If more than 500 iterations are required for convergence, 
something in the circuit has probably been incorrectly connected or 
one or more nodes needs to be initialized with the .NODESET or .IC 
statements. 

The best use for ITL1 is in improving the dc operating point conver¬ 
gence characteristics. For all circuits, start the simulation with ITL1 
set to 500. Since ITL1 is the upper limit on the allowed number of iter¬ 
ations, simulations which do not require excessive iterations to deter¬ 
mine the dc operating point are not affected by raising ITL1. 

The number of iterations SPICE uses in computing the dc operating 
point, dc sweep, or transient analysis can be observed by using the 
ACCT option discussed earlier in this chapter. 

.OPTIONS ITL2=X. The ITL2=X parameter defines the upper limit on 
the number of iterations allowed at each step in a dc sweep analysis. 
The default value of ITL2 is 50. 

The best use for ITL2 is in aiding convergence during a dc sweep 
analysis. For circuits with high-gain switch points, such as flip-flops, 
comparators, op-amps, and triggers, set ITL2 to 200. 

.OPTIONS ITL3=X. The ITL3=X parameter sets the lower limit on the 
number of iterations for the Iteration-Count timestep control algo¬ 
rithm. When the Iteration-Count method is being used in a transient 
simulation, SPICE monitors the number of iterations at each time- 
point in the simulation. If the number of iterations is ever less than or 
equal to ITL3, the timestep control algorithm automatically doubles 
the step size before the next timepoint is calculated. 

The default value for ITL3 is 4 iterations and was determined empir¬ 
ically. 2 Neither Nagel nor the author found any significant improve¬ 
ment for the Iteration-Count timestep control algorithm by either 



SPICE Options 175 


raising or lowering this value. For most circuits, ITL3 should remain at 
the default value. 

.OPTIONS ITL4=X. The ITL4=X parameter sets the upper limit on iter¬ 
ations for either the Iteration-Count timestep control algorithm or the 
Local Truncation Error timestep control algorithm. If the number of 
iterations at a given timepoint is greater than ITL4, SPICE discards 
the current timepoint, cuts the timestep by a factor of 8, then reat¬ 
tempts the solution at the new timepoint. The default value of ITL4 is 
10 iterations. 

The best use for ITL4 is enhancing the convergence characteristics 
and increasing the speed of a transient simulation. ITL4 has a dra¬ 
matic effect on transient simulations. By raising ITL4, you raise the 
number of iterations allowed before SPICE discards the current time- 
point. Because of the asymmetric manner in which SPICE changes the 
size of the timestep, reducing the number of times the timestep is cut 
drastically increases simulation speed. The section on transient non¬ 
convergence in Chap. 3 explains more about why raising ITL4 will 
improve both the speed and the convergence characteristics of a tran¬ 
sient simulation. 

For transient simulation, set ITL4 to 40. Use the values of NUMTTP, 
NUMRTP, and NUMNIT (see the ACCT option section for a discussion 
of NUMTTP, NUMRTP, and NUMNIT) and adjust ITL4 accordingly for 
improved speed and efficiency in transient simulation. 

.OPTIONS ITL5=X. The ITL5=X parameter is another safety net placed 
in SPICE to limit the total number of iterations used for transient sim¬ 
ulations. At each timepoint in a transient run, SPICE counts the num¬ 
ber of iterations required for solution at that timepoint. If the total 
number ever exceeds ITL5, SPICE will stop the transient simulation, 
print the result up to that point in the simulation, then print a mes¬ 
sage in the output file exclaiming the need to raise ITL5 to continue a 
longer simulation. The default value of ITL5 is 5000 iterations. 

Depending on the circuit, 5000 iterations may require only a few 
minutes of computer time. Many simulations require more than 5000 
iterations during the transient simulation. Users have the choice of 
raising ITL5 to the approximate number of iterations required; simply 
setting ITL5 to 0 suppresses the limit on the total number of allowed 
iterations. 

.OPTIONS ITL6=X. The parameter ITL6=X is both a flag and a param¬ 
eter value. The default value of ITL6 is 0. If ITL6 is set to a nonzero 
value, SPICE discards the normal dc operating point calculation and 
replaces it with the source-stepping algorithm. The value declared for 
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ITL6 serves as the iteration limit for each step in the source-stepping 
algorithm. Many people wrongly believed ITL6 was related to the size 
or number of steps in the process; not so. The size of the step is fixed 
and cannot be changed. ITL6 only determines the number of iterations 
allowed at each step. 

The best use for ITL6 is in assisting in the dc operating point solu¬ 
tion for circuits where the .NODESET statement cannot be used. 
Using .NODESETs is the preferred way to find the dc operating point, 
but .NODESETs cannot be applied to nodes within a subcircuit defini¬ 
tion. For many large circuits, bistable nodes (such as flip-flops) cannot 
be accessed with the .NODESET statement. In these cases, raising 
ITL6 to 400 often helps achieve convergence. 

Using the iteration limits wisely is one of the pieces of the puzzle to 
producing fast, accurate, convergent simulations. Another part of the 
puzzle is the transient analysis options. These are covered in the next 
section. 

Transient analysis options 

For most designers, transient analysis is the most often used analysis 
type in SPICE. Transient analysis is also the most complicated analy¬ 
sis mode and therefore the most error-prone type of analysis in SPICE. 
To produce accurate transient simulation, the user must understand 
the limitations of the transient analysis algorithms including numeric 
integration, timestep control, and the six SPICE options which are 
applicable only to transient analysis. 

.OPTIONS LVLTIM=X. The parameter LVLTIM=X defines which time- 
step control algorithm SPICE uses during the transient simulation. By 
default LVLTIM is set to 2, which calls the Local Truncation Error 
(LTE) timestep control algorithm. Although slightly faster, the LTE 
timestep control algorithm is more error-prone than the Iteration- 
Count method. Without applying a limit to the maximum internal 
timestep, the LTE method may yield poor results when simulating 
asynchronous, sinusoidal, or inductive circuits. For these circuits, set¬ 
ting LVLTIM to 1 selects the Iteration-Count (IC) timestep control 
algorithm. The IC timestep control algorithm automatically limits the 
size of the internal timestep and may be much more reliable on these 
types of circuits. For a more detailed discussion of LVLTIM=X and 
timestep control algorithms, see Chap. 5 in this text. 

.OPTIONS METHOD=‘yyyy\ The METHOD =‘YYYY’ option defines the 
numeric integration method used to calculate capacitor currents and 
inductor voltages. Valid settings are METHOD =‘GEAR’ or METHOD^ 
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‘TRAP’. The default is TRAP which, not surprisingly, selects the trape¬ 
zoidal numeric integration method. The trapezoidal numeric integra¬ 
tion method is a relatively fast, accurate method, but trapezoidal 
integration is not foolproof. The trapezoidal method suffers from an 
annoying tendency to oscillate around the correct solution, especially 
on switching circuits or long transient simulations. The Gear method 
integration does not oscillate and tends to be more stable (stays close 
to the exact solution) over long transient simulations, but at the 
expense of longer simulation run times. 

.OPTIONS MAXORD=X. The MAXORD=X parameter defines the maxi¬ 
mum order for Gear’s multiorder integration method. When a user 
specifies METHOD=GEAR, SPICE uses the Gear second-order method. 
Higher-order Gear methods are available, including Gear’s third-, 
fourth-, fifth-, and sixth-order methods. To utilize one of the higher- 
order methods, the option MAXORD must be set to the highest order 
which will be used. MAXORD=3 specifies the third-order Gear method, 
MAXORD=4 specifies the fourth-order Gear method, MAXORD=5 
specifies the fifth-order Gear method, and, finally, MAXORD=6 speci¬ 
fies the sixth-order Gear method. 

The MAXORD=X parameter has a default value of 2 and may be set 
between 2 and 6. 

It should be noted, in theory, the higher-order Gear methods should 
produce less error at each timepoint, thereby resulting in less total 
timepoints and a faster transient simulation; due to the extra overhead 
involved with the higher-order Gear methods, however, none of the 
higher-order methods is significantly faster than the Gear second- 
order method. 2 

For a more thorough discussion on numeric integration methods and 
which type of integration to use for each circuit type, readers are 
directed to Chap. 4 of this text. 

.OPTIONS MU=X. The MU=X parameter was never documented in the 
SPICE2G.6 User’s Guide. 1 Like several other options, MU serves as 
both a flag and a parameter value. The default value of MU is .5. 

To understand MU, the equations for trapezoidal integration and 
backward-Euler numeric integration should be compared. There is a 
striking similarity between the two methods. If one carefully writes 
the equations with a single, well-placed variable, the value of the 
variable may be used to change the equation from trapezoidal inte¬ 
gration to backward-Euler integration. In SPICE, MU is that vari¬ 
able. By setting MU to .5, the default setting, the equations reduce to 
the trapezoidal integration method. By setting MU to 0, the equations 
reduce to the backward-Euler integration method, and by setting 
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MU somewhere between 0 and .5, a blending of the two methods is 
used. Readers are directed to Chap. 4 for a more thorough discussion 
of integration methods and when to use backward-Euler or trape¬ 
zoidal integration. 

.OPTIONS TRTOL=X CHGTOL=X. The TRTOL=X and CHGTOL=X 
parameters are two of the more misunderstood parameters. Because 
the names end in TOL, much like RELTOL, ABSTOL, and VNTOL, 
many people believe these are error tolerances also. Unfortunately this 
is not the case. Both the TRTOL and CHGTOL options were installed 
during the development of the Local Truncation Error (LTE) timestep 
control algorithm, and both are related to the LTE algorithm only. 

During the development of the LTE algorithm, the predicted value of 
local truncation error was found to be seven times larger than the 
exact value of truncation error. 2 To compensate for this discrepancy, 
the parameter named TRTOL was added to the equation for truncation 
error, and, not surprisingly, TRTOL has a default value of 7. 

While TRTOL has a direct relationship on the overall stepsize used 
during transient analysis, increasing TRTOL increases the stepsize, 
while decreasing TRTOL decreases the stepsize. Because TRTOL was 
added to compensate for a poor estimate of truncation error, this 
parameter is best left at the default value of 7. 

At the same time, it was observed that on certain circuits the Local 
Truncation Error timestep control algorithm would “lock up” by pro¬ 
ducing a step size that was infinitesimally small. This small stepsize 
occurred whenever the charge on a capacitor or flux through an induc¬ 
tor was less than the error tolerances of the program (RELTOL, 
ABSTOL, and VNTOL). To prevent this from occurring, a parameter 
named CHGTOL was added to the equations as a lower limit on capac¬ 
itor charge or inductor flux. Whenever the circuit produces a maximum 
capacitor charge or inductor flux which is less than CHGTOL, SPICE 
uses the value of CHGTOL in the LTE equation to predict the next 
timestep. CHGTOL is simply a means of preventing the LTE timestep 
control method from failing. 

CHGTOL has a default value of le-14 coulombs. Lowering this value 
does not significantly increase accuracy and only increases the possi¬ 
bility of lengthening the total transient simulation time. For these rea¬ 
sons, CHGTOL is best left at the default value. 


Numerical pivoting options 

The last two numeric control options were also the last two options 
added to SPICE. One of the last major changes made to the program 
was the addition of a modified numerical pivoting algorithm. Pivoting 
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reorders the matrix to reduce the number of operations required to 
solve the matrix. 

.OPTIONS PIVTOL=X PIVREL=X. The PIVTOL=X and PIVREL=X 
parameters both relate to the numerical pivoting algorithm in SPICE. 
The PIVTOL parameter defines the smallest numeric value (matrix 
entry) which is considered an acceptable matrix entry. If an entry is 
less than PIVTOL, a numerical overflow condition (such as a divide by 
zero) could occur. To prevent this, SPICE checks each of the matrix 
entries to ensure no pivot entry is less than PIVTOL. The PIVREL=X 
parameter defines the ratio between the largest entry in a given col¬ 
umn of the conductance array and PIVTOL. 

Numerical pivoting is a largely experimental science. Optimum piv¬ 
oting is different for different circuits. The default value of PIVTOL is 
le-13, the default value for PIVREL is .001, and both were determined 
empirically based on 64-bit double-precision conductance values. For 
this reason, both PIVTOL and PIVREL should be left at their default 
value. 

Numeric control options summary 

To achieve fast, accurate, convergence simulations, SPICE users must 
understand and know how to set the 18 parameters which form the 
numeric control option group. Quality simulations require significant 
user input. Part of the input is supplied in the form of these 18 numeric 
control option parameters. Chapters 3, 4, and 5 of this text describe in 
detail how each parameter should be set for a given circuit. Study these 
chapters and become familiar with the numeric control options. Time 
invested in these 18 parameters will produce the greatest return in qual¬ 
ity simulation when compared with other aspects of circuit simulation. 

MOS Geometry Options 

The next four parameters form the MOS geometry options. All four are 
applicable to MOS transistors only. The geometry of an MOS transistor 
gate has a large impact on the I-V and C-V characteristics of the 
device. Because, on any integrated circuit, many transistors share 
common geometries, four options were added to SPICE to reduce the 
effort in describing the transistor geometries in the circuit netlist. 


.OPTIONS DEFL=X DEFW=X 

The first two of these are the DEFL=X and DEFW=X parameters. The 
DEFL option describes the default gate length SPICE uses in the cal- 
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culation of the transistor drain current. In many IC houses, the mini¬ 
mum gate length is defined by the process of fabrication, and design¬ 
ers modify the width of the gate to adjust the drive current of the 
transistor. If the DEFL=X parameter is set to the minimum allowed 
gate length, designers need only specify the transistor gate width on 
the element line. If a designer chooses to use a transistor with a gate 
length other than the default, simply adding the L=X description to 
the element line overrides the default value declared with DEFL. 

The gate width parameter DEFW=X defines the default width of the 
transistor gate. Like the DEFL=X parameter, DEFW=X may be overrid¬ 
den by specifying W=X on the element line. For both DEFL=X and 
DEFW=X, the units of length and width are in meters. Remember to add 
the U (micro) scale factor if you describe these quantities in microns. 
Both the DEFL and the DEFW parameters default to 100 microns. 

.OPTIONS DEFAD=X DEFAS=X 

Just as the gate length and width may be described by default settings, 
so may the area of the drain and source be assigned default values. The 
area of drain and source are used in the capacitance calculation for the 
drain and source capacitance if the CJ model parameter is declared. 
The DEFAD=X defines the default value of the drain diffusion area and 
DEFAS=X defines the default value of the source diffusion area. Like 
the length and width parameters, both the drain area and source area 
are in units of meters squared. 

Miscellaneous Options 

The last four options fall into a category called the miscellaneous 
options. 

.OPTIONS TNOM=X 

The first of these is the TNOM=X parameter. TNOM defines the nomi¬ 
nal temperature of analysis. In most cases this is room temperature, or 
the nominal junction temperature. Values for TNOM should be 
expressed in units of degrees Celsius. By default, TNOM is set to 27C. 
While this is warm for room temperature, the assumption is that, even 
at room temperature, a semiconductor junction is generating some 
heat. TNOM represents the temperature of the device junctions. 

All of the temperature-dependent terms in SPICE change according 
to a change in temperature away from TNOM. TNOM should be the 
temperature used as a basis for temperature-dependent model param¬ 
eter or element value extraction. 
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.OPTIONS CPTIME=X LIMTIM=X 

The CPTIME=X and LIMTIM=X parameters are two more options 
which were installed when SPICE was originally developed and were 
added to limit the amount of computer resources a given simulation job 
could use. 

The CPTIME=X parameter defines the maximum CPU time which 
could be used to simulate a circuit. Because of the continued reduc¬ 
tion in cost of computer resources, most versions of SPICE have dis¬ 
abled this option by setting the default value extremely high. By 
default, CPTIME is set to 1 billion seconds. (For those of you who are 
reaching for your calculators to see how many days that translates 
into, let me save you some time. One billion seconds is approximately 
32 years!) 

Like CPTIME, LIMTIM is another safety net added to save com¬ 
puter resources. The LIMTIM=X parameter defines the maximum 
time allowed to generate the output. This parameter has been dis¬ 
abled in many versions of SPICE. By default, LIMTIM is set to 
2 seconds. 

.OPTION LVLCOD=X 

The last option is named LVLCOD. The LVLCOD=X parameter was 
used in early CDC versions of SPICE. The default for LVLCOD is 0. If 
LVLCOD was set to 1, early versions of SPICE called a subroutine 
named CODGEN. 

During the iterative procedure, with each new iteration, SPICE 
must solve the system equations array (the computer equivalent of 
performing a Gaussian elimination on a linear set of equations). 
SPICE solves the arrays using several nested fortran do-loops. 
Because fortran do-loops tend to be relatively inefficient on large, 
sparse matrices, the authors of SPICE decided to add a routine which 
would bypass the fortran do-loops and replace the solution with direct 
machine instructions. The routine which did this was named 
CODGEN. CODGEN wrote machine instruction to solve the system 
equations directly. While the addition of CODGEN reduced simulation 
time by 10-30 percent, the routine was dropped from later versions of 
SPICE. CODGEN was removed from SPICE for two reasons: first, 
because compiler efficiency and computer speed have both increased 
dramatically and, second, because the machine instructions were writ¬ 
ten specifically for a CDC computer. This made porting SPICE to an 
IBM or Cray or Vax or any other computer a nightmare. For both of 
these reasons, the CODGEN routine has been eliminated from most 
versions of SPICE. 
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Suggested Option Settings 

For most circuits, performance may be substantially improved by 
resetting several of the default option values. The following are sug¬ 
gested settings for all circuits. 

.OPTIONS ACCT LIST ITL1=500 ITL2=200 ITL4=40 ITL5=0 LIMPTS=10K 

Other options which need to be set include GMIN=X, RELTOL=X, 
ABSTOL=X, VNTOL=X, LVLTIM=X, and METHODS ABCD', but all of 
these must be set to values appropriate for your circuit. Chapters 3, 4, 
and 5 discuss setting specific options for different types of circuits. 

Summary 

The 34 options in SPICE define the information printed in the output 
file, the numeric algorithms used in the program, and iteration limits 
of the program. Learning to use and knowing how to set the proper 
options for a given circuit define the difference between a knowledge¬ 
able SPICE user and a novice, and the difference between fast, accu¬ 
rate, convergent, high-quality simulation results and garbage. SPICE 
is not as autonomous or reliable as many believe. SPICE simulations 
require significant user input to produce accurate results. The inputs 
SPICE requires include accurate netlist connectivity, accurate device 
model parameter sets and element values, and accurate option param¬ 
eter settings. 
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A Comparison of Vendor-Offered 
Simulator Options Settings 


Option 


Default Settings 



SPICE2G.6 (& Rspice) 

Hspice 

IS_Spice 

Micro-Cap IV 

Pspice 

LIMPTS 

201 

2001 

201 

n/a 

n/a 

NUMDGT 

4 

4 

4 

n/a 

4 

GMIN 

IpS 

IpS 

IpS 

IpS 

IpS 

RELTOL 

.001 

.001 

.001 

.001 

.001 

VNTOL 

luV 

50uV 

luV 

luV 

luV 

ABSTOL 

IpA 

InA 

IpA 

IpA 

IpA 

ITL1 

100 

200 

100 

100 

40 

ITL2 

50 

50 

50 

50 

20 

ITL3 

4 

3 

4 

4 

n/a 

ITL4 

10 

8 

10 

10 

10 

ITL5 

5000 

0 

5000 

1000000 

n/a 

ITL6 

0 

n/a 

0 

n/a 

n/a 

LTLTIM 

2 

1 

2 

n/a 

n/a 

METHOD 

TRAP 

TRAP 

TRAP 

n/a 

n/a 

MAXORD 

2 

2 

2 

n/a 

n/a 

MU 

.5 

.5 

.5 

n/a 

n/a 

TRTOL 

7 

7 

7 

7 

n/a 

CHGTOL 

.OlpC 

.OOlpC 

.OlpC 

.OlpC 

.OlpC 

PIVTOL 

.01p 

.001p 

.001p 

n/a 

n/a 

PIVREL 

.001 

.001 

.001 

n/a 

n/a 
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Gear integration, 88, 93, 96 
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On the disk which accompanies this text are the RSPICE 
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postprocessor which may be used to view and plot 
simulation results. Both of these programs use extended 
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text are free programs with a reserved copyright. 
Individuals may use, copy, and distribute these program 
free of charge. Users may wish to register their copies with 
RCG Research for technical support and upgrade notices as 
the programs evolve. 
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solutions for Chaps. 2, 3, 4, 5, and 6. Readers are 
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main commands required to run and execute any program in it. 

(search in Google for further information or tutorial about). 
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